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An Investigation of Laminar, Transitional, 
and Turbulent Heat Transfer on Blunt-Nosed 
Bodies in Hypersonic Flow! 


ROBERT J. CRESCI,* DONALD A. MacKENZIE,** PAUL A. LIBBY *** 
Polytechnic Institute of Brooklyn 


Summary c. = coefficient of specific heat at constant pressure 
H, = form factor 
Laminar, transitional, and turbulent heating rates have been h = enthalpy 
measured by means of the shrouded model technique. The h = h/h 
se 
Reynolds number was varied over a ninefold range; the enthalpy haw = he + 0.45 u,2, adiabatic wall enthalpy 
ratio (stagnation to wall) varied from 2.3 to approximately 1.5. K = 1—(1 — o!/3)((1 — &)/(1 — &)] 
Two different pressure distributions were imposed on the model k = coefficient of thermal conductivity 
which consisted of a spherically capped cone. M = Mach number 
The experimental data are compared to the laminar hypersonic Nyu = qulCp)eeRo/Re(hs, — hy), Nusselt number 
boundary-layer theory and shown to be in good agreement on the Ny.’ = local Nusselt number based on reference conditions 
conical portion of the model. On the spherical portion the data [defined by Eq. (13)] 
ss oe. me Of thts Ctscrepancy can be attributed to radiation to Nr’ = local Reynolds number based on reference conditions 
the nose of the model. 7 
[defined by Eq. (14)] 
The fully developed turbulent heat-transfer data are compared Ne = Neg,}/? 
to two theories: (1) a relatively simple turbulent theory which is Nuno = pte O/mm Reynolds number based on momentum 
based on recent theoretical work and which takes into account hitmen 
the upstream history of the boundary layer, and (2) the flat-plate ney 
n = coordinate normal to wall 
reference-enthalpy theory, which depends on only ‘‘local’’ con- p = pressure 
ditions. Although both theories are in reasonable agreement with p —s Ip 
the data, the latter method is simpler and somewhat more ac- vbcous : si 
Vw = heat transfer per unit area per unit time 
Ry = reference length (radius of spherical nose) 
For transitional flow the theory mentioned first can be readily aa ree : 
ro = radial distance from axis of symmetry 
modified in order to permit reasonable estimates of transitional i 
ae heat transfer to be obtained. On this basis it is possible to esti- a ; 
= distance along the surface measured from the stagna- 
mate laminar, transitional, and fully developed turbulent heat 2 : 
tion point 
= transfer under hypersonic blunt-body conditions. ; = ft 
ae The behavior of transition Reynolds number based on momen- ‘ ~ 
tum thickness is also discussed and shown to be in quantitative ve os eee 


agreement with recent shock-tube measurements. 


hs, 


t This research was conducted under the sponsorship of Con- 
tract No. AF 33(616-)6118 administered by the Aeronautical = Prandtl number 
Research Laboratory, WADC, USAF. w = wall shearing stress 

* Research Associate. Pre = Pate, 

** Research Assistant. 

*** Professor of Aeronautical Engineering and Assistant Subscripts 
Director of the Aerodynamics Laboratory. e = conditions external to the boundary layer 


u 
B = similarity parameter 
Symbols Y = ratio of specific heats 
= diu/d§, stagnation-point velocity gradient 
Cr = 27r,,/p,u,? local skin-friction coefficient 
Received June 12, 1959. p = mass density 
= Pe/ Pre 
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= stagnation conditions 
= wall conditions 


w = 
Superscript 
cy denotes conditions evaluated at the state correspond- 


ing to the reference enthalpy [defined by Eq. (12)] 


(1) Introduction 


T= DETERMINATION of large turbulent heat-transfer 
rates under conditions of strong favorable pressure 
gradients is required in several technological problems 
of current interest. Two examples of such problems 
are the turbulent heat transfer in the nose of hypersonic 
missiles and of aircraft in relatively low-altitude flight 
and in the throat region of hypersonic wind tunnels. 
Under these conditions the flow is characterized by large 
values of the ratio of stagnation to wall enthalpy, by 
relatively low local Mach numbers external to the 
boundary layer, and by a density at the wall much 
greater than in the stream external to the boundary 
layer. The present interest in this flow is indicated by 
the number of recent publications providing theoretical 
analyses thereof. References 1 through S have ap- 
peared since 1957. 

In all of these theoretical analyses a combination of 
assumptions must be made before a prediction of the 
heat-transfer rate is obtained. These assumptions per- 
tain, for example, to a relation between skin friction and 
a local boundary-layer thickness, and/or velocity and 
enthalpy profiles, and/or a relation between skin fric- 
tion and heat transfer. Unfortunately, the present state 
of knowledge concerning turbulent boundary layers 
does not permit a unique choice of assumptions; indeed, 
this ignorance, along with the large number of possible 
combinations of the several assumptions, accounts for 
the large number of publications presenting theoretical 
analyses thereof. It should be noted that some of the 
analyses in references | through 8 involve analytic com- 
plexity which is not warranted in view of the lack of 
substantiation of the basic assumptions; such com- 
plexity does not lead to increased confidence in accuracy 
of prediction. 

It is, furthermore, unfortunate that presently avail- 
able experimental data obtained under the requisite 
conditions do not supply the information necessary to 
establish the proper combination of assumptions. The 
data are not only limited in number but pertain ex- 
clusively to the measurement of heat transfer on blunt 
bodies under hypersonic flow conditions, or in regions of 
large enthalpy ratios but negligible pressure gradients. 
Since only the heat-transfer rate which is the result of 
a combination of assumptions can be subjected to ex- 
perimental verification, it is not possible to test the 
validity of each individual assumption. Consequently, 
simple theories with apparently little substance can be 
in better agreement with experimental data than more 
systematic theories. 

Consider the available experimental data. In ref- 
erence 9 Libby and Cresci presented laminar, transi- 
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tional, and turbulent heat-transfer data obtained by 
means of the shroud technique.'° The body tested was 
a portion of a sphere and was subjected to a pressure 
distribution corresponding to an expansion from the 
stagnation point to a local Mach number of approxi- 
mately 0.75. The ratio of stagnation enthalpy to wall 
enthalpy varied during the tests from roughly two to 
1.4 while the test Reynolds number was varied over a 
sixfold range. The experimental data were compared 
to the predictions of several hypersonic analyses; it was 
found that a simple analysis based on a flat-plate, in- 
compressible heat-transfer formula applied locally was 
in better agreement with experiment than the more 
rational analyses. 

These results were qualitatively confirmed by the 
data presented by Beckwith and Gallagher.'! Heat- 
transfer measurements were made on a sphere at two 
free-stream Mach numbers and with an enthalpy ratio 
(stagnation to wall) of about two. A twofold Reynolds 
number range was achieved. It was concluded that 
“. . the approximate level in the heat-transfer co- 
efficient as well as the variation with distance was fairly 
well predicted by a simple theory based on a turbulent 
heat-transfer formula for flat plates.” 

Recently, additional data involving transitional and 
turbulent flow have been provided in reference 12. 
The test conditions corresponded to enthalpy ratios 
from approximately 1.6 to 2. The Reynolds number 
was varied over a sevenfold range, transition being ob- 
served in the highest third of this range. The test re- 
sults for high Reynolds numbers were found to be in 
reasonable agreement with the predictions of a method 
given by M. R. Dennison, provided transition occurred 
“at a forward location on the body.’ For more rear- 
ward transition points the method is relatively in- 
accurate. 

Additional turbulent heat-transfer data with large 
values of the enthalpy ratio, relatively low Mach num- 
bers external to the boundary layer but in a region of 
essentially constant pressure are given in references 3 
and 13. These data were obtained by use of the shock- 
tube technique. Both reports concluded that a flat- 
plate, incompressible formula with the gas properties 
evaluated at conditions external to the boundary layer 
is in good agreement with the experimental data. 

A partial but incomplete basis for the agreement be- 
tween experiment and simple flat-plate theories for the 
pressure gradient case was provided by Vaglio-Laurin.* 
The transformation of the partial differential equations 
for the mean motion of an adiabatic, compressible tur- 
bulent boundary layer to incompressible form was 
demonstrated by Mager.'* Vaglio-Laurin showed, 
among other things, that by applying the same trans- 
formations and making the same assumptions concern- 
ing the turbulent shear stresses, the equations for the 
turbulent boundary layer with large enthalpy ratios and 
relatively low local Mach numbers could be reduced to 
an incompressible, constant pressure form. With this 
reduction as a basis, Vaglio-Laurin applied the logarith- 
mic velocity profile approach to obtain predictions of 
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Fic. 1. Test arrangement. 


turbulent heat-transfer rates which are in good agree- 
ment with the data of reference 9 and with those pre- 
sented below.* 

Concomitant with this reduction is the substantiation 
of the applicability of Reynolds analogy, to the case of 
simultaneous pressure gradient and heat transfer pro- 
vided p, p, ™ (u,/u)*, where p and uw are the mass density 
and velocity component parallel to the wall, respec- 
tively, and where the subscript e denotes conditions ex- 
ternal to the boundary layer. As was pointed out by 
Libby,” this approximation, which was used by Lees'® 
in developing a successful laminar hypersonic bound- 
ary-layer theory, corresponds to a boundary-layer form 
factor (/7,) of minus one. The approximation /7,~ —1 
not only provides the basis for the use of Reynolds 
analogy under pressure gradient conditions but elimi- 
nates in the more rational analyses the apparent dis- 
crepancy noted in reference 9 between the Reyn- 
olds analogy and what was termed therein the 
modified Reynolds analogy. In the latter analogy the 
Crocco relation between velocity and enthalpy was 
used in satisfying on the average the energy equation 
for turbulent flow while the Reynolds analogy assumes 
implicitly that the Crocco relation replaces the energy 
equation. However, the heat-transfer predictions based 
on the modified Reynolds analogy and form factors of 
either zero or unity were considerably lower than the 
experimental data of reference 9, while those based on 
Reynolds analogy and the same form factors were gen- 
erally higher than the data. If the form factor is minus 
one, the two analogies are identical and result in pre- 
dictions still somewhat high (10 to 20 per cent with 
respect to the data), but in considerably better agree- 
ment than those of previous, more rational methods 
considered in reference 9. 

It will be realized that the theoretical basis for the 
flat-plate formulas is still incomplete; if the inverse of 
the transformations discussed above were applied so as 
to obtain the heat transfer in terms of the physical 
variables, a functionally different form for the heat 
transfer is obtained from that given by a flat-plate 
formula applied locally. 

In this report, additional laminar, transitional, and 
heat-transfer data obtained by the shroud technique 
are presented. The model consisted of a spherically 
capped cone; the local Mach number external to the 
boundary layer varied from zero at the stagnation point 


* The data were made available prior to publication. 


to roughly three at the downstream end of the cone. 
The Reynolds number was varied over a ninefold range 
while the enthalpy ratio was varied from approximately 
2.2to 1.5. The laminar data are compared to the hyper- 
sonic theory of Lees'® with the stagnation point value 
modified according to Fay, Riddell, and Kemp."” The 
transitional and fully developed turbulent data are 
shown to be in good agreement with an analysis which 
is based on reference 4 but which is more convenient for 
the treatment of transitional heat transfer than that 
based on the logarithmic velocity law used therein. 

It should be noted that for many practical problems 
the ability to predict the heat-transfer rates in the 
transition region is of considerable importance since the 
transition region can be of significant spatial extent on 
the surface and can occur on a re-entering body for 
significant time periods in the trajectory. Analyses pre- 
sented in references 1-8 deal only with fully developed 
flow and do not permit treatment of the transition re- 
gion. In reference 9 the suggestion of Persh'* was ap- 
plied with both the Reynolds analogy and the modified 
Reynolds analogy and was shown to yield results in 
rough agreement with experiment. The analysis pre- 
sented- here similarly permits the prediction of transi- 
tional heat-transfer results. 

The authors are pleased to acknowledge the helpful 
discussions pertaining to the comparison of theory and 
experiment with Prof. Roberto Vaglio-Laurin, and the 
invaluable assistance of Samuel Lederman in the test- 
ing and data reduction. The authors apologize for not 
subjecting all of the available theories, for example, 
those in references 1—S, to comparison with the data pre- 
sented here; this is left to the individual authors. 


(2) Model Design, Testing Procedures, and Test 
Conditions 


The experimental results presented here were ob- 
tained by application of the shrouded model technique 
described in reference 10. The shroud was used in con- 
junction with the hypersonic facility of the Polytechnic 
Institute of Brooklyn described in references 19 and 20. 
The test arrangement is shown in Fig. 1. 

A schematic drawing of the model tested is shown in 
Fig. 2 with the thermocouple and pressure tap locations 


_-— PRESSURE TAPS 


. 
THERMOCOUPLES : 
#1-13 5 775 
THERMOCOUPLES 
# 30-130 i 
2_{3,30 5 |60| 7 | 8a} 12 (13,1 5,130 
301601 91 


THERMOCOUPLE # 7_| 8a 9 
SURFACE COORD.-5 | 0.800.85093 
9 


PRESSURE. T. 2/3 $i/6/7/86 10 | 12 13 
SURFACE COORD. #0080 1201710260 370 540 650 SBT OTT 541.621 1.88 213 


Fic. 2. Model design including thermocouple and pressure tap 
locations. 
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Fic. 3. Pressure distributions. 


indicated. The coordinate (5)°is-defined by ‘the dis- 
tance along the surface (s) measured from the axis and 
nondimensionalized with respect to the nose radius 
(Ro). The pressure taps and thermocouples were lo- 
cated along two mutually perpendicular diameters, 
more pressure taps being concentrated in the stagnation 
region in order to permit determination of the velocity 
distribution there more accurately. The model was 
made of type 304 stainless steel. 

The heat-transfer measurements were obtained 
by means of a transient, one-dimensional heat con- 
duction technique described in reference 21. The 
method consisted of recording the surface temperature 
of a one-dimensional plug, partially insulated from the 
surrounding model, as a function of time. The heat- 
transfer rates were then computed as a function of time 
by an EASE analog computer. The surface tempera- 
ture history was fed into the computer as one boundary 
condition; the inside surface was assumed to be in- 


4 sulated. The thermal properties of the metal were con- 
y sidered constant except in several check cases discussed 
below. 
s The pressures were measured by transducers and re- 
corded throughout the test on strip chart recorders. 
4 
20 

8 

= 
x 

02 LAMINAR: THEORY 4 
[€Q.(1)] 
0104 06 08 10 20 40 


Fic. 4(a). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple 41, 5 = 0 
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In the stagnation region, where the variation in pres- 
sure is small, differential transducers were connected 
between adjacent taps to measure pressure differences. 

In order to investigate the effect of a change in pres- 
sure gradient, the model was tested in two different 
axial positions with respect to the shroud. This re- 
sulted in the two pressure distributions shown in Fig. 3 
and denoted henceforth as positions Aand B. The pres- 
sure is nondimensionalized with respect to the stagna- 
tion pressure on the model and plotted as a function of 
(5). The data of Fig. 3 were obtained during several 
tests with the stagnation pressure held constant. Over 
the entire range of stagnation pressures the pressure 
distribution in terms of ) was found to be independent 


’ -of'thestagnation pressure. From Fig. 2 it will be noted 


that several pressure taps were placed out of the me- 
ridional plane in which the pressure distribution was 
measured. The pressure at these points was found to 
be axisymmetric within the measuring accuracy for 
pressure. No attempt was made in these tests to 
achieve a pressure distribution simulating actual flight 
conditions, the purpose of these tests was to investigate 
the accuracy of various analyses in the prediction of 
laminar, transitional, and turbulent heat-transfer rates. 
As a result, the measured distributions correspond to 
an acceleration in the stagnation region that is much 
lower than that predicted by Newtonian theory. On 
the conical portion of the model, on the other hand, the 
measured Mach number is larger than would be ob- 
tained in hypersonic free flight. 

The heat-transfer data were obtained by two different 
methods of testing; in one, the Reynolds number was 
held essentially constant during the test and the varia- 
tion of heat transfer with time was determined. In the 
other method the Reynolds number was varied con- 
tinuously during the run by changing the tunnel stagna- 
tion pressure. In this type of test, information concern- 
ing boundary-layer transition can be obtained. 
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Fic. 4(b). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple #2, 5 = 0.15. 
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The data are presented in terms of a Nusselt number 
defined by 
Nwu = — hw) 


and a Reynolds number parameter 


where Nr = Ro / se 


and Pse = Ps/ 


Although a continuous history of heat transfer and 
pressure were obtained, the final parameters were com- 
puted and plotted for discrete time intervals 4 seconds 
apart. 

In order to prevent the shroud from reaching exces- 
sive temperatures, internal cooling was provided by 
means of an annular jet of high-pressure, cold air 
located at the minimum section of the shroud. Monitor- 
ing thermocouples were inserted on the exposed surface 
of the shroud and the temperature history was recorded. 
During the tests it was found that the shroud tempera- 
ture was within +150°F. of the surface temperatures 
recorded on the conical portion of the model. 

In previous wind-tunnel tests using the shroud tech- 
nique (cf. reference 9), it was noted that the surface 
roughness of the model varied considerably during the 
tests due to air contamination. To reduce this effect, a 
double layer of fine mesh screens was installed in the 
tunnel approximately 4 diameters upstream of the 
model and shroud assembly, as shown in Fig. 1. These 
would be expected to provide a more uniform flow into 
the shroud, as well as a reduction in surface deteriora- 
tion. During this test series, the roughness of the model 
increased somewhat. Initially, the surface roughness 
was less than 5 microinches (r.m.s. as measured by a 
Brush Surface Indicator) over the entire model. The 
first eight tests were run with the model in position 
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Fic. 4(c). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple #3, § = 0.35. 
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“A,” the next five runs in position “B.’’ At the end of 
the test series, the roughness of the spherical portion of 
the model was measured to be between 40 and 50 micro- 
inches and on the conical skirt between 15 and 20 
microinches. No checks of the roughness were made 
between the various runs; as a result the variation in 
transition Reynolds number could not be correlated 
with the varying roughness. Since the data appear 
consistent within the experimental accuracy for equal 
Reynolds numbers in different tests, the roughness 
variation resulted in no observable effect on the transi- 
tion behavior. 

All tests were performed at a stagnation temperature 
in the range of 1,500°R. to 1,700°R. The model stag- 
nation pressure was varied from approximately 20 to 
180 psia. It is of interest to note that the highest 
Reynolds numbers in these tests corresponded to simu- 
lation of full-scale Reynolds number on a body with a 
nose radius of 2 ft. at a flight Mach number of 20 and 
at an altitude of 85,000 ft. Evidently the simulation 
was not complete, inasmuch as the flight values of the 
enthalpy ratio and of the absolute temperature were not 
simulated. 

No correction for radiation between the model and its 
surroundings was made. This appears to be clearly 
justified on the conical portion of the model since the 
temperatures of the shroud were close to those on the 
model. The spherical portion of the model, however, 
was subjected to radiation from the screens and nozzle 
in which the shroud was installed. The temperatures of 
these elements were higher than the model throughout 
the test runs, and as a result an incremental heat trans- 
fer due to radiation was imposed on the nose region of 
the model. It was not possible, however, to make more 
than rough estimates of this increment because of the 
complex geometry and the lack of detailed temperature- 
time histories throughout the region influencing the 
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Fic. 4(d). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple #3a, § = 0.35. 
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Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple #4a, 5 = 0.56. 


Fic. 4(e). 


nose of the model. The increment appears to be 5 to 15 
per cent of the total heat transfer. 


The surface temperature measurements are believed 
accurate to within +2 per cent. Stagnation tempera- 
ture was measured with an electrically heated probe 
which had a maximum error of approximately 1 per 
cent. These limitations include recorder accuracy, 
reading error, and the possibility of the thermocouple 
junction being slightly below the surface. The overall 
heat-transfer data are believed to be accurate within 
+15 per cent on the skirt of the model while on the 
spherical nose region this figure can be somewhat higher 
(+20 per cent). These figures are based on the self- 
consistency of the data obtained and include errors 
caused by the computation of heat transfer from the 
temperature histories. 
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Fic. 4(g). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple 46a, § = 0.70. 


JOURNAL OF THE AERO/SPACE SCIENCES—JUNE, 1960 


FULLY DEVELOPED —. VA 
TURBULENT ANALYSIS |, 
* 
1.0 
08 
04 \ 
| 
02 LAMINAR THEORY — 
(EQ.(1)] 
04 O6 O8 10 20 40 
Fic. 4(f). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple #5, § = 0.64. 


(3) Heat-Transfer Results 


The experimental heat-transfer data have been pre- 
sented in terms of the Nusselt number distribution with 
Reynolds number at a given station on the model in 
Figs. 4(a)-4(p). There the data from various test runs 
are denoted by different symbols. The open and closed 
symbols correspond to the data for pressure distribu- 
tions A and B, respectively. 

The data presented in Figs. 4(a)—4(p) were obtained 
during the variable Reynolds number tests. The data 
for fixed Reynolds numbers were in good agreement 
with the data presented herein and are therefore not 
given. This result indicates that within the experi- 
mental accuracy there is no effect of enthalpy ratio 
over the range extant in these tests when the heat 
transfer is presented in terms of the Nusselt and 
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Fic. 4(h). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple #7, § = 0.80. 
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Fic. 4(i). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple #8a, § = 0.85. 


Reynolds numbers used here. The insensitivity of 
these similarity numbers to enthalpy ratio provides the 
main justification for their selection; this point has been 
discussed extensively in reference 22. 

Inspection of Figs. 4(a)—4(p) indicates the typical 
heat-transfer behavior with Reynolds number at a given 
station. At low Reynolds numbers Vy, ~ Nr’, cor- 
responding to laminar flow; as the Reynolds number in- 
creases, transition occurs, and the Nusselt number in- 
creases sharply with Reynolds number; as the Reynolds 
number increases further, Ny, ~ N,»’’°, indicative of 
fully developed turbulent flow. 

In terms of this behavior it is possible to estimate the 
transition Reynolds number based on momentum thick- 
ness (denoted Nes) at each station for which transition 
occurred. The momentum thickness Reynolds num- 
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Fig. 4(k). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple #10, 5 = 1.06. 
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Fic. 4(j). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple #9, § = 0.93. 


ber was computed according to the theory of Lees'*® as 
discussed below. In Figs. 5(a) and 5(b) the experi- 
mental values of (Ny), ate compared to the variation 
of Neg with § for constant values of Nr. From this 
comparison the point on the body where transition was 
initiated can be observed. For example, in Fig. 5(a) 
which corresponds to pressure distribution A, transition 
appears to have been initiated first near the sphere-cone 
junction at Nz ~ 106 and (Ney) & 380 and to have 
moved forward onto the spherical portion as Nx in- 
creases. From Fig. 5(b) a somewhat more progressive 
movement forward of the transition point with increas- 
ing Reynolds number may be observed. It should be 
noted that some asymmetry in transition behavior can 
be noted from a comparison of the transition point 
location and the thermocouple locations given in Fig. 2. 
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Fic. 4(1). Distribution of Nusselt number with Reynolds num- 


ber. Thermocouple #11, § = 1.30. 
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Fic. 4(m). Distribution of Nusselt number with Reynolds num- 


ber. Thermocouple #lla, § = 1.30. 


It is of interest to compare the transition results pre- 
sented in Figs. 5(a) and 5(b) with the results of a shock 
tube investigation of transition presented in reference 
23. In the shock tube the enthalpy ratio (stagnation to 
wall) was considerably higher than in the present tests, 
being roughly 10 to 30. The model tested was a hemi- 
sphere-cylinder. The results indicated that the transi- 
tion point moved progressively forward along the cylin- 
der and occurred at a constant Nr» ~ 500-600 as the 
stream Reynolds number increased. However, when 
transition occurred at the hemisphere-cylinder junction 
a small increment in flow Reynolds number resulted in 
the appearance of transition in the sonic region of the 
body at Nry ™ 200-300. This corresponds to a 
“flashing’’ forward of transition point. The results pre- 
sented here are quantitatively in good agreement with 
the results of reference 23. The transition values of 
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Fic. 4(0). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple 413, § = 1.91. 
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Fic. 4(n). Distribution of Nusselt number with Reynolds num- 
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er. Thermocouple 412, § = 1.60 


Nrg are in the range 200 to 500 and are lower on the 
spherical portion of the model. Indeed there is some 
similarity in transition behavior with increasing 
Reynolds number inasmuch as transition does not move 
progressively forward. 

A more complete discussion of the heat-transfer re- 
sults will be given in Section (5) where there is a com- 
parison between theoretical analyses and these data. 


(4) Theoretical Analyses 


For the comparison of laminar heat-transfer data with 
theoretical predictions, the theory of Lees,'® modified 
for the more accurate stagnation point theory of Fay, 
Riddell, and Kemp,” has been used. In terms of the 
Nusselt and Reynolds numbers used here, there results 
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Fic. 4(p). Distribution of Nusselt number with Reynolds num- 
ber. Thermocouple 413a, § = 1.91. 
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Fic. 9(a). Distribution of Nusselt number along surface; pres- 
sure distribution A, Ve = 2.0 X 10%. 
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where the Prandtl number has been assumed to be 
0.71. 

Corresponding to Eq. (1), the Reynolds number 
based on momentum thickness is given by the theory of 


Lees as 
3 1/2 
Ne» | f | 


R 


where = / ue} 
f (pu/ peu.) [1 — (u/te) |dn 
0 


the momentum thickness, and » = coordinate normal 
to the wall. Eq. (2) is useful for the evaluation of the 
transition behavior of the boundary layer. 

In evaluating the right-hand side of Eqs. (1) and (2) 
from the known geometry and pressure distribution, the 
flow tables of reference 24 can be conveniently applied 
to a reasonable approximation.* For the test condi- 
tions in the tests reported here, an average isentropic 
exponent of 1.32 was used. It will be observed that the 
factor ¢,,!/* = (p./pseMts,)'/* appearing explicitly in Eqs. 
(1) and (2) will disappear when an isentropic flow rela- 
tion between velocity and pressure is assumed, there 
then results, for example, 


NyuNr~ (Prbtw/ 


It should be kept in mind that Eqs. (1) and (2) are 
based on the “‘cold wall’ approximation [p,/p (u/u,)*] 
which does not apply on the conical portion of the 
model, there the density at the wall is roughly equal 
to that in the external flow. If, however, the similarity 
parameter defined in terms of the variables of Eq. (1) 
as 


‘ 
paar)? 


* The right-hand side of Eq. (1) is relatively insensitive to the 
value of the isentropic exponent in the range 1.2 < y < 1.4. 


is on the order of 1 2 for all values of 5, then the analysis 
which results in Eqs. (1) and (22) is approximately 
valid. It will therefore be of interest to compute 6 
= B(5). It is to be noted that if the external flow be- 
comes uniform (wu, ~ constant), 8 = 0 and Eqs. (1) 
and (2) give approximate relations for the heat-transfer 
and momentum thickness distributions in constant 
pressure flows. Thus, the theory of Lees is frequently 
used to estimate the heat transfer on blunted slender 
cones and results in good agreement with experimental 
results (cf. reference 25, for example). 

For the comparison of theory and experiment in 
transitional and fully developed turbulent flows, a 
procedure similar to that denoted as Method VI in 
reference 9 has been followed. However, the analysis of 
reference 4 has indicated that a modification of the skin- 
friction law used in reference 9 should be made. This 
modification along with the assumption of a form factor 
Il, ~ -—1 discussed above leads to the following 
analysis: In anticipation of its use in transitional 
flows, the momentum equation for the axisymmetric 
boundary layer can be written in the form 


dNprs/ds = (Nr/ gue!) (pits = 
Nap (3) 


where c; = 27,,/p,t,”, the local skin-friction coefficient. 
Now the analysis of reference 4 suggests for fully de- 
veloped flows that the extension of the Prandtl incom- 
pressible skin-friction law to the compressible ‘“‘cold 
wall” case under consideration is 


c,/2 = 00.013) Neg 4) 


Egs. (3) and (4) with suitable initial conditions, usually 
taken to be at the stagnation point for an analysis of 
fully developed turbulent flow, give the skin-friction 
and momentum thickness distributions. Again on the 
basis of reference 4, Reynolds analogy can be applied 
to the pressure gradient case when the ‘‘cold wall’’ ap- 
proximation is valid. If the Reynolds analogy is modi- 
fied to account empirically for nonunity Prandtl number 
effects, there is obtained 


= — Nw) pelte(C4/2) (5) 


where o = Prandtl number and h,, is the adiabatic wall 
enthalpy defined by the equation 
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Fic. 9(b). Distribution of Nusselt number along surface; pres- 
sure distribution A, Nr = 1.7 X 108. 
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An analysis based on Eqs. (3) through (6) and ap- 
plied to a region of constant cylindrical radius (7) and 
constant properties external to the boundary layer is 
consistent with the experimental results of references 
3 and 13. For this case there results 

For regions of relatively low local Mach number /,, 
~ and to a close approximation the 
experimental data will agree as well with Eq. (7) as 
with the similar relation involving u,'/° explicitly. 

It should be noted, however, that the skin-friction 
law of Eq. (4) differs from that usually assumed (cf. 
Methods III and IV of reference 9 and reference 3, for 
example) by the factor 7. Therefore, the relation be- 
tween the Nusselt and Reynolds numbers is altered; for 
the fully turbulent case Eqs. (3) and (4) can be inte- 
grated in closed form with Ney = 0 at 5 = 0 and used 
in Eq. (5) to give 

Nyy 


ds 
0 


where K = (haw — Iw)/(se — Nw) 


(8) 


For transitional flow a numerical integration of Eq. 
(3) along with a skin-friction law equivalent to Eq. (4) 
is required. In accordance with the suggestion of 
Persh'* and including the modification suggested by 
reference 4, the skin friction is taken to be 


c,/2 = — B (9) 


where B is selected such that c,/2 is continuous at 
transition from laminar to turbulent flow. If Eqs. (3) 
and (9) are combined and Eq. (3) integrated from the 
transition point, then the skin-friction and momentum 
thickness distributions are known. Eq. (5) can then be 
applied to determine the heat-transfer distribution. For 
this purpose the following relation derivable from Eq. 
(5) and from the definitions of Ny, and Ne is con- 
venient: 
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Fic. 9(c). Distribution of Nusselt number along surface; pres- 
sure distribution B, Ng = 3.0 X 108. 
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Fic. 9(d). Distribution of Nusselt number along surface; pres- 
sure distribution B, Ng = 2.0 X 108. 


Since this analysis is based on the ‘cold wall’ ap- 
proximation leading to a modified skin-friction law, its 
applicability to the conical portion of the model tested 
in this investigation may be questioned. However, as 
in the case of laminar flow cited above, the application 
of the analysis to the case of uniform external flow (u, ~ 
constant) leads to an estimate of the boundary-layer 
behavior in constant pressure flows. The pertinent 
fluid properties have been evaluated at conditions ex- 
ternal to the boundary layer in accordance with the ex- 
perimental results of references 3 and 13. Therefore 
the estimates obtained from the above analysis when 
the local Mach number is high will not be in accord with 
experiments. However, this investigation is not di- 
rected to this case, which is of practical interest only 
for pointed bodies or for slender bodies with relatively 
small nose radii. 

Finally, as a matter of interest, the fully developed 
turbulent heat transfer has been compared to the pre- 
diction of flat-plate reference-enthalpy method (ef. 
Method I of reference 9), denoted henceforth as the 
FPRE method. For completeness the pertinent equa- 
tions for this method are repeated here; in terms of the 
Nusselt and Reynolds numbers used here, with a con- 
stant Prandtl number and with the recovery effect in- 
troduced by the factor A, there is obtained 


 0.0300'/3 
AY 


V 4/5 2/5 
4VR 


where p’ = p’/p,,, and @’ = w’/u,, and where the primed 
quantities are evaluated at the state corresponding to 
the reference enthalpy defined by* 


h’ = 0.5h, + 0.220'/*h,, + h.(0.50 — 0.220'/*) (12) 


It will be of interest to consider a local Nusselt num- 
ber - Reynolds number relation which forms the basis 
of Eq. (11). Define 


Nyu!’ = — (13) 
Nr’ = (14) 


Then the basic formula for the FPRE method with con- 
stant Prandtl number is 


*In reference 9 the Prandtl number was set equal to unity in 
defining h’. 
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Fic. 10. Distribution of Nusselt number along surface (data 
from reference 9). 


= (15) 


(5) Comparison Between Theory and Experiment 


The comparison between the theoretical analysis of 
Section (4) above and the data presented here has been 
carried out for laminar and fully developed turbulent 
flow and for four cases of transitional flow. 

Consider first the agreement between the laminar re- 
sults and the predictions of Eq. (1). In evaluating the 
right-hand side of Eq. (1), average stagnation and wall 
temperatures, and a constant isentropic exponent of 
1.32 were assumed. The variation of 7 was determined 
from the local temperature external to the boundary 
layer corresponding to the average stagnation tempera- 
ture and from the tabulated viscosity values of reference 
26. The velocity gradient C was determined from the 
accurate pressure measurements in the stagnation re- 
gion and from the assumption p ~ 1 —C?5*/2¢,, for 
i. 

The predictions based on Eq. (1) have been shown in 
Figs. 4(a)—4(p) and in Fig. 6. Also shown in the latter 
are the bands of data in terms of Ny, Np~"” at a 
given value of 5. Examination of these figures indicates: 
(1) the theoretical and experimental values are in good 
agreement over the conical portion of the model for 
both pressure distributions, and (2) on the spherical 
portion the experimental results are from 5 to 30 
per cent higher than the theoretical results with the 
mean experimental values from 20 to 25 per cent higher 
than theoretical results. On the assumption that the 
theoretical results were correct the possible sources of a 
consistent experimental error were investigated at some 
length; these are discussed briefly below. 


(a) Radiation From Screens and Nozzle 


As mentioned in Section (2) an increment in heat 
transfer to the nose region of the model due to radiation 
from the screens and nozzle walls upstream of the 
shroud is known to be present. This increment was 
estimated to be between 5 and 15 per cent of the con- 
vective heat transfer depending on the geometry, the 
temperatures, and the emissivities of the radiating ele- 


ments. It is possible that this accounts for a major 
portion of the deviation between the theoretical and 
experimental results. This is indicated also by the good 
agreement between the two sets of results obtained over 
the conical portion of the model where the radiation was 
negligible. However, the estimated increment was not 
considered sufficiently accurate to justify correction 
of the data. 


(b) Variation of Thermal Properties With Temperature 


The thermal conductivity (%,,) and heat capacity (c) 
of the stainless steel were known to vary by roughly 20 
per cent over the temperature range 500 to 1,200°R. 
It was therefore considered possible that some error in 
the data reduction from surface temperatures to heat 
transfer could be associated with the variation of ther- 
mal properties. This possibility was examined by 
noting that the thermal diffusivity was relatively in- 
sensitive to temperature over the cited range. Thus 
the consideration of variable conductivity and heat 
capacity could be reduced to one of constant properties 
by the transformation 


7 * 
0 


where # is the new dependent variable. The measured 
surface temperature history was transformed into #(/) 
and the heat transfer corresponding to variable k,, and 
c determined. It was found for several typical cases 
that this analysis resulted in a change in experimental 
values of several per cent—that is, within the accuracy 
of the analog computer used for data reduction. There- 
fore, the variability of material properties could not 
account for a significant deviation. Thus, all of the 
data presented here were obtained by solutions of the 
heat conduction equation involving constant properties. 


(c) Heat Conduction Along the Surface 


The analysis of the heat conduction in terms of a one- 
dimensional, unsteady system assumes that conduction 
along the surface into the heat-transfer plug at its small 
area of contact at each surface is negligibly small. 
Without a detailed two-coordinate, unsteady analysis 
it is difficult to assess this error. Although rough esti- 
mates of the possible effect of a small surface conduc- 
tion were made and resulted in changes of a few per 
cent, no convincing corrections were obtainable. 

As mentioned in Section (4), it is of interest to com- 
pute the parameter 6 as a function of §. The distribu- 
tions of 8 = (5) for the two pressure distributions are 
shown in Fig. 7 and are seen to vary from 0.5 to 
1.4; over this range the heat-transfer parameter 
Zw /(1 — gw) (ef. reference 16) varies by less than 10 per 
cent and therefore the heat-transfer predictions of Eq. 
(1) should be accurate over the entire model. 

Consider next the agreement between the fully de- 
veloped turbulent heat-transfer results and the pre- 
dictions of Eq. (8) and of the FPRE method. The 
same average temperatures, isentropic exponent, and 7 
variations were used in the turbulent analyses as were 
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used in the laminar case. The predictions of Eq. (8) 
have been shown in Figs. 4(a)—4(p) for the two pressure 
distributions. In addition, Figs. S(a) and S(b) give the 
distribution of Ny,,/N,’® over the model for the two 
pressure distributions; these figures, therefore, repre- 
sent the mean of the data corresponding to Vy, ~ Np’? 
at the values of § for which clear indications of fully de- 
veloped flow are obtained. Also shown in Figs. 8(a) 
and 8(b) are the predictions of Eq. (8) and of the FPRE 
method. 

Consideration of these figures indicates that on the 
conical portion of the model both theoretical predic- 
tions are in good agreement with the experimental re- 
sults, the FPRE method giving a somewhat better 
agreement. On the spherical nose the predictions of 
Eq. (8) and of the FPRE method are higher than the 
experimental data by approximately 20 and 7 per cent, 
respectively. On the basis of these data for fully de- 
veloped flow, it appears that for engineering calcula- 
tions the FPRE method applied locally is not only 
simple but quite accurate. 

Asa final comparison between theory and experiment, 
four cases involving transitional flow were analyzed by 
numerical integration of Eq. (3) with the skin-friction 
law of Eq. (9) assumed, and with the heat transfer com- 
puted from Eq. (10). The distributions of external 
flow properties and the average temperatures used pre- 
viously were again applied. The numerical integra- 
tions were carried out by the Kutta-Runge method with 
AS = 0.2. The Bendix G-15 computer at PIBAL was 
used. 

The results in the form of the Nusselt number dis- 
tribution over the body for a given value of the Rey- 
nolds number are shown in Figs. 9(a)—9(d)._ Examina- 
tion of these figures confirms the previous discussion 
with respect to the laminar and fully developed turbu- 
lent flow. In the transition region the agreement be- 
tween theoretical prediction and experimental result is 
qualitatively good and quantitatively fair. Estimates 
of the per cent error in this region are meaningless be- 
cause of the severe spatial gradients of heat transfer 
and the assymmetry of transition which is shown in 
these figures. It should also be noted that the heat 
transfer resulting from the transition analyses in some 
cases exceeds that given by the fully developed turbu- 
lent flow analysis. This discrepancy is associated with 
inaccuracy in the integration procedure and could be 
eliminated if smaller values of AS were used. 

To complete the comparison of the prediction of 
fully developed turbulent heat transfer given by Eq. (8) 
the data of reference 9 are presented in Fig. 10. Also 
shown is the result given by the FPRE method. In- 
spection of this figure indicates that the previous con- 
clusions with respect to these two predictions are sub- 
stantiated by this earlier data. It should be noted, 


however, that the results of Eq. (8) are in better agree- 
ment with the data than the methods which were given 
in reference 9 and which take into account the bound- 
ary-layer history. 

The relative accuracy of the FPRE method suggests 
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method with local Reynolds number. 


a presentation of data directly in terms of Eq. (15). 
A similar representation of data in terms of external 
flow properties and with the difference between the 
adiabatic wall and stagnation temperatures neglected 
was given in reference 22. Fig. 11 shows all the data for 
fully developed turbulent heat transfer presented here 
as well as from reference 9. The satisfactory agreement 
for a wide variety of conditions may be seen. This re- 
sult may be combined with that of reference 3 to indi- 
cate that local flat plate values satisfactorily correlate 
available data over a wide range of enthalpy ratios. 
However, it will be noted that the data of reference 3 
corresponds to constant pressure; thus additional data 
for large enthalpy ratios and significant pressure 
gradients are required before the general validity of this 
useful result is substantiated. 


(6) Conclusions 


The investigation discussed above has provided 
laminar, transitional, and turbulent heat-transfer data 
over a wide range of Reynolds numbers and with en- 
thalpy ratios on the order of two. The shrouded model 
technique has been employed. On the basis of a com- 
parison of the experimental results with theoretical pre- 
dictions, it can be concluded that the flat-plate ref- 
erence-enthalpy method provides a simple and yet 
relatively accurate means for estimating heating rates 
in fully developed turbulent flow. A rational method 
of permitting the inclusion of transition behavior and 
transitional heat-transfer rates has been suggested and 
is shown to be fairly accurate. The predictions for 
fully developed turbulent flow from this method are 
also reasonably accurate and relatively simple to ob- 
tain. 
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A Study of Spray Generated by Seaplane 
Hulls 


SHIZUO KIKUHARA* 
Shin Meiwa Industry Co., Ltd 


Summary 


The mechanism of spray generated by seaplane hulls is ex- 
plained by the same physical conception as in the slender body 
theory of aerodynamics The relation between two- and three- 
dimensional spray is made clear and demonstrated by experi- 
ments both in a towing tank and in a two-dimensional spray tank. 
The latter is essentially a U-tube of rectangular cross section and 
in it a two-dimensional vertical flow with free surface can be pro- 
duced in a glass passage. Two-dimensional spray shapes gen- 
erated by a wedge fixed in the glass passage can be observed 
photographically. An approximate theory of spray generation 
in two dimensions is presented for the case of straight wedges, 
which shows a fairly good agreement with the spray tank experi- 
ments. A groove-type spray suppressor was designed and tested 
in a towing tank. The effect of the new spray suppressor in- 
creases the load carrying capacity of a hull by about 40 to 50 per 
cent as far as the height of spray is concerned. 


Introduction 


a SoTTorF! showed by model experiments in a 
water tank in 1938 that the seaplane float having a 
high length-to-beam ratio L/B had excellent hydro- 
dynamic properties, a number of research efforts have 
been carried out in the U.S. to develop this configura- 
tion for seaplane hulls. The results of the research 
work have brought a remarkable progress in seaplane 
design. 

The obsolete principle of hull design in which the 
beam loading coefficient C,, = W/w B* was the control- 
ling factor was completely revised and a new coefficient 
k = W/w BL’ was proposed instead of C,, as the main 
design factor, where W = gross weight of seaplane, B 
= beam width, Z = length of hull, and w = specific 
weight of water. 

It was also shown by model experiments that the 
hydrodynamic properties could be kept unchanged pro- 
vided k was kept constant. Since Cy, = k(L/B)*, Ca, 
can be made larger than before by adopting high length- 
to-beam ratio L/B, thus greatly reducing the volume of 
seaplane hull. 

Although the hydrodynamic properties, such as the 
impact load, the stability, and the water resistance, 


Received April 28, 1959. 

* Senior Engineer 

The author wishes to express his sincere thanks to Prof. T. 
Moriya and Prof. M. Yamana of the Department of Aeronautics, 
University of Tokyo, and Prof. I. Tani of the Aeronautical Re- 
search Institute, University of Tokyo, for their helpful sugges- 
tions and encouragements throughout the preparation of the 
paper. Thanks are also extended to the members of the Hydro- 
dynamics Staff of Shin Meiwa Industry Co., Ltd., for their colla- 
boration in the experimental work. 


415 


have been clarified by the research work done within 
the last decade, the spray problems, especially the 
source and control of spray, have been left unsolved, so 
that the effect of spray on hulls still presents a problem 
involving considerable trouble. 

Since the size of hull for a given gross weight is us- 
ually determined so as to make the height of spray be- 
low a certain limit acceptable in practical operation, the 
volume, the structural weight, and the air resistance of a 
seaplane hull could be further reduced, provided the 
height of spray is lowered. 

The purpose of the present paper is to study the 
mechanism of spray generation and find some method of 
spray suppression. 

The study was carried out both experimentally and 
theoretically. The experiments were performed in an 
ordinary towing tank and also in a two-dimensional 
spray tank especially designed and built for the pur- 
pose. A newly developed groove-type spray suppres- 
sion device was tested in the towing tank and the effect 
was found to be remarkable. 

A theory of spray generation is proposed on the basis 
of the slender body theory in aerodynamics and is com- 
pared with the two- and three-dimensional spray test 
results. 


The Mechanism of Spray Generation 

Typical Form of Spray 

Fig. 1 is a typical shape of spray generated by a con- 
ventional hull bottom. There are two kinds of spray. 
The encircled A is called the velocity spray. It creeps 
upward with high velocity along the bottom surface 
and is shot out sideward relative to the moving hull. 
The velocity spray is usually considered harmless since 
it can easily be turned down to the water surface by a 
small amount of chine flare. The encircled B is the 
blister spray and is shot out from near the point of in- 
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Fic. 1. Typical shape of spray. 
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Fic. 2. Cy = 4.50 (near hump), k = 0.0376; Cy = V/VgB. 


Fic. 3. Cy = 1.80 (low speed), k = 0.0376. 
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Fic. 5. Comparison of the cross-sectional shape of blister spray 


at corresponding values of v and //v. 


TABLE 1 
Combinations of V and @ and Values of v 


6.5° | 8.6° | 11.5% | 15.5° 
3 m/sec. 0.45 

4 0.45 |0.60|,0.80 |,1|.08 
5.33 0.60} 0.80} |.08 
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Fic. 6. Two-dimensional spray tank. 
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Comparison of spray cross sections in two and three 
dimensions. 
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Fic. 8. Comparison of spray cross sections in two and three 
dimensions. 


417 


Fic. 9. 


3-Dimension 


—— 3-Dimension 0.8! 
2-Dimension O8| 


74 
74 


v d 
O8| 146 
08! 150 


Comparison of spray cross sections in two and three 
dimensions. 
am 
d 


64 
63 


1.08 
1.08 


2—Dimension 


1.07 
1.08 


146 
147 
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Spray in two dimensions (flat bottom), v = 0.81 m./sec. 
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tersection of the chine and the water surface and flies 
obliquely rearward relative to the moving hull in the 
form of a broad membrane or a cluster of water drop- 
lets. At relatively low speed the blister spray hits 
propellers or is sucked into the inlet duct of jet engine. 
As the speed increases the blister spray moves rear- 
ward, flies higher and longer, and hits the wing flap 
and the tail surfaces. The blister spray is the main 
object of the present study. 

The relative position of the two kinds of spray is 
clearly shown in Fig. 2, in which Cy = V/V gB, where 
g = acceleration of gravity. The velocity spray leaves 
the chine above the water surface. The forward end of 
the blister spray is near the point where the chine line 
cuts the piled up water surface; this point will be called 
the spray root. 

An example of the blister spray at comparatively low 
speed is shown in Fig. 3. The solid water separates 
from the hull side wall immediately behind the spray 
root and gradually comes back to it, leaving the blister 
spray some distance apart from the side wall. At a 
glance it seems as if the blister spray appears suddenly 
above the water surface, but this is not the case. In 
reality, it is most reasonable that the water particles of 
blister spray are shot out from the spray root with the 
local high velocity, fly in the air for a time, and then 
fall back on the water surface. 


Application of the Slender Body Theory to the 
Mechanism of Spray Generation 


According to the slender body theory in aerodynamics 
the flow around a slender body moving in a fluid with a 
small angle of attack may beconsidered two-dimensional 
in a plane perpendicular to the direction of motion. 
The two-dimensional flow is approximately similar to 
that caused by the movement of the body of that cross 
section in the direction perpendicular to the axis. 
The physical conception may be applied to explain 
the mechanism of the spray generation, since the change 
of the cross-sectional shape of the usual type of hull is 
slow longitudinally and also the trim angle is small. 

Now let us consider the case of a hull running on the 
water and observe the phenomena occurring in a fixed 
plane perpendicular to the direction of motion. The 
keel of the hull passing the plane penetrates vertically 
downward into the water surface at a ratev = V sin @, 
where @ is the angle between the tangent to the keel 
line and the still water surface and V is the velocity of 
the hull. At the moment of chine immersion—that is, 
at the moment when the spray root passes the plane— 
water particles near the chine in the plane are shot out 
into the air, because they are given high local velocity. 
(This point will be discussed more fully later.) The 
local high velocity at the spray root is caused by the 
rate of penetration v, the other component velocity 
V cos 6 parallel to the keel giving little change in ve- 
locity. Thus, v is an important factor of spray. 

If the cluster of water particles shot out is assumed to 
fly in the perpendicular plane, its shape and position in 
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that plane changes as time goes on. Meanwhile, a 
station B (Fig. 4) of the moving hull located at a dis- 
tance / behind the spray root reaches the plane at a 
time //V later than the spray root. Therefore, the 
shape and position of the cluster of water particles at 
that time contribute to the shape of the spray cross 
section at station B. Thus, the shape at a station de- 
pends on //V. 

In an ordinary towing tank the station at which the 
phenomenon stated above occurs moves forward with 
the velocity V. Consequently, it seems as if the hull 
moves with a stationary spray accompanied along its 
side. 

Now we can determine the relation between the two- 
and three-dimensional spray shapes. Consider a two- 
dimensional wedge of the cross section identical with 
that at the spray root set in a two-dimensional vertical 
flow with a free surface rising with a velocity v = V sin 
6. The shape of two-dimensional spray observed at a 
time ¢ after chine immersion corresponds to the shape of 
spray cross section occurring at a station / = Vt behind 
the spray root of a three-dimensional hull moving with 
V and 6. In other words, the three-dimensional spray 
shape can be obtained by locating the two-dimensional 
spray shapes at successive times ¢ at distances x = Vt 
behind the spray root. 


Measurement of Spray Generated by a 
Cylindrical Model 

To give an experimental proof of the postulate for 
spray generation described in the preceding section, a 
cylindrical hull model, having a length of 4.9 ft. and 
the same cross section as that of the station 500 of model 
I (see Fig. 28), was tested in a towing tank. According 
to the postulate, the cross-sectional shapes of the blister 
spray measured at stations having the same value of 
//V should be the same for various combinations of V 
and 6 provided the value of v is kept constant. Since 
the draft d is / tan 6 for the present model, //V = d/v 
for small values of 6. Therefore, the cross-sectional 
shape of the spray should be the same at stations having 
the same value of draft d. 

The combinations of V and @ in the test were selected 
as given in Table 1. 

The spray was measured at three stations for each 
run and the results are summarized in Fig. 5, in which 
the shapes of spray for corresponding sections are drawn 
on the same figure for the purpose of comparison. The 
coincidence is very good except for @ = 15.5°. Thus it 
may be said that the postulate concerning the mech- 
anism of spray generation is applicable for small values 
of #6. For values of 6 larger than 12° ~ 15° the assump- 
tion that water particles move in the perpendicular 
plane is no longer valid, their forward component ve- 
locities becoming too large to be neglected. 


Measurement of Two-Dimensional Spray 


Two-Dimensional Spray Tank 
It has been considered desirable to find some means 
to observe the two-dimensional spray shapes for a body 
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Fic. 12. Spray in two dimensions (deadrise 25°), v = 0.81 
m./sec. 
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spray photographs. (Top) Deadrise 
Fic. 13. Spray in two dimensions (deadrise 50°), v = 0.81 0°, v = 1.08 m./sec. (Center) Deadrise 25°, v = 0.71 m./sec. 
m./sec. (Bottom) Deadrise 50°, v = 0.91 m./sec. 
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of arbitrary cross section. For this purpose a new test 
apparatus was designed and built, which was called the 
two-dimensional spray tank. It is essentially a U-tube 
of rectangular cross section. Fig. 6 shows the general 
arrangement of the spray tank. 

By sending compressed air into the top space of the 
steel tank, water is pushed up into the glass passage so 
that a two-dimensional upward flow with free surface is 
obtained. The cross section of the glass passage is 2 in. 
by 20 in. and the height is 28in. By virtue of the high 
resistance of a perforated plate in the water passage, 
the velocity becomes constant in a short time. This 
velocity is controlled by changing air pressure. A half 
beam model is fixed on the side wall of the channel, and 
the shape of spray at any instant can be observed photo- 
graphically. Models having the same beam as the 


three-dimensional models were used throughout the 
present test. 


Comparison Between Two- and Three- 
Dimensional Sprays 


Two-dimensional spray photographs were taken in 
the spray tank with a model having the same cross 
section as that of the cylindrical model of the preceding 
section. For the same value of v, the draft d was used 
to determine the time ¢ = // V = d/ V, which corresponds 
to the three-dimensional case in Fig. 5. 

The test results are shown in Figs. 7 to 10, in which 
two-dimensional spray photographs and their sketches 
are shown together with the corresponding three-di- 
mensional shapes of spray cross section obtained by the 


3310 
U= 0.6 
v= 


, 
| 
. 
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Fic. 17. Spray curve (straight wedge deadrise 50°). 
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Fic. 18. Initial stage of spray generation. 
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Velocity and direction on wedge and free surface. 
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Surface shape, deadrise 25° and 50°. 


Fic. 20. Surface shape at chine immersion (deadrise 50°) 
in spray tank. 
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towing test in the preceding section. The agreement is 
quite good, and it may be concluded that the theory is 
applicable with sufficient accuracy for most practical 
purposes. 


Test Results of a Few Sections 


Two-dimensional spray photographs of a few sections 
taken in the spray tank are shown in Figs. 11 to 15. 

Figs. 11 to 14 are for deadrise angles of 0° (flat bot- 
tom), 25°, and 50°, Fig. 14 being a set multi-exposure 
photographs of the three models. Sharper bottoms 
generate lower sprays. 

The effect of chine form is shown in Fig. 15 for a 
deadrise of 25°. The shape of spray is not much al- 
tered by a slight chine flare (model S-25 + 3B), but it 
becomes very much scattered by a large flare (model 
S-25 + 6B), however, the height of spray is not lowered. 
When the chine is rounded (model S-25 — 9B) the 
lower part of the spray is drawn close to the side wall, 
suggesting one of the design principles of the groove 
type spray suppressor (see below). 


Spray Curve 
Some of the results of the spray tank tests are sum- 
marized in the form of the “spray curve’’ which shows 


_ Blister Spray 
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Fic. 22. Change of cross-sectional shape of two-dimensional 
spray with time. 
Cv-073 vou Experiment 
£ Theory 
Undisturbed W.L. at t=0.090Sec 
0.058 
0.026 
25° 
Fic. 23. Cross-sectional shape of blister spray (deadrise 25°). 
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the relation between the height of spray / and the draft 
d (see Figs. 16 and 17). The values of h and d were 
measured from the true still water surface corrected for 
the effect of the side wall of the spray tank. 

The height of three-dimensional spray can be pre- 
dicted by the spray curve and the foregoing theory, 
provided the planing conditions, V and @, are known. 

The spray tank has made it possible to study the 
spray problem more closely than before. 


An Approach to Two-Dimensional Spray 
Analysis 


The problem is to predict analytically the shape of 
spray generated by a straight wedge of finite width 
penetrating vertically with a constant velocity into a 
perfect fluid (water in the present case). 

Since the problem is of unsteady nature with free sur- 
face, it is too complex to be treated by a strict analytical 
method. (The spray tank was developed for this 
reason.) 

An attempt has therefore been made to solve the 
problem approximately with a view to understanding 
the mechanism of spray generation more thoroughly. 


Cv=0.73 © Experiment 
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Cross-sectional shape of blister spray (deadrise 50°). 


Fic. 25. Design principle of spray groove. 
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Fic. 26. Spray in two dimensions with groove-type spray suppressor. v = 0.81 m./sec. 
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Initial Stage of Spray Generation 


Since the spray breaks out of the free surface im- 
mediately after the chine immersion, its shape and 
height should be determined mainly by the conditions 
of the free surface at that moment. The shape and 
velocity distribution of the free surface before and at 
chine immersion can be determined theoretically by 
virtue of the law of similitude developed by H. Wagner* 
and J. D. Pierson.‘ Examples of calculated results for 
deadrise angles of 25° and 50°, are shown in Fig. 18 and 
Figs. 19(a) and 19(b). 

The thin water film along the bottom in Fig. 18 and 
Fig. 19(b) will be called the spray tip. The part of 
the water surface where the pile-up changes into spray 
tip describes a sharp curve. This part will be called 
the spray root in two dimensions, in which high local 
velocity is occurring as shown in Fig. 19(a). 

Among various forces acting on the wedge before 
chine immersion the unsteady state force is by far the 
largest and is equal to the time rate of change of mo- 
mentum in the fluid. The momentum is proportional 
to the rate of penetration V, of the wedge multipled by 
the virtual mass; that is, V, r/2 pC’, C being the wetted 
width. 

The work done by the force is accumulated in the 
water as kinetic energy, which increases as the depth of 
penetration increases. After the chine immersion, the 
unsteady state force becomes zero, since C becomes con- 
stant and the rate of penetration is also assumed con- 
stant. Therefore, the kinetic energy in the fluid reaches 
its maximum value at the instant of chine immersion, 
and will become the motive force to generate the spray, 
the other forces being assumed to be small. 

From the keel contact to chine immersion the bound- 
ary conditions which materialize the theory do not 
change and the free surface maintains the similar form. 
At the moment of chine immersion the boundary condi- 
tions change abruptly; there is no further broadening of 
the bottom surface which restricts the motion of the 
water particles, and, therefore, the free water surface 
begins to change its shape into the spray. 


Velocity Spray 

The spray tip will become the velocity spray in three 
dimensions. According to the preceding theory the 
shape of the spray tip should be triangular, but in the 
towing tank it is gathered into a cluster of irregular form 
by the surface tension, as shown in Fig. 20. As the 
velocity is constant in the spray tip it may be assumed 
that the velocity of the cluster is the same as that of the 
spray tip. 

The cluster flies along a path determined by the ini- 
tial velocity under the influence of gravity force and air 
resistance, as shown in the multi-exposure photographs 
in Fig. 14. Values of the initial velocity were deter- 
mined photographically from the spray tank experi- 
ments and compared with the theoretical velocities ob- 
tained by the law of similitude in Fig. 21. The theory 
give velocities about 20 per cent higher than the experi- 
mental values for the deadrise 8 = 25° and 50°. 
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Fic. 27. Spray curve in a two-dimensional spray tank. 
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Fic. 28. Lines of model I. 
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Fic. 29. Spray suppressor of model II. 
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Blister Spray 


The water sheet lying on the surface of the spray root 
at the moment of chine immersion will become the 
under surface of the blister spray in three dimensions. 
The present analysis attempts to trace the path of the 
water particles in the sheet mentioned above. 

Let U,, U, and @,,, 4,, be the initial velocities and di- 
rections of the particles at the points by) and Cy on the 
surface in Fig. 22. These particles were assumed to be 
independent of each other and fly from by and c along 
the path determined by U and @, under the influence of 
gravity, surface tension, and air resistance. Strictly 
speaking, they are not independent, and the influence of 
the adjacent fluid can not be disregarded. However, 
because they have very high local velocities and are 
shot out in the air immediately after chine immersion, 
the foregoing assumptions may be granted so far as the 
path of the local surface sheet of the spray root is con- 
cerned. 

Since U,> U,and 6,, < 6,. [Fig. 19(a)|, the sheet bc 
turns upside down soon after the particles are shot out 
and expands its area in the course of flight as the under 
surface of the blister spray. It is necessary to deter- 
mine the mass of the flying water sheet to go through 
the computation. 

Owing to the difficulties involved, however, a trial 
and error method was adopted: the shape and position 
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were computed assuming several values of the mass or 
the mean thickness of the water sheet chosen arbitrarily, 
and the one which gave the best overall agreement with 
the test results was chosen as the most probable value. 

The initial velocities of the particles were reduced to 
80 per cent of the theoretical value for the reason de- 
scribed previously about the velocity spray. 

The effect of surface tension was taken into account 
as the loss of kinetic energy due to the work done by the 
expanding sheet but was found negligible for both 8 = 
25° and 50°. 

The effect of air resistance on the position of the 
flying sheet with expanding area was calculated as the 
loss of the velocity as time goes by and was found to 
have a considerable effect for 8 = 25° but little effect for 
8B = 50°, because for the latter case the rate of expan- 
sion and the velocity of the sheet are small. 


The calculated shapes of the water sheet are shown in 
Figs. 23 and 24, together with those copied from the 
multi-exposure photographs. The mean thickness of 
the water sheet was assumed 0.02 in. at ¢ = 0.032 sec. 
for 8 = 25° in the present case. They resemble each 
other in shape and position to some extent. 

The analysis described above is based on too many 


assuniptions to be called a theory; nevertheless, it 
seems to be of value for understanding more thoroughly 
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Fic. 30. Size of hull and gross weight. 
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the mechanism of spray generation and the process of 
its subsequent development. 


A Groove-Type Spray Suppressor 

Design Principle of Spray Groove 

The mechanism of spray generation discussed in the 
preceding sections suggests that there will be no blister 
spray if we can dispose of the high-speed water particles 
in the spray root before they are shot out high to form a 
spray sheet 

The cross-sectional shape of a blister spray is in 
general as shown in Fig. 25(a). The high-speed water 
is separated from the sharp chine. If the chine is 
rounded as in Fig. 25(b) the water sticks on the surface 
and goes up a distance along the side wall and can easily 
be caught by a groove as in Fig. 25(c) or Fig. 25(d) be- 
fore being shot out in the air. The water will flow 
through the groove rearward without any spray or 
come out through a slit so low as to be harmless. The 
lid of the groove may be folded in the air to minimize the 
air resistance. 


Test Results in Two Dimensions 


Fig. 26 shows the two-dimensional spray photographs 
of two types of the spray groove (B-1 and B-2) and 


Fig. 27 shows the spray curve of the model B-1. The 
heights of the sprays of these models are much lower 
than those of the conventional models without spray 
groove shown in Fig. 12 and Fig. 16. 

Let the cross-sectional area of the spray groove be A 
that is just necessary to eliniinate the spray in two di- 
mensions. Then AV is the minimum quantity of 
water per unit time which should be disposed of to 
eliminate the blister spray in three dimensions. If the 
sectional area of the grooves is less than A some amount 
of spray will appear and the height of spray will be- 
come higher as the area of the groove becomes smaller. 
In actual design it is not necessary to eliminate the spray 
completely. It will suffice for practical purposes to 
suppress the height of spray so low as to be harmless for 
the specific design configuration. A two-dimensional 
spray tank will be a convenient apparatus to determine 
the necessary groove area for any shape of hull section. 


Test Results in Three Dimensions 


Two models, I and II, were prepared. The line 
drawings of the models are shown in Figs. 28 and 29. 
Model I has a standard bottom shape with L/B = 13.7, 
and model II has the same general lines as model I ex- 
cept that a groove-type spray suppressor is fitted along 
the forebody chine line. 
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Fic. 31. Cross-sectional shape of blister spray. 
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Fic. 32. Cross-sectional shape of blister spray. 


ue. 
to 
de- 
int 
the 
the 
the 
to 
for 
in- 
in 2 
he 
of 
ec. 
ch 
ny 
it 
ily 
4 
| 
N TAIT TTT stasoo 
LW | LITA 
1 
++ +--+ + 
stais0o |, , | | LAT 
ml 


428 JOURNAL OF THE AERO/SPACE SCIENCES—JUNE, 1960 


Sta. 300 600 708 1200 1691 
Mode! | Without Spray Suppressor. 


Envelope of Spray Height. 


Fic. 33. Comparison of spray projection and envelope. 


Both models were tested for a gross weight W up to 
14 kg., that is, k = W/wBL? = 0.0376, where B is the 
beam width, Z is the length of hull and w is the specific 
weight of the water. In usual practice the value of k 
= 0.025 is considered the limit of loading in order to 
hold the spray height below a certain value acceptable 
in practical operation. 

Fig. 30 shows the general appearance of sprays for k 
= 0.0376 at a speed near the hump. The cross- 
sectional shapes of the spray were measured during 
test runs by a tool of the form of a long toothed comb. 
Some of the test results are shown in Figs. 31 and 32 and 
summarized in Fig. 33. The effect of the spray groove 
is clearly seen in these figures. 

The envelope of the spray height of model II with 
spray groove for k = 0.0376 is lower than that of model 
I without a spray groove for k = 0.0268 (Fig. 33). The 
latter value of k being around the practical limit of 
loading, it may be said that the load-carrying capacity 
of a hull is increased by about 40 per cent by the spray 
groove so far as the spray height is concerned. 

Fig. 34 compares the above-mentioned values of k 
with those of actual seaplanes hitherto designed and 
manufactured. There are four points plotted which 
represent typical seaplanes during the period of 1940 to 
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Fic. 34. Spray in three dimensions. 


1955, together with the points k = 0.0268 for model I 
and k = 0.0376 for model II. 

The actual seaplanes were all designed around the 
value of k = 0.027, while for the model II with the 
spray groove the value of k = 0.0376 may be permissible. 
It is not too much to say that the idea of spray groove 
might provide a powerful means to mark a step toward 
the possibility of designing heavily loaded hull to reduce 
the size, structural weight, and air drag, thus resulting 
in a substantial improvement of the performance of 
water-based airplanes. 
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Supersonic Inlet Dynamics 


H. R. FRAISER* 
Chance Vought Aircraft, Incorporated 


Summary 


An approximation of the differential equation for compressible 
duct flow is presented. The equation is linear and of the second 
order. The duct transfer function and response characteristics 
are obtained by applying small-perturbation theory to the dif- 
ferential equation. The resulting equations describe duct 
natural frequency as a function of duct areas and volumes, and 
damping ratio as a function of the slope of the steady-state mass 
flow, pressure-recovery curve. 

The calculated response agrees, to a first approximation, with 
measured response as obtained from tests of a fixed-geometry, 
sugar-scoop inlet model with bypass for matching airflows. Test- 
ing was done in the 10 X 10 and 8 X 6 ft. supersonic tunnels 
at NASA Lewis Flight Propulsion Laboratory. Further agree- 
ment was obtained during flight tests of the F8U-3 airplane. 


Symbols 
Pi = total pressure, lbs. per sq-ft. 
p = static pressure, lbs. per sq ft. 
T: = total temperature, °R. 
‘i = static temperature, °R. 


= specific weight, lbs. per cu.ft. 

= Mach number 

= velocity, ft. per sec. 

= volume, cu.ft. 

= airflow rate, lbs. per sec. 

= acoustic velocity, ft. per sec. 

= acceleration, ft. per sec. per sec. 

= area, sq.ft. 

= length, ft. 

= dynamic head, lbs. per sq.ft. 

= gas constant, lbs. (mass) ft./lb. (force), °R. 
= gravitational constant, ft. per sec. per sec. 
= ratio of specific heats, Cp/Cy = 1.4 

= aconstant 

= natural frequency, rad. per sec. 

= natural frequency, cycles per sec. 

= temperature ratio, 7;/T 


| 


= P,/2117 
(fz)m = afunction of Mach number 
Ws = air mass in the control volume, Ibs. 
Subscripts 
c a BP e 

0 = free stream 
i = inlet 
c = control probe station 
a = by-pass station 
e = exit and/or choke point 
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bp = by-pass exit 

1 = section ahead of bypass 
2 = section aft of bypass 

s = capacitive or stored term 
j = inertance term 

R = resistive term 


Superscripts 

, = first derivative with respect to time 
= second derivative 

= an average value 


Introduction 


SS supersonic aircraft designed for 
Mach numbers greater than 1.5 must utilize some 
type of variable geometry in the engine air duct in order 
to match engine airflow to inlet airflow and obtain the 
optimum propulsive efficiency. This variable geometry 
usually includes a by-pass arrangement and, in more 
sophisticated designs, variable supersonic diffusers. 
Supersonic-diffuser geometry control may be obtained 
from an open-loop servo system which programs inlet 
geometry as a function of some airplane flight param- 
eter. By-pass control is more critical and must be 
obtained from some duct-airflow parameter such as duct 
Mach number or normal-shock position. An error of 
1 per cent in by-pass control can result in 3 to 5 per 
cent reductions in net propulsive force. Thus, the 
servo system gains must be as high as stability require- 
ments will allow. Furthermore, by-pass response 
must be very fast in order to follow engine transients. 
If the by-pass system fails to follow engine transients, 
inlet ‘‘buzz’’ may result, and/or the engine compressor 
may stall as the result of “buzz” or a deterioration in 
engine face pressure distribution. It is evident that, 
as part of the closed control loop, inlet duct dynamics 
must be known in order to design a high-performance 
control system. 

A block diagram for a typical by-pass control system 
is shown in Fig. 1. The free-stream transducer senses 
free-stream pressure ratio (po/P;,) and schedules duct 
pressure ratio (p,/P;.) as a function of free-stream 
pressure ratio by means of a one-dimensional, nonlinear 
cam. The duct transducer senses the existing duct 
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Fic. 1. Block diagram of a by-pass control system. 


pressure ratio. The two output signals are compared 
and the error signal drives the by-pass doors in the direc- 
tion to reduce the error. By-pass door movement feeds 
back through the duct to the duct transducer. It is 
the purpose of this paper to present a method for deter- 
mining the transfer function of the inlet duct which 
forms the feedback loop of the system just described. 


Analytical Evaluation of Duct Dynamics 


The inlet duct to be analyzed is a fixed-inlet duct 
having a variable bypass which is controlled by duct 
Mach number ahead of the by-pass station. The 
engine will be assumed to be operating at a constant 
corrected airflow (W.V 0./6,). The desired transfer 
function for such a system is AM,/AA,, where M, is 
duct Mach number at the control station and A,, is 
by-pass exit area. The relation for this transfer func- 
tion is obtained by using the ‘‘lumped constant’’! con- 
cept where it is assumed that the air mass within a 
control volume may be considered as a single homog- 
enous slug, the properties of which can be represented 
by values taken at some average point in the control 
volume. The continuity equation then becomes, for a 
transient condition, 


AM, 


W,.-W.=% (1) 


where W, = inlet airflow rate, W,. = airflow rate at the 
engine choke point, and w, = the first derivative with 
respect to time of the slug in the control volume. The 
pressure drop from inlet to exit is 


Pu — Py = APs, + (2) 


Where AP,, is the steady-state ‘‘pipe’’ loss and AP;, is 
the differential required to accelerate the slug in the 
control volume during a transient. Using the assump- 
tions listed in Appendix (A), relations (1) and (2) can 
be combined and expressed as a function of A/,. The 
resulting equation, as derived in Appendix (A), is 


KK 4, A, Ky, K 
K,14M, 


fin) ae 


+ K.,Kr(fy)ue — ( + 


M, + KA fu) ue — 
K.( fx) KopA fr) (3) 


Inspection of the terms containing 1/, shows that, for 
small perturbations about a fixed point, all Mach func- 
tions may be considered constant except (fu-) and 
(fw). For small excursions about a fixed point Eq. (3) 
may be linearized. The Laplace transform of the 
resulting relation is 


{45° + [B + CAyp + D+ E — F(fale) S + 
G (Afw/AM.) m-} AM, — (4) 


Rearranging this to obtain the desired transfer func- 
tion, 


MAyp, AS? + [B+ CAnp + D+ E — Fl faye) mcd S + 


This is a second-order, linear equation which can be 
written in terms of damping ratio and natural fre- 
quency as follows: 


AM, K (6) 


(S?/wo”) + + 1 
Where ¢ = damping ratio (dimensionless), and w) = 
natural frequency (rad./sec.). 

To account for the time lag between the by-pass 
doors and the control point, Eq. (6) should be multiplied 
by e~* where a is the lag time in seconds calculated on 
the basis of sonic propagation in steady flow. Eq. (6) 
then becomes 


AM, _ 
AAyp, (S?/wo”) + (26S/0) + 1 
Now damping ratio can be written as 
_ B+ CAy + D+ E = 
2V AG(Afo/AM.) ate 


(7) 


H 
(5) 
And, for the air where y = 1.4, damping ratio is 
VoA 
f= 0.343 4| 
V(Afy/AM,) 
+4 
AV, AV; 
1.4V,M, Ke(fu) | Vi 
V2 B.** 
(fx) me 


The determining factor in Eq. (8) is (Fac). 
(Fa)Me = Pu/Pu = (9) 


(Faye) = = 
(10) 


Referring to Fig. 2, which is a typical plot of inlet pres- 
sure recovery versus duct Mach number for steady con- 
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ditions, (a/c) 4, is the slope of the plotted curve and can 
be obtained by graphical measurement of the slopes. 
Thus, (Fz,,) can be obtained from the steady-state inlet 
performance curves. 

Inspection of Eq. (8) reveals that damping ratio 
will increase as (/4/,) increases in the negative direc- 
tion, and will decrease as (F4,,) increases in the positive 
direction. This phenomenon can be explained as 
follows: If (Fu,-) is negative, then: (a) w,/AW, is 
negative, and (b) AP,,/AW; is negative. Consider 
an inlet which is operating at a negative value of (Fy/,), 
for instance, point D on Fig. 2. Let there be a step 
increase in duct-exit flow area (this is done by opening 
the plug on a model or increasing r.p.m. on an engine. 
p, and P,, drop immediately, making AP,, positive 
and w, negative, both of which tend to increase W,,[{W, 
= W; — (—w,)|. As the low-pressure wave propagates 
forward to the inlet, the normal shock will move aft 
in the duct and P,;,; decreases, making AP,, negative. 
If the inlet is supercritical (choked) W; will remain 
constant. Otherwise, W; increases but lags W, by the 
amount w,. As AP,, becomes negative, the accelerat- 
ing force (AP;,) becomes smaller and the rate of in- 
crease in W, and W; becomes smaller and finally goes 
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to zero in an overdamped response. 

On the other hand, if (Fz,,) is positive: (a) w,/ATV, is 
positive, and (b) AP,,/AW, is positive. 

Considering now an inlet which is operating at a 
positive value of (Fu,,), for example point C on Fig. 2, 
let there be an increase in exit area to a value which 
will give a steady-state operation at point B. Immedi- 
ately p, and P,, drop, making AP,, positive and w, 
negative. The low-pressure wave travels forward to 
the inlet and W; starts to increase. As W; starts to 
increase, w, and AP, ; become positive because (Fu/,) is 
positive. The result is that W; leads W,, (W; = 
W, + w,) and AP,, increases so that the process ac- 
celerates. This will continue until W,; overshoots 
point B and travels over to the negative slope of the 
curve where it comes to a stop and reverses itself. 
Again, it will overshoot point B as it returns. 

Thus, it can be seen that the determining factor on 
damping ratio is the slope of the inlet recovery airflow 
curve. This characteristic is similar to that which 
causes compressor surge on the left side of the pressure- 
ratio peak on a compressor map. 

Referring to Eq. (A-6) of Appendix (A), natural fre- 
quency can be written: 


wy = VG(Afw/AM,)/A ~ rad./sec. (11) 
For air where y = 1.4 this becomes | : 


fo = = X 
| cps (12) 


— 


Examination of Eq. (11) reveals that the major influ- 
ence on natural frequency is the duct and engine 
geometry, and that duct Mach number has a minor 
effect. The effect of Mach number is to reduce natural 
frequency as Mach number increases. For the case 


where \/, is zero, Eq. (11) reduces to: 
fo = (V ygRT/28)V A/ (13) 


Eq. (13) is very similar to the equation for the natural 
frequency of a Helmholtz resonator [fp = (c/2r) X 
WV The Mach number term 

V (Afo/AM,) 
in Eq. (12) can be replaced by (ég/tm, the ratio of 
propagation time in quiescent air to propagation time 
in a moving stream. (g/t, is plotted against Mach 
number in Fig. 10. If tg/tm is to be used in deter- 
mining natural frequency, the average duct Mach 
number must be used rather than control-point Mach 
number. 

The work of reference 3 tends to discredit the 
Helmholtz theory of explaining duct resonance or in- 
stability. However, there is a basic difference in the 
model of this report which is believed to make the 
Helmholtz theory valid. That is, this model has by- 
pass ducts which act effectively as sudden expansions 
from the section ahead of the by-pass ducts to the 
section aft of the by-pass ducts so that the forward 
section acts as the neck or resonant column while the 
aft section acts as the damping mass. 

With the relations derived for damping ratio and nat- 
ural frequency, the response characteristics of an inlet 
may be obtained from its steady-state performance 
curves. Fig. 2 is the steady-state operating line, at 
the design Mach number, of a fixed-geometry sugar- 
scoop inlet model having by-pass ports for engine-inlets 
matching. The geometry of the inlet is as follows: 
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Fic. 3. Calculated dynamic response characteristics of a 0.2 
scale inlet model at design Mach number. 


A, = 0.25 sq.ft. A = 0.31 sq.ft. 
= Ve = Zleuk. Ve = 2.7 cuit. 


Using the relations just derived, the damping ratio and 
natural frequency are obtained and plotted in Fig. 3. 
The damping ratio varies from zero to 1.2 while natural 
frequency is fairly constant at approximately 18 cycles 
per sec. During the steady-state tests from which 
the data of Fig. 2 were obtained, the inlet was found 
to be unstable over a small range of supercritical mass- 
flow ratios corresponding to the calculated range of 
zero damping as noted on Fig. 3. The inlet normal 
shock was observed moving back and forth in the inlet 
throat during this instability. The inlet was, how- 
ever, self-starting. That is, the normal shock could 
be pulled back into the duct and stable operation ob- 
tained by moving the exit plug out so that the back 
pressure was reduced. The frequency of ‘“‘flutter”’ 
oscillations as measured by continuous recording of 
engine face static pressure was approximately 18 cycles 
per sec. with the doors open and 12 cycles per sec. with the 
doors closed. At much lower mass-flows (not shown 
on Figs. 2 or 3) the inlet “buzzed” at a measured 
frequency of 28 cycles per sec. with the doors closed. 
Evidently, “flutter,” which is an instability caused by 
internal flow phenomena, occurs at the duct funda- 
mental frequency while ‘‘buzz,’’ caused by external 
flow phenomena, occurs on a harmonic of the funda- 
mental frequency. 


Determination of Duct Dynamics 
From Transient Testing 


After the steady-state tests previously referred to 
were complete, one by-pass port on the model was 
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equipped with a solenoid-operated butterfly valve so 
that step changes in by-pass area could be made. 
Statham differential pressure transducers were installed 
in the model so that static pressures at the control 
station could be recorded continuously. Total pres- 
sure was obtained from an integrator rake, while static 
pressures were obtained from flush ports on the duct 
wall. All pressures plus butterfly position were re- 
corded simultaneously on a Minneapolis-Honeywell 
“Visicorder’’ with the paper speed set at 25 in. per sec. 

Transient tests were performed at the design Mach 
number and at several Mach numbers above and below 
design by putting step changes into A,, and recording the 
response of static and total pressures. At each Mach 
number, transient data were obtained for approximately 
fifteen different mass flows near the peak of the re- 
covery line. Testing was concentrated around the re- 
covery peak as the optimum operating point is near, 
if not on, the peak. Several representative traces are 
shown in Figs. 4 through 7. Point A is a trace obtained 
while the inlet was operating in the “‘flutter’’ region 
(Point A on Fig. 2). Both the butterfly valve and 
the duct mass flow plug were held constant during this 
trace. The amplitude of the oscillation is 4 in. Hg 
(+4 per cent of P;,) and the frequency is 17 cycles 
per sec. Points B and C are traces taken at subcritical 
points. The inlet is stable but underdamped at these 
points. Point D is a trace taken at a supercritical 
point. The inlet is overdamped here. Natural fre- 
quencies were measured directly from the traces, but 
could only be obtained where the inlet was under- 
damped or unstable. Damping ratios were obtained 
as follows. 

Static and total pressure at the control point were 
obtained from steady-state instrumentation prior to 
each transient test. Variations in static and total 
pressure during a transient were measured from the 
traces each 1/240 sec. Then P,./P,, (or Mach number) 
was plotted against time. From these plots, damping 
ratio was obtained. For the underdamped cases, the 
following relation was used: 


__ In (amplitude one/amplitude two) 


Where the outcome of this relation was doubtful, 
usually near critical damping, the amplitude of the 
first overshoot was compared to the known overshoot 
of a second-order linear response. For overdamped 
cases the damping ratio was estimated by measuring 
nondimensional time constant t/wo, assuming the nat- 
ural frequency to be the same as that measured from 
the underdamped cases. Further estimates of the 
overdamped damping ratios were made by visual com- 
parison of the plots with plots of second-order linear 
response to step inputs. 

The measured damping ratios and natural frequen- 
cies are plotted on Figs. 8 and 9 for the design Mach 
number and one Mach number below the design. 
Also plotted are the calculated response characteristics. 
The calculated damping ratios follow the same trends 
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Fic. 4. Pressure trace of unstable duct at point A on Fig. 2. 
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Fic. 5. Transient response trace at point B on Fig. 2. 
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Fic. 6. Transient response trace at point C on Fig. 2. 
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Fic. 7. Transient response trace at point D on Fig. 2. 
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Fic. 8. Comparison of measured and calculated dynamic re- 
sponse at design Mach number. 
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Fic. 9. Comparison of measured and calculated dynamic re- 
sponse at design Mach number minus 0.2. 
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Fic. 10. Ratio of wave propagation time in a moving stream to 
wave propagation time in quiescent air. 


as the measured values, but scatter in the measured 
damping ratios prevents determining the absolute 
values. The scatter is due mostly to aerodynamic 
‘noise’ in the traces. However, the calculated values 
and measured values do appear to agree to within 20 
per cent. Possibly the increasing discrepancy between 
measured and predicted damping as pressure ratio de- 
creases (supercritical duct) is the result of ignoring the 
effect of normal-shock position on the air-mass volume 
in the duct. An analysis of this effect is presented in 
reference 4. Note that the discrepancy is not so pro- 
nounced at the Mach below design (Fig. 9). At this 
condition, the inlet is unstarted and the normal shock 
never enters the duct. Too few measurements of 
natural frequency were obtained to determine whether 
the predicted variations with Mach number are cor- 
rect, but the general level of natural frequencies are 
comparable. 

During flight tests of the full-scale airplane, several 
records of duct “‘flutter’’ were obtained which showed a 
frequency of 4.5 to 5.5 cycles per sec. The calculated 
values were 3.5 and 5.0 cycles per sec., based on choke 
points at the engine exhaust nozzle and turbine nozzle, 
respectively. No measure of damping ratio was ob- 
tained. The performance of the complete control 
loop was satisfactory. The system was stable at all 
flight conditions and system response rates were near 
the predicted values, all of which implies that duct 
dynamic characteristics were approximately as pre- 
dicted. Unfortunately, the flight-test program was 
terminated before by-pass control tests were complete. 


Conclusions 


(1) Duct response characteristics calculated by the 
described methods agree to within approximately 20 
per cent with measured response characteristics of a 
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fixed-geometry supersonic inlet duct as obtained from 
transient testing in a supersonic wind tunnel. 

(2) Damping ratio appears to be practically inde- 
pendent of duct geometry. The determining factor 
on damping ratio is the slope of the inlet steady-state 
recovery-airflow line. For the negative slopes, the 
inlet is overdamped and for positive slopes the inlet is 
underdamped. On the peak of the recovery curve 
where its slope is zero, the damping ratio is approxi- 
mately 0.6 for the inlet tested. 

(3) Natural frequency is determined primarily by 
the duct geometry. Duct Mach number has a second- 
ary effect which reduces frequency as Mach number 
increases. Geometry effects are such that the fre- 
quency will decrease as duct dimensions increase. 

(4) The ‘‘sugar-scoop” type of inlet has, at or near 
its design Mach number, a small region of instability at 
slightly subcritical mass flows. The phenomenon 
occurs as the duct internal normal shock moves out of 
the inlet throat and appears to be a result of duct in- 
ternal-flow phenomena. This instability, or ‘‘flutter”’ 
is generally milder than the more commonly known 
“buzz” phenomena. However, at Mach numbers 
higher than design, the flutter amplitude increases 
rapidly. 

(5) “Flutter”? occurs at the duct fundamental fre- 
quency which can be predicted. “‘Buzz’’ appears 
to occur at a harmonic of the fundamental frequency. 
On the model tested, buzz frequency was approximately 
twice the fundamental frequency. 


Appendix (A)—Derivation of the Transfer 
Function for a Supersonic Inlet Duct 


It is desired to write the transfer function, A\/,/AA,,, 
for a supersonic inlet duct which uses bypasses to match 
the inlet and engine airflows. The control parameter 
is duct Mach number ahead of the bypass. 

The following assumptions will be used: 


(1) By-pass exits are choked. 

(2) Prop = 0.9P;¢ 

(3) W, is constant. 

(4) Flow is one dimensional and adiabatic. 

(5) Incremental flow variations occur isentropically. 
(6) Free-stream conditions are held constant. 


A one-dimensional ‘“‘lumped-constant” analysis will 
be used, wherein the air mass ahead of the by-pass will 
be considered as one homogeneous slug having average 
values at station c, and the air mass behind the by-pass 
will be another homogeneous slug having average 
values at station e. 

The effect of normal shock movement on the air- 
mass volume ahead of the by-pass is ignored for several 
reasons. The primary interest is in operation around 
the optimum point which is usually critical or slightly 
subcritical. For these conditions, the normal shock 
is in or forward of the throat. The sugar-scoop inlet 
retards shock movement forward. The duct contains 
a boundary-layer removal slot just inside the throat 
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which tends to retard shock movement aft during 
supercritical operations. Thus, for the conditions of 
primary interest, the shock movement is small and 
volume changes are correspondingly small. This 
simplification results in some error at supercritical 
points but is felt to be justified. 
The basic relations for transient flow are: 

Continuity 


W, = We + + (A-1) 
Pressure Drop From Inlet to Exit 


Fe == == AP», + AP,, 


AP» = steady-state resistance or ‘‘pipe”’ 
loss (A-2) 
AP, = differential required to accelerate ‘to 
the air mass in the control 
volume 


The terms in Eq. (A-1) may be rewritten in terms of 
mach number and total pressure as follows: 
W, = (WV = 
(fw) ~ (fw) KP, (A-3) 
W. = (fo) = 
K.P (M, assumed constant) (A-4) 
Wop = (fw) op(ArpP 11) = (A-5) 
W, = Ws, + Ws, (A-6) 
w, = (d/dt)(pV/RT) = | 
(V/RT)p — (A-7) 
Using the isentropic relations between p and 7, and 
between p and P;, w, may be written as a function of 
Mach number and P,. 
Ws = (V/yRT)) (1.411 P,/B*) (A-8) 
w,., = —  (A-9) 
w,, = K,,( P,,/B.2*) (M, is constant and M, is zero) 
(A-10) 


Now Eg. (A-1) may be put in this form 


— (4M — 
= 0 (A-11) 


A relation is needed now between P,,, P:;, P:,, and Px. 


Pi = Py (A-12) 


P,, = Pi, — APtr;, — APy,, [from Eq. (A-2)] 
(A-13) 


Pr = Px, AP tri-a BP (A-14) 
= = (A-15) 


Ke is the steady-state resistance loss factor (AP;/q) 
and can be estimated from published empirical data.’ 

AP,, can be determined as follows, ignoring duct area 
changes which are small: 


Force = Mass X Acceleration 


F = (w,/z)a ~ 
Ap; = (w./gA)a 
a=t=cM, w, = pv/RT 
Ap; = (pVc/gRTA)M = (pV/A)V y/gRTM 
Ap) = = 
But what we need is AP,;;; 
p = 

Apia = (Pu /Bi**) — (Pu/B,**) 
so, /Bi**) — = 
For subsonic Mach number and small perturbation, 


35 8,35 ~ 
B; Be B 


so AP, = 
AP1,,, = (A-16) 
AP ty. = Kiq oP (A-17) 


Combining Eqs. (A-12) through (A-17) yields 


Py = Pull — — Be? Md 
(A-18) . 


Pu P,,{l = (Ke,M2, (A-19) 


Because Eq. (A-i1) contains P,, and P;, terms, a rela- 
tion is needed for P,,. For constant Mo, the following 
relations can be obtained from inlet steady-state per- 
formance characteristics: 


P.,/P = 6/6) = f(M,) ~ f(M.) 
Pu (fo) tp 
Pi, = AP,,/At = (AP,,/4M,) (AM,/At) = 


(AP;,/AM) M, 
P,, = P.M. = (fa) (A-20) 
P,/P = (fa/fe)m-Me (faye) meMe (A-21) 


(fa)m, is the derivative with respect to M, of (f.)a, 
(P,,/P;,) and can be determined by graphic measure- 
ment of the slopes of the steady-state recovery curve 
at constant free-stream Mach number (Fig. 2). 
Combining Eqs. (18) and (20) yields 
(1 Kr, M2) /B.** (fae) (A-22) 
Pi, Pu = (fa)acPuMe = (faje)acP iM. (A-23) 
Now, by substituting the above relations for P,,, P:,, 
P,,. Pi. and P,, into Eq. (11), we get 
[K,,K;, Me + + 
pK Be? + + 
Mc { + 
fur) } (faye) me) Me + fu) we — 
KA fx) me = KopAop( fr) me (A-24) 


Inspection of the M, functions in Eq. (24) shows all ex- 
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cept (fa/-) and (fw) to be practically constant for small 
perturbations about a fixed point. 

Substituting constants A, B, C, etc., for the con- 
stant terms in Eq. (24), we can write 


F( faye) + G( fw) (A-25) 


Assuming small excursions about a fixed point, we write 
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AAM, + [B + CAypy + D + EJAM, — 
F( faye) — FM-A( faye) ae + 


Dropping the nonlinear terms gives this result: 


AAM, + [B + CA,, + D + EJAM, - 
F( faye) m-AM. + GA( fw) = (A-26) 


The Laplace transform of Eq. (26) is 
[B+ CA,, + D+ E- F( faye) mel S + 


Rearranging, G AM, = 

AA, AS? + (B+ CA, +D+ E — + G(Afw/AM.) 
3 This is the form of a second-order linear system where 

B+CA D+ E — F( fae 
Damping Ratio = ¢ = me 
2V AG(Afw/AM,) 
Natural Frequency = wy) = V (G/A)(Afw/ AM.) 
Static Gain = H/G(Afw/AM,) 
Now we can solve for and w. 
» Afw/AM) ] (8.25 
= (A cA 1/ V2V 142)0-926) (Afe/ AM.) 
= = 8.12 —_, ( ) ~ 
: The damping ratio is 
fu) me K,, x) me\ 
Be 5 
VA 2.13A.Vase 

J 
4 1.4ViM. . Kal 
where 


(fu) = (2M./B.**) — 


(faye) 
(fw) 
(Afw/AM-) 


(fx) = 1 — (KeM.?/6,**) 
(WV 


V [(1/B*) — | 


References 


1 Delio, Gene J., and Schwent, Glennon V., Transient Behavior 
of Lumped-Constant Systems for Sensing Gas Pressures, NACA 
TN 1988, December, 1949. 

? Henry, John R., Designs of Power Plant Installations Pressure 
Loss Characteristics of Duct Components, NACA Advanced Re- 
stricted Report L4F26, June, 1944. 


3 Trimpi, Robert L., An Analysis of Buzzing in Supersonic Ram 
Jets by a Modified One-Dimensional Nonstationary Wave Theory, 
NACA TN 3695, July, 1956. 

‘ Hurrell, Herbert G., Analysis of Shock Motion in Ducts 
During Disturbances in Downstream Pressure, NACA TN 4090, 
September, 1957. 


436 
| 
de 
th 
lo 
pr 
su 
pr 
sil 
sh 
for 
ap 
an 
res 
off 
thi 
up 
or 
lin 
tio 
dei 
th 
of 
sui 
an 
sel 
of 
det 
| 
na 
Di 
nat 
Ad 
nai 
ora 


6) 


Ds 


Theory of Supersonic-Propeller Aerodynamics’ 


DONALD EARL ORDWAY* AnD RICHARD W. HALE*™* 
Cornell University 


Summary 


A supersonic propeller with blades attached to an infinite cylin- 
der as a hub is studied. The forward speed may be subsonic, but 
the relative speed at each section is supersonic. The lightly 
loaded blades are represented by a surface distribution of ap- 
propriate ‘‘modified’’ sources in a fashion similar to ordinary 
supersonic thin-wing theory. These sources are found by ap- 
proximating the exact potential for a constant-strength compres- 
sible source traveling along a helical path. The usual relation- 
ship between the source strength and boundary condition is 
found; and subsequently the source distribution is given, to the 
appropriate order, in terms of the blade geometry. 

Tip effects are considered by extending the theory of Evvard 
and Krasilshchikova. The present investigation, however, is 
restricted to those planforms for which no vortex sheet appears 
off the tip. For points in the tip region, the potential is obtained 
through the appropriate distribution of ‘‘modified’’ sources in the 
upwash region off the tip. By transforming to a curvilinear, non- 
orthogonal coordinate system coincident with the modified Mach 
lines described by the infinities of the potential, an integral equa- 
tion for the required source distribution in the upwash region is 
derived. Without having to solve this equation, it is shown that 
the potential for a point in the tip region can be obtained in terms 
of an integration of known source distributions over the blade 
surface only. 

The case of a twisted flat plate of particular planform is treated, 
and a sample calculation is made of the pressure distribution at 
selected radial positions within the noncommunicating portion 
of the blade, as well as over the entire tip region. 

Though this analysis is carried out explicitly for the supersonic 
propeller, it could also be extended to calculate various rotary 
derivatives for highspeed flight vehicles. 


Symbols 
a = speed of sound 
Cre, Cre = coordinates of leading and trailing edge, re- 
spectively 
Cp = pressure coefficient 
D = distance from source to field point 
D = propeller diameter 
D, = hub diameter 
tf = source strength distribution functions 
K = M?/4Mrm7* 
K = the ratio of | cz | for numerical case considered 
k,q = source strengths 
1 = average chord to semidiameter ratio 
M = forward Mach number 
M = rotational Mach number 
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Mr = relative-free-stream Mach number—ie., 
(M? + M?)!2 

mr = (M,? 1)! 

Pn = a(t — 

Fs = radius of helix, or radial distance along blade 
from axis of rotation 

r.¢ = dimensionless, curvilinear, nonorthogonal co- 


ordinates in the xy plane of the tip-fixed 
x, ¥, g system, whose axes are given by the 
modified Mach lines through the origin of the 
local coordinates 


1D, Sb = field point in the upwash region, expressed in r, 
s coordinates 

rw, Sw = field point on the blade, expressed in r, s coor- 
dinates 

S = surface area 

¢ = time variables 

T(x) = projection of tip profile on the tip-fixed xy plane 

u,v, W = perturbation velocities in the x, y, 2 directions, 
respectively 

Uj, W; = induced inflow velocities in the x and z directions, 
respectively 

X, Y,Z = space-fixed Cartesian coordinate system 

X’, Y’, Z’ = translating and rotating propeller-fixed Cartesian 
coordinates 

X,¥Y,Z = dimensionless values of X’, Y’, Z’—ie., X = 

z, 7,0 = dimensionless cylindrical coordinates in the 
Z system 

= dimensionless, propeller-fixed, Cartesian co- 


ordinates, with the origin located at X’ = 0, 
Y’ = r,, Z’ = 0, and the x and y axes lying, 
respectively, along the tangent to the negative 
relative-free-stream direction at the origin and 
the Y’ axis 


Xp, Vp) Zp = source coordinates in the local coordinate system 

x’, y’,2' = local coordinates whose origin is located at the 
point x», ¥p, Zp in the x, y, z system 

a = angle of attack 

& = reference surface with respect to x, y, z 

ve = section mean camber line 

= coordinates of the high and low pressure blade 
surfaces, respectively 

tn, = position of the distributed source with respect 
to the local coordinates x, y, z 

20 = section thickness 

?, ¢ = perturbation velocity potential 

Q = angular frequency, assumed constant 

Introduction 


Ww THE ADVENT of modern high-speed aircraft, 
the propeller was generally considered to be out- 
moded as a means of propulsion. The problems under- 
lying this position were primarily twofold, one stemming 
from structural considerations and the other, from 
aerodynamic efficiency. Research over the past sev- 
eral years, however, has demonstrated that these prob- 
lems can be solved adequately. Combined with im- 
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Fic. 1. Source traveling along an arbitrary path in a compres- 
sible medium otherwise at rest. 


proved turbine engines, supersonic propellers offer 
definite promise for long range, high subsonic speed 
aircraft. 

This study develops a refined aerodynamic theory 
for a supersonic propeller with a finite number of thin 
blades. A propeller is called supersonic if either its 
forward speed is supersonic or its forward speed is 
subsonic, but its relative speed at each blade section is 
supersonic. The blades are assumed to be lightly 
loaded and attached to an infinite cylinder as a hub. 

There exists a considerable amount of research on the 
aerodynamic theory of supersonic propellers. As far 
as is known the initial work was by Frankl in 1942.': ? 
He described each blade by a straight vortex line with 
variable circulation along the radius and its associated 
spiral vortex sheet. The resulting velocity potential 
is expressed in the form of compressible doublets. 

Four years later, Hilton*® investigated a supersonic 
propeller with supersonic forward speed. He essen- 
tially used Ackeret airfoil theory for section lift and 
drag and maximized the section efficiency. The overall 
efficiencies in a few examples were about 70 per cent. 
Later, Weinig,’ Smith and Giarratana,° Borst,® Gilman 
and Crigler,’?’ Nilson and Smith,* and El Badrawy® 
carried out various studies utilizing the same two- 
dimensional section approach. 

In a paper published originally in 1950, Fabri and 
Siestrunck” performed a somewhat analogous calcula- 
tion. However, they included a modified lift-cancella- 
tion technique to provide hub and tip corrections and 
examined an estimation of the change in the zone of 
influence by using the local Mach angle at each radius. 
Considering only subsonic forward speeds, they also 
computed the velocities induced by the free vortices 
from a generalization of the conventional propeller 
theory calculation. Later, Borst!! postulated the effect 
of the modified Mach cone on these induced velocities. 

About this same time, three new analytical ap- 
proaches appeared. In the first, Burns! investigated 
a windmill operating at supersonic forward velocities, 
but at low rotational tip speeds, by expanding the po- 
tential in powers of the tip rotational Mach number. 
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In the second, Vogeley'*® outlined a compressible- 
actuator-disc theory for a subsonic propeller. After- 
wards in 1953 Delano and Crigler'! elaborated on this 
idea extensively, including the possibility of supersonic 
flow immediately behind the disc. In the third, 
Klawans and Vogeley” presented a cascade-general- 
momentum theory of a supersonic propeller annulus, 
treating the entire forward speed range from subsonic 
to supersonic. 

Several other aspects of supersonic propellers have 
been investigated; for example, general considera- 
tions,'*~'* sound and external pressure fields,?’~*4 
unique configurations,': structural problems,‘ 
and special applications.'*: ** In addition quite a num- 
ber of measurements have been made on overall char- 
acteristics, some of which are given in references 33-35. 

For this study, we are particularly concerned with a 
more accurate theoretical prediction of the pressure 
distribution on a supersonic propeller blade, taking 
into account the swirl of the incoming flow to a higher 
order in the linearized theory. In the first section, the 
exact potential for a helical acoustic source is developed 
and then appropriately approximated. These sources 
are distributed over a mean reference surface in the 
next section. Sections 3 and 4 relate the source dis- 
tributions to the blade geometry and develop solutions 
for the “noncommunicating’’ regions. In Section 5, 
the problem of the ‘“‘communicating’’ tip region is 
formulated, and, by a transformation to modified 
Mach coordinates, it is solved in the subsequent three 
sections. The final section presents the results of a 
typical example. To utilize our results for the calcu- 
lation of overall characteristics and to compare with 
the results using sectional data, the reader is referred to 
the recent comprehensive work by Smith.'* 


(1) Derivation of the Potential 


If the usual assumptions of small-perturbation theory 
are made, the potential ® for a source traveling along 
an arbitrary path in a compressible medium otherwise 
at rest is given by 
Y,Z, = 

1 q(Tn*) (1) 

| a Or,* H 
T is the time at which the potential is evaluated at the 
X, Y, Z field point in a space fixed Cartesian co- 
ordinate system; g(7) is the source strength as a func- 
tion of the dummy time variable r; and D(r) is the 
distance from the source to the field point (see Fig. 1). 
The 7,*’s, then, correspond to those times at which 
signals are emitted, such that propagating at the 
undisturbed speed of sound a, they will just arrive at 
the field point at time 7—1.e., they are the solutions of 
the equation 


D(r*) = a(T — (2) 


an 
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The term 4 
|1 + a 0D(r,*)/d7,*| | 
/ \ 
is the absolute value of the well-known Doppler factor. 3] ry — a 
In general, D(r*) # 0 and the infinities or Mach sur- | / be 
faces of the source, if they exist, are given by the zeros > E 
of this factor. | 
The generalized source potential of Eq. (1) was | / 5 isa 
derived in reference 36 by extending Prandtl’s technique a 9 =- 45° x 
for the rectilinear path.* This result has also been ob- nT ( } / 
tained independently by others, such as Blokhintsev,* | 7 ae 
Lagerstrom,*® and Westley.“ Moreover, the analogous 
LZ 
problem for an electromagnetic field due to a moving °. | > 3 4 5 6 
point charge was treated much earlier by Lienard,*! p 
Wiechert,** and Schott.* Fic. 3. Signals received at points along a radial line located at 
xX = 0,0 = —45°, from a source at ¥ = 0, 7 = 1, 6 = O in the 


To represent a propeller by an appropriate distribu- 
tion of sources, we will now specialize Eq. (1) to the 
case of a constant strength source moving along a helical 
path, such as shown in Fig. 2. Let X,, Yu, Zy de- 
scribe this path. Then, if / is the forward Mach num- 
ber, 7, is the radius of the helix and © is the angular 
frequency, 


= Mar 
Y, = fr, cos Qr (3) 
Zo = 7, sin Qr 


assuming that we arbitrarily fix the helix such that it 
intersects the Y axis at r = (0. Therefore, in this case 
we have 
D(r) = [(X — Mar)? + 

(Y — r, cos Qr)? + (Z — r, sin Qr)?]'” (4) 


and the 7,*’s are found from Eq. (2). Eq. (1) gives 


Ag 
a(l — M?*)r,* + M( Y sin Qr,* — Z cosQr,*)}-!| (5) 


where k denotes the constant source strength and , 
the rotational Mach number of the source or 7,Q/a. 
If r, > 0 or 2 > 0, this degenerates to the usual recti- 
linear source—e.g., see Rott.*4 

Of course, a propeller-fixed coordinate system is most 
Therefore, let us rewrite Eq. (5) with 


s 


convenient. 


Z 


Fic. 2. Constant strength source in motion along a helical 


path of radius 7, with forward Mach number J/ and angular 
frequency 


source-fixed cylindrical coordinates. 


respect to a translating, rotating coordinate system 
X’, Y’, Z’ in which the source is fixed at the point 
X'=0, Y’ =7,,Z' = 0. From Fig. 2, with r = T = 
t, the time at which the potential is evaluated in the new 
coordinates, 


X = X’+ Mat 
Y = Y’ cos — Z’ sin 2 (6) 
Z = Y’sin + Z’ cos Mt 
Substitution of Eqs. (6) into Eq. (5) yields 
MX — sin Mp, + Zcos Mp,)}—| (7) 
where X = X’/r,, etc., and p, = a(t — 1,*)/rs, ie., 


the distance a signal that arrives at the field point 
travels compared to the radius of the helix. Using 
Eqs. (4) and (6) in Eq. (2), we find the p,’s are given 
as solutions of 


p= \(X + Mp)? + P+ 2+ 
1 — 2(X cos Mp — Zsin Mp)}' (8) 


Thus, in principle, p, has become just a parameter and 
the potential at a point fixed in the moving system be- 
comes independent of time. Or, in other words, Eq. 
(7) together with Eq. (8) determines the potential for a 
source at the point 0, 7,, 0 in the Cartesian coordinate 
system X’, Y’, Z’, with a steady incoming, rotating 
stream of forward Mach number ./ in the negative X’ 
direction and angular velocity 2 about the X’’ axis. 
Now, in general, Eq. (8) cannot be solved for p be- 
cause of its transcendental nature. Therefore, we wish 
to make an appropriate approximation so that we can 
determine the potential explicitly. In the rectilinear 
case of M = 0, Eq. (8) reduces, after squaring both 
sides, to a quadratic in pif M # 1. Solving the quad- 
ratic, we find for 4 > 1 that there are two solutions, 


(M? — 1)pi2 = 
—MX + {X? (wt 1)? + 


inside the rear Mach cone defined by 
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Fic. 4. Comparison of Mach lines formed by the intersection 
of the Mach cone with the z = O plane as given by the approxi- 
mate and exact solutions of the present theory, Fabri, and 
Siestrunck,'® and the usual rectilinear flow solution. 


Furthermore, both of these signals contribute equally 
to the potential as shown by substitution into Eq. (7). 
In the region outside the rear Mach cone known as the 
region of forbidden signals, the potential is zero. But 
in our case with (7? + M *) > 1, we will usually have 
more than two signals; for example, see Fig. 3 which is 
a plot of Eq. (8) with ? = 0.5, M = 4,0 = —45°, 
and Y, Z replaced by polar coordinates 7, 6 or Y = 
F cos 6, Z = sin 6. However, in the neighborhood of 
the source, we should expect a type of modified Mach 
cone behavior, as typically indicated by the closed 
“loop” of Fig. 3. Though the potential may not be 
zero outside this cone if the advance Mach number is 
subsonic, it will correspond to signals emitted far down- 
stream and, in general, may be neglected. The excep- 
tions, of course, are the regions very near the singular 
Mach surfaces of that portion of the Mach cone after it 
“opens-up” forward. These infinities, though, result 
from the acoustic approximation and should, in the 
real fluid, be fairly well damped out. An indication of 
this is given by examination of the photographs of 
Hilton.*° Hence, we will: (i) consider only the contri- 
butions to the potential from the two smallest signals, 
and (ii) restrict ourselves to a neighborhood of the 
source on the order of the propeller chord-to-semi- 
diameter ratio, say O(/). 

In order to carry this approximation out, it is easiest 
if we transform to a local Cartesian coordinate system 
x, y, 2 which has its origin at the source, its y axis 
coincident with the Y axis and its x axis in the negative 
relative free-stream direction at the source. If we let 
Mr, = [M? + M2]/ ? denote the relative-free-stream 
Mach number at the source, then 


X = x(M/Mr) — 
Y=i1+y (10) 
Z = x(M/Mr) + 2(M/Mz) 


Substitution of Eqs. (10) into Eq. (8) and expansion of 
the sine and cosine terms for M = 0(1), and p < 1 
Mp < gives 
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pb = + Mrp)? + 
yt + Meyp? + (11) 


By setting M = Oin Eq. (8), we see that this is of the 
same form as for the rectilinear source, except for the 
M2yp? term. Hence, we can determine the two lowest 
signals correct to O(/*). Therefore, if we rewrite Eq. 
(7) with respect to the local coordinates, expand the 
sine and cosine again and substitute the approximate 
signals, we find that each contributes equally to the 
potential and 


P(x, y, 3) ~ —(R/2ar.) X 


— + 2°) — M2y(x? + y? + 2°) + 
(12) 


inside the modified Mach cone defined by 
—x = {mz%(y? + 22) + My(x? + y? + (13) 


where m7”? = (M,*? — 1). Outside this modified cone, 
@ =. This result was also obtained independently by 
Westley.” 

Examination of the shape of the modified cone reveals 
a swirl inward toward the axis of rotation. Fig. 4 
showing the intersection of the cone with the z = 0) 
plane illustrates this. In addition the exact solution 
determined numerically is given. Because of the im- 
portance of the intersection or modified Mach lines in 
defining the area of integration for a distribution of 
sources, the agreement is very promising. On the 
other hand, the approximation of Fabri and Siestrunck, '” 
which assumes that the slope of the Mach line is every- 
where given by the local relative Mach number does 
not agree so well. The usual Mach lines x = +mrpy 
result if the higher order terms of Eq. (13) are neglected. 

If Eq. (13) is extended beyond its range of validity 
for x, y, 2 < 1, the modified cone eventually opens 
forward as anticipated for a subsonic advance Mach 
number. However, since this occurs whether 1/ < 1 
or M > 1, this tendency is only accidental. 

Further justification of the approximate potential of 
Eq. (12) is given in reference 36 by the comparison of 
exact and approximate velocities at typical field points. 
For the typical case of M? = 0.5, i = 4.0, the dis- 
crepancy is negligible. 

It should be carefully noted that the lowest order 
terms inside the bracket of Eq. (12) correspond to the 
usual rectilinear supersonic source. The correction 
term — M2y(x? + y*? + 2°) is of one order higher and 
will be the order carried throughout. In fact, the 
criterion of carrying one higher order will be a guide for 
evaluating a consistent order in the subsequent work. 


(2) Distribution of Sources 


In order to treat the propeller problem, it is necessary 
to distribute sources over some mean surface corre- 
sponding to the plane ¢ = O in ordinary linearized, 
supersonic wing theory. Therefore, we must be able 


‘a 440 
; to 1 
fro1 
OL 
cal 
nev 
Vy, 2 
the 
will 
que 
ls b 
shiq 
new 
(10 
wine 
we 
our 
x’: 
} 
On 
nev 
the 
and 
The 
(12 
. 
insi 
P(x 
| 
| 
| 
4 but 
bra 
sou 
fort 
(x 
of 
rev 
sine 
mu 
} 
giv 
4 for 
the 
sur 


12) 


SUPERSONIC-PROPELLER AERODYNAMICS 441 


to write down the potential for a source located away 
from the origin of our local coordinates x, y, 2. 

Suppose we have a source at the point xp, Yp, Zp = 
O(1) corresponding to a point X,’, rp, 6, in the cylindri- 
cal system. Now, if at the new source we construct a 
new set of local axes x’, y’, 2’ in the same fashion as x, 
y, 2 were constructed at X’ = 0, Y’ = r,, Z’ = O, then 
the potential (x, y, 2; Xp, Vp, Zp) for this new source 
will have the same functional form as the old. Conse- 
quently, we can replace x by x’, y by y’, z by 2’, and 
r, by rp in Eq. (12), and, hence, if we find the relation- 
ship of the primed coordinate system to the unprimed— 
i.e., x’ = x'(x, Xp, Vp, Zp, Tp, etc., we will have 
the desired potential. This relationship may be found 
by constructing a new cylindrical system such that the 
new source is at the point 0, r,, 0, and then using Eqs. 
(10) as well as the simple relations of the hub-fixed 
Cartesian coordinates to the cylindrical coordinates. If 
we neglect O(/*) terms as required for consistency with 
our approximate potential, we find 


x’ = (x — xp) + (M?/My?)xp (vy —y + 
(MM/Mz*) — + 
= (y — yp) + + M2y)/Mr*] X 
{| M(x — x») + M(z — + O(/) 
2! = (3 — — (MM/M72)y,(x — xp) + 
(M/Mz*)(yp — y)(Mxy + Mzp) + 
On the other hand, the rotational Mach number at the 
new source and r, is needed only up to O(/) terms. In 
the above calculations, we find r,/r; = 1 + yp + O(/*) 
and, hence, 


M’ = (r,/r.)M = + y,) + (15) 


Il 


| 
| 


Therefore, using the form of the potential given by Eq. 
(12) and the relations of Eqs. (14) and (15), we obtain 
inside our modified Mach cone: 


P(x, Y, 2; Xp, Vp, Zp) ~ X 
{(x — xp)? — my*[(y — yp)? + (2 — — 
M*(y + yp) + — — 
M*(y — yp)(x? — — 
2MM(x — xp)(yZp — yp2) + (16) 


Eq. (16) will not be examined in detail at present, 
but it should be noted that the O(/*) terms inside the 
bracket coincide with the familiar rectilinear supersonic 
source displaced from the origin. We also see that the 
form is not the same as Eq. (12) with x replaced by 
(x — xp), ete.; this is caused by the skewed orientation 
of the primed and unprimed local coordinate system 
revealed in the O(/?) terms of Eqs. (14). Furthermore 
since the forward and rearward modified Mach cones 
must coincide, the potential is symmetric in x, y, z and 
Xp, Vp; 2p- 

Now, our aim is to distribute sources of the type 
given by Eq. (16) to represent a propeller blade, either 
for thickness or camber. Since in general the edges of 
the blade are mostly “supersonic,” the top and bottom 
surfaces of the blade do not communicate within the 


limitation of this approximate source potential, except 
at the hub and the tip which we will discuss later. 
Also to represent a lifting surface by sources instead of 
doublets, the reader is referred to Evvard*® or Krasil- 
shchikova.* 

Appropriately, let us choose a mean-blade surface to 
be defined by the x axis at each 7,—i.e., at each r, the 
surface lies in the local free stream direction. We 
take x = £, y = n, 2 = $ to denote a point on the mean- 
blade surface with respect to the local coordinates x, 
y, z. To find a point on this surface, we note that if we 
place the origin of our primed local coordinates on the 
radial line r, to the origin of our unprimed system, then 
points along the x’ axis generate our mean surface. 
Therefore, setting x,, 2p, y’, 2? =0,x = = yp 
and z = ¢, in Eqs. (14), we have from the third equation 


b> = (MM/M7?)nt + (17) 


Hence, the surface relative to the local coordinates has 
to this approximation a constant twist per unit length 
of MM/M,*. 

Let dS be an element of area of the surface defined 
by Eq. (17). Then the total potential $(x, y, z) due 
to a source distribution per unit area f over the mean 
surface is obtained from Eq. (16) by setting k = fdS, 


= yp = nands, = = This gives 
2m 
m,?MM 
mz*[(y — 9)? + 2] + ant — 
M,’* 


My + [(y — 2)? + 2) — 2) X 


—1/2 

(x? — + 2MM(x — + out fdtdn (18) 
where, consistent with our approximation, @.S has been 
replaced by r,’dédn since for ¢&; = O(/*) as given by Eq. 
(17), it can be shown that dS = r,*didn{1 + O(/*)]. 
The area of integration S is determined by the projec- 
tion of the intersection of the forward Mach cone from 
the point x, y, with the distribution surface. Eventu- 
ally, the field point will be placed on this surface, and 
so z = O(/*). Therefore, the z terms can be dropped 
consistently from Eq. (18) later. Or more simply, we 
will have to consider only (x, y, 0). Furthermore, as 
can be seen from Eq. (16), the area S, and the area 
determined by the intersection of the same cone with 
the z = 0) plane will be identical to the order retained. 
Now if we were concerned with only the direct prob- 
lem—i.e., the source distribution / is specified, we could 
evaluate the potential, calculate perturbation velocities 
and blade shape. However, our main interest is the 
indirect, or inverse, problem: given the blade geome- 
try, find f and the perturbation velocities and, hence, 
the loading. Therefore let us consider the determina- 
tion of f in this case. The boundary condition for 
inviscid flow requires only that at the surface the 
normal velocity component be zero. For ordinary 
linearized supersonic wing theory, this means that the 
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Fic. 5. A typical blade configuration attached to an infinite 
cylinder as a hub. Flow is viewed from a blade-fixed reference 
such that the free stream approaches with a forward Mach num- 
ber M and an angular frequency ®. 


ratio of the normal perturbation velocity to the free 
stream velocity must equal the local slope in the stream 
direction for both the upper and lower wing surfaces. 
Or from Eq. (18), we have, neglecting the O(/*) terms 


lim li ff 2 
— = — — lim — — t)?— 


OZ 
mr*[(y — + n)dédn 


where the left hand side is known and / has to be found. 
Rather than carry out the operation of the right hand 
side, the inversion is usually given in rigorous form by 
resorting to Green’s theorem. But this requires 
knowledge of the differential equation which the ele- 
mentary source solution satisfies. Since this is no 
longer known exactly for the approximate potential of 
Eq. (16), the differentiation and limiting process must 
be carried out in a straightforward manner. This is a 
lengthy calculation and is given in detail in reference 
36, including the ordinary case as well. The result is 
lim (0/02) ®(x, y, 2) = 


+ (r,/2)f(x, y)[1 + OW?)] (19) 


Thus, we see it is consistent with our approximate po- 
tential to employ the usual slope-source-strength rela- 
tionship because there are no terms of O(fl). The 
O(fl?) terms arise from: (1) distribution of the sources 
not in a plane but on the twisted surface expressed by 
Eq. (17), and (2) presence of a vertical velocity com- 
ponent at zs = 0 even if the distribution were in the xy 
plane (differentiate the integrand of the RHS of Eq. 
(18) with respect to z and set z = 0). An ordinary 
supersonic source does not have this property. The 
reason that the modified source possesses it is seen by 
referring back to the nonsteady coordinates—.e., the 
two signals which contribute to the approximate poten- 
tial did not lie in the plane when they were emitted, and 
so the above property is the result of the previous his- 
tory of the path of the source. (3) Rotation of the 
modified forward Mach cone due to the change in the 
local free-stream direction along the x axis. (4) In- 
clination of the distribution surface to the computed 
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velocity component 0®/0z|,_9, except at the point 
0, 0, 0. 


(3) Relationship of Blade Geometry to Source 
Distribution 


We shall consider the blades of the propeller of di- 
ameter 2 to be attached to an infinite cylinder of 
diameter 2, as a hub. The forward speed may be 
either subsonic or supersonic, but if subsonic, the hub 
diameter is assumed to be large enough to render the 
relative flow over the entire blade supersonic. 

The blade geometry is specified in the following way: 
Let the blade be centered about a radial line from the 
axis of rotation and let the distance along this line from 
the axis be designated by 7, as before. At r, we define 
the blade cross section in the ordinary thin-wing-theory 
fashion, separating camber and thickness with 


€(2r./D, 
o(2r,/D, 


Il 


as the ¢ ordinate of the camber line and the thickness 
distribution of the airfoil, respectively. Then, the ¢ 
ordinates of the surfaces consistent with our approxima- 


tions are simply 


= ‘ 
in = (21) 


The I surface will correspond to what is generally called 
the high-pressure surface of the section and II, to the 
low-pressure surface, for positive thrust in the forward 
direction and positive rotation in the same sense 
as before. To complete the description, we take 
Cre(2r./D) and Crpg(2r,/D) as the coordinates of the 
leading and trailing edges. Thus, we have a typical 
configuration as shown in Fig. 5 in which the incoming 
flow is visualized in the propeller-fixed coordinates. 

As we have said, we wish to represent the blade by a 
distribution or distributions of our modified sources as 
given by Eq. (18) in which as usual in the small per- 
turbation technique, we have placed the sources on the 
€ axis at each r,. Thus, if we require (1, fi = O(/*), 
Eq. (18) with z = 0 is a consistent representation of the 
blade since we only need to know the perturbation 
velocities on the surface to determine the propeller 
loading—i.e., the basic potential of Eq. (16) is unal- 
tered by the choice of z or z, if they are both of O(/*) as 
they will generate only O(/‘) terms which must be 
neglected. 

Also, we would like to satisfy the boundary conditions 
consistent with our approximations by considering the 
velocity components at ¢ = 0 in the é plane at each r,. 
Of course, there may be a radial component parallel to 
the nf plane, but with o and ¢, so limited and no radial 
“free stream”’ of O(1), this is of a much higher order and 
can be neglected. Furthermore, we can neglect the 
changes in direction and magnitude of the local free 
stream as we move along the & axis since they are 
O(aM for Che > &, x > Crp, i.e., &£,x = O(l). For 
generality, the inflow velocities uv; and w; in the — and ¢ 
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directions, respectively, due to both vortex wakes and 
interference between blades, are included. In general, 
these velocities only exist when the advance speed is 
subsonic and thus signals are capable of being propa- 
gated upstream, though for \/ > 1 there still may be 
blade interference for certain blade geometry and rota- 
tional speeds. The assumption is made that u; and w, 
are much less than a.\/7 and considerably smaller than 
u and w, the perturbation velocities produced by the 
proposed distribution of sources. What this means is 
that the perturbation velocities on any blade resulting 
from satisfying the boundary conditions there for an 
undisturbed free stream are in the main much larger 
than those induced in the approaching free stream by 
other blades and their shed vortices, as well as the 
blade’s own vortex sheet. This assumption is reason- 
able for calculating overall thrust, etc., but for finding 
optimum blade geometry, the inflow effects must be 
taken into account as in usual lightly loaded propeller 
theory. 

Now, from the boundry condition requiring no flow 
normal to the blade surface, we can write for the I sur- 
face 


wi(é, 0, w(é, 0, 1) 


— (22) 
+ O(?)] — u(é, 0, — 0, 


where the O(/?) term represents the change in the local 
free stream mentioned previously as we move along 
the — axis. Since ¢; < O(/*) and the derivatives 
ow ,/Ofr, Ow/OF1, and Ou/Of; are O(aM zl) or 
higher, we can expand the above equation about ¢ = 0) 
and neglect O(/*) terms. Using Eq. (19) to a consistent 
order, we obtain 


w(é, 0, 0+) Al, 0)/2 
and, thus, 


w,(&, 0, ()) + [fi(&, 0), 2] 
— (0f1/0£) — 0, 0) — 0, + 
(23) 


Let w;, u; = Ofor the moment. We still could not de- 
termine f1 directly from Eq. (23) as u is given as an 
integral of fr. But we do see 


fi/aMy = O(0¢1/0é) = O[20/(Cre Cre) | 


or O(a) where a is the angle of attack of the section. 
Therefore, u must be of this order since it is of the same 
order as w, and so 


0)/2 = [1 + O(0S1/08) ] 


Hence, for 0f1/0é or a, say, of O(/*), it is consistent to 
take with w; + 0, 


0)/2 = —aMy7(0s1/0E) — 0, 0) (24) 


For / and a of the same order we should use the ordinary 
supersonic source—i.e., neglect the O(/*) terms in the 
singularity. But this is just what we would expect, 


since the angles of attack and thickness ratios we are 
concerned with are about 0.1, and chord-to-semi- 
diameter ratios of this magnitude would mean very 
narrow blades for which the Ackeret theory must be 
very nearly correct. In fact, in this régime a surface 
distribution is unnecessary. However, from structural 
considerations, chord-to-semidiameter ratios of 0.3 are 
more probable. 

Now to determine the potential at a particular (£, r,) 
and, hence, the perturbation velocities, we will have to 
integrate over the sources intercepted by the forward 
modified Mach cone—.e., at each r,, we need to know 
fi(é, n). To do this, we note that at a different radius 
say r,’, the form of Eq. (24) is the same with r, replaced 
by r,’ and by But from Eggs. (14), &’ = & + 
and r,’ = r,(1 + 9), so 


firlé, n)/2 —wilf, = 
+ + )/D, —] (25) 


Likewise, 


fulé, = n, 0) + 
+ n(M?/ (1+ 9)/D, (26) 


(4) Solution of Certain Problems 


It is now possible without further calculations to de- 
termine the potential, hence, loading, for certain 
regions of the following cases: 

(i) Consider a single blade with thickness only and 
negligible inflow velocities. Then, w;, ¢, = 0 and from 
Eqs. (21), {1 = o and {1 = —o. Hence, using Eqs. 
(25) and (26), we have fi = fi where f; say, is given 
by Eq. (25) with w, equal to zero and ¢; replaced by co. 
The appropriate potential form for determining the sur- 
face pressure is given by Eq. (18) with z = 0 as we have 
noted, or, 


y, 0) = —(r./2m) ff — — me — 


My — — & + y? — n)dtdn (27) 


where S has become the area bounded by the projection 
of the leading edge of the blade on the & plane and the 
forward branch of the curve obtained by setting 
{(x — &)? — ...}!”? = 0, i.e, the projection of the area 
of the blade intercepted by the modified forward 
Mach cone from a field point located on the distribution 
surface. The solution would be valid over the whole 
blade except for a region near the hub where the inner 
wave of the modified forward Mach cone from the trail- 
ing edge of a section would first intersect the hub rather 
than the leading edge. See Fig. 6. 

If the integral of Eq. (27) is evaluated in its present 
form, the use of elliptic integrals will be required for the 
n integration since the integrand contains a reciprocal 
square root of a cubic in 7. In subsequent sections, 
however, we will find a way to avoid this complication 
to a consistent order. 

(ii) Consider a single blade again with negligible 
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Fic. 6. A typical blade plan form including extent of hub inter- 


ference region where the present solution is not valid. 


Fic. 7. Diagram indicating extent of tip communicating region 
for a blade of arbitrary plan form. 


inflow velocities, but with camber only or w;, ¢ = 0. 
Then, since we have shown that in the neighborhood of 
our source there exists a modified Mach cone, it is con- 
sistent to extend the idea, familiar in ordinary super- 
sonic thin-wing theory, of representing a lifting surface 
by a source distribution in the regions where the upper 
and lower surfaces do not communicate.’ These 
regions are found simply by looking at the modified 
rear Mach cones from the leading edge and determining 
where their intersections with the blade planform are 
just tangent to the leading edge. In general, the 
leading edge will be all supersonic except near the tip 
so that there will be only one communicating region as 
shown in Fig. 7. This will be defined more precisely 
later. Also, of course, we must exclude the hub region 
as before. 

So in the remaining region, ff = —f1: from Eqs. (25) 
and (26) with ¢; = {1 = ¢,and w, = 0. But with no 
communication between surfaces we can consider, say, 
the I surface, forget momentarily the given II surface 
and create a fictitious II’ surface such that it and the I 
surface form a hypothetical section with thickness 
only. Hence, we can find ®(x, y, 0+). Similarly, we 
can find ®(x, y, 0-), but since fr = —/fn, this is just 
the negative of (x, y, 0+). 


(5) Generalized Tip Problem 


The preceding section concerned itself with the non- 
communicating portions of a single blade by extending 
the conventional, supersonic thin-wing theory tech- 


niques through the use of the appropriate potential 
given by Eq. (25). The fact that upper and lower blade 
surfaces communicate in the region of the tip presents 
us with an inherently three-dimensional problem which 
would appear to obviate the use of the previously de- 
rived source potential. However, the technique offered 
by Evvard® and Krasilshchikova*® for the analogous 
case of ordinary rectilinear supersonic flow provides an: 
approach which will permit us to accomplish this with- 
out recourse to a doublet potential and its ensuing 
complexities. 

The concept introduced by Evvard and Krasil- 
shchikova was to consider a diaphragm or pseudo-wing 
of zero thickness to exist in the upwash field off the tip. 
If it ‘fellows’ the stream contours, the diaphragm will 
not alter the flow pattern. Note that the existing 
wing need be “‘extended”’ in this fashion only out to 
the foremost Mach line describing the upstream 
extent of the perturbed flow. This device then returns 
the problem to one where no communication at all 
exists between upper and lower surfaces of the blade 
providing the shape of the pseudo-wing, or more con- 
veniently, the required source distribution to replace 
this fictitious wing, is known. 

In the ordinary rectilinear-supersonic-flow case, 
Evvard and Krasilshchikova demonstrate, by examin- 
ing the behavior of the potential across the diaphragm 
(or across the upwash region in the real problem), that 
this source distribution can be obtained analytically. 
If the tip is a leading edge, in contrast to a trailing edge 
where a vortex sheet would appear off the tip, there is 
not only continuity in pressure but also in potential 
across the diaphragm. Writing the potential for points 
on the bottom and top surfaces of the diaphragm and 
equating them yields an integral equation of Abellian 
form for the required source distribution. Without 
having to solve this integral equation, they are able to 
write the potential for a point in the tip region of the 
wing in terms of an integration over the wing surface 
only, using a modified source distribution over part of 
the region of integration. 

In the case where a vortex sheet exists off the tip, the 
problem becomes considerably more difficult because 


Fic. 8. Sketch of modified Mach lines from points on a blade of 
arbitrary plan form. 
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of the unknown jump in potential across the sheet. 
Consequently, in our analogous problem we shall re- 
strict ourselves to the former case. This will therefore 
require us to specify the family of planforms for which 
our analysis is appropriate. 

Let us define the tip region of a single blade as it ap- 
pears in the rotating problem and to the order we are 
retaining here. Suppose we have a blade of arbitrary 
planform as shown in Fig. 8. Now, we identify the 
intersections of the forward and rear modified mach 
cones and the mean blade surface as the forward and 
rear modified Mach lines, respectively. Then, by con- 
sidering rear modified Mach lines emanating from the 
forward edge of the blade, we see that at a certain 
point O the outer rear modified Mach line OB is first 
tangent to the boundary of the planform and defines, 
to our approximation, the initial upstream disturbance 
off the tip. Similarly, the point C is the position along 
the rear edge of the blade at which a forward modified 
Mach line is first tangent to the planform boundary. 
These two points of tangency are evidently the points 
that separate the supersonic and subsonic edges of the 
blade. We will call the edge OC the tip. Hence, by 
definition, it is a subsonic edge, being either ‘‘leading”’ 
or “‘trailing’’ depending on whether the flow at a par- 
ticular point is directed on or off the blade. When we 
refer to the leading or trailing edge of the blade, as op- 
posed to the tip, we shall mean the supersonic forward 
and rear edges, respectively. The region between OB 
and the tip OC is the upwash field. The tip region then 
is defined as the portion of the planform between the 
tip and the inner rear modified Mach line OA. There- 
fore, only for points in this region do the forward modi- 
fied Mach lines include areas of the upwash region, 
which represents communication between upper and 
lower surfaces. 

To simplify our analysis, we have said that no vortex 
sheet is to appear in the upwash region off the tip. 
From the above description this clearly means that the 
tip must be a leading edge, requiring thereby that no 
streamlines can leave the tip but must instead leave 
from the supersonic trailing edges of the blade. In the 
limiting case the tip can follow a streamline of the un- 
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Fic. 9. Coordinate system and notation used in tip analysis. 
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Fic. 10. Relation between Cartesian and curvilinear coordinate 
systems. 


disturbed flow, thus appearing as the only symmetric 
planform possible under our restriction. 

For the remaining analysis, we set the origin of the 
local coordinates x, y, z at the blade tip, that is at 7, = 
D/2. We can determine quantitatively the equation 
of the undisturbed streamline which passes through the 
origin of the coordinates by finding the intersection of 
the mean blade surface given by Eq. (17) with a cylin- 
drical surface of radius 2/2 whose axis is coincident 
with the axis of rotation. If we neglect O(/*) terms as 
consistent with the order we are neglecting in the poten- 
tial, we find that the equation of this streamline pro- 
jected in the xy plane is given by 


y= —(1/2)(M?/ M,*)x? + O(/*) (28) 


If y = 7 (x) represents the projection of the tip in the 
xy plane of the tip-fixed coordinates, then 7(0) = 0 by 
definition, and in order for the tip to be a leading edge, 
T(x) is required to be a monotonic function of x having 
the property that (see Fig. 9) 


T(x) < —(M?/2M7*)x?, 0< x < Cur 
T(x) > —(M?/2M7?)x?, Cre <x <0 


The tip-fixed coordinate system which has been in- 
troduced is found to be most suitable for this problem. 
By virtue of its position on the boundary between the 
propeller surface and the upwash region, it best allows 
a coverage of both areas through one set of coordinates 
without violating the criterion that the variables be of 
O()). The only apparent change required as a result 
of this choice is to replace r, by D/2 in the potential 
given by Eq. (27), but it should also be noted that the 
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(S) OF S = Sp (r) 
Fields of integration used to evaluate the source 
strength distribution in the upwash region. 


Fic. 11. 


variables are now nondimensionalized by 2/2 rather 
than 7,. 


(6) Transformation of the Potential to Modified 
Mach Coordinates 


Setting the bracketed term in Eq. (27) to zero yields 
the equations of the projections of the modified Mach 
lines on the xy plane, but in the following we shall refer 
to these projections simply as the modified Mach lines. 
The process of finding an analytic solution for the tip 
problem is vastly simplified if a set of equations can be 
found which transforms the potential into a coordinate 
system along these modified Mach lines. This is the 
technique employed by Evvard, and its analogy in the 
rotating problem provides just the mathematical simpli- 
fication required to render a solution possible, at least 
without undue difficulty. 

Although the choice of axes and a system of measure- 
ment are somewhat arbitrary, it becomes expedient to 
choose modified Mach lines through the origin of the 
tip-fixed xy system as axes of the new nonorthogonal, 
curvilinear system. We shall use much of Evvard’s 
terminology in the following work, starting here by 
designating these axes r and s as shown in Fig. 10. 

We determine the values of r and s for any field point 
x, y in the Cartesian system by measuring along the 
axes of the 7, s coordinates, arc lengths from the origin 
to the points at which the modified Mach lines emanat- 
ing from the field point intercept these axes. 

Solving for the zeros of the bracketed term in Eq. 
(27), we determine the location £, 7 of the sources hav- 


ing modified Mach lines passing through the field point 
x, y. Or, by the reciprocity discussed previously, this 
is equivalent to the modified Mach lines emanating 
from a source at the field point x, y. Carrying the solu- 
tion out, we find 


= Fi2(n; x, y) = + 
1+ My 1) 
1+ 
-@ 


where = F,(n; x, y) corresponds to the family of 
modified Mach lines inclined to the left facing up- 
stream, and = F.(n; x, y) corresponds to those in- 
clined to the right. Then, if rw, sw in the 7, s coordi- 
nates refers to the field point x, y in the Cartesian sys- 
tem, we write from our above description 


{1 + [(d/dn) Fo(n; 0, 0) 1/2dy 


3 
Il 


(30) 


0 
f (1 + [(d/dn) F((n; 0, 0) 
m 


where 7; and 7 are the 7 coordinates of the points where 
the modified Mach lines from x, y or rw, sw intersect 
the r and s axes, respectively. From Eq. (29) we see 
that these quantities are given by the relations 


0, 0) F, (2; x, (31) 
Fi(m; 0,0) = Fo(m; x, 


If we solve for rw and sw in Eqs. (30) using Eqs. (29) 
and (31) and maintain one higher order than in Evvard’s 
development, we obtain to the appropriate order the 
transformation equations 


Il 


Sw = (My7/2mr)(x + mry) + 7/8m;?) Xx 
{ (2 — (1/myz?) jx? + — (2/my)xy} + 
rw = (Mz/2mz)(x — — (M?M7/8m7") X 
{ [2 — (1/mz?) jx? + y? + (2/mr)xy} + 


(32) 


Note that the O(/) terms are exactly the transformation 
equations in the Evvard theory. 

The inverse relations are needed in order to write the 
potential in the modified Mach coordinates. This 
can be done most easily by adding and subtracting 
Eqs. (32) and solving the resulting equations simul- 
taneously for x and y. From this we get, to the proper 
order 


(mr/Mr)(sw + rw) t+ 
(M?/2M7*mr)(sw? — rw?) + O(/) 
y = (1/Mr)(sw — rw) — (33) 


(M?/2M7°)(sw + rw)? + 
(M?/M,7*m7*)rwsw + 


Now in the tip-fixed Cartesian coordinate system, 
the potential given by Eq. (27) assumes the form 
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in which z = Oisimplied. The transformation of this 
potential into the nonorthogonal, curvilinear, modified 
Mach coordinate system involves transforming both 
the hyperbolic radius and the elemental area. 

Using Eqs. (33) and the analogous equations relating 
§, » to r, s we find that the hyperbolic radius is equal to 


2m Mr) ((s — sw)(r — rw)|'? x 
(1 + — rw +s — + 


(rw, Sw) = —(D/4e! 
g(rw, Sw) (D/4a My) V(s— sw)(r 


This form of the potential reveals the important advan- 
tage of using the modified Mach coordinates: we note 
that the first term in the denominator of the integrand 
is the hyperbolic radius in Evvard’s transformed poten- 
tial, and is the only singular term as long as sy, ... 
= O(/). Therefore, we can now expand the second 
term to obtain the following simplification of the poten- 
tial essential for its application to the tip problem, 


g(rw. Sw) = —(D/4rM 7) 
V(s — sw)(r — rw) 4M 
(sw —-rwts drds_ (34) 


If a similar separation and expansion were attempted 
using the potential expressed in Cartesian coordinates, 
the expansion would not be justified in the vicinity of 
the modified Mach lines. 

Eq. (34) clearly reveals our modification of one higher 
order to the Evvard theory in two ways: one, in the 
alteration of the Mach lines, and the other, in the ef- 
fective source strength of the distribution f. 


(7) Potential for a Point in the Tip Region in 
Terms of Known Source Distributions 


In Section (3) we introduced the notation I and II 
for the high- and low-pressure surfaces of the blade. 
Here we wish to extend this terminology to the pseudo- 
wing or diaphragm which we have constructed off the 
tip, despite the fact that it supports no pressure differ- 
ential. 

We recall that the fictitious wing has zero thickness. 
In the second part of Section (4) we showed that in 
this case fr. = —fu. If we let f*(r, s) represent the 
source strength distribution for the I surface of the 
diaphragm, then —f*(r, s) represents the distribution 
for the II surface. To obtain the integral equation 
that relates f*(r, s) to the known source distribution on 
the blade, we write the potential at an arbitrary point 


— &)? — — 9)? — M*%(y — (x? — & + y? — 9?) ]'” 


and the elemental area by the Jacobian, 
didn = |drds 


This is just Evvard’s result as the O(/) terms within the 
bracket which we must carry identically cancel. 

These results enable us now to write the source po- 
tential in the rs system. If we let ®(x, y) = ¢(rw, sw), 
and f(£, ») =f(r, s), we have to the correct order, 


fr, s)drds 


— rw) V1 + — rw +5 — 7) 


rp, Sp in the upwash field for both diaphragm surfaces 
by using Eq. (34) and Fig. 11. With no vortex sheet 
in the upwash field as we have stipulated, there is not 
only continuity in pressure but also in potential. 
Hence, equating ¢; and gy gives us the desired in- 
tegral equation. 

As in the Evvard case, it is not necessary to solve the 
integral equation. By an extension of his technique to 
account for the modification of the effective source 
strength, it is possible to show that the inner or r 
integrals are identical. Consider then Fig. 12 which 
shows the areas of integration for the potential at an 


Fic. 12. Fields of integration used to evaluate the potential at 
an arbitrary point in the blade tip region. 
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arbitrary point rw, sw in the blade tip region. Note 
that the area sw on the blade enclosed by the forward 
modified Mach lines from rw, Sw is divided into two 
subregions, Sw; and Sw», by the s constant line, s = 
S(rw), emanating from the intersection of the r = rw 


’ = — 
gr, Sw) 4a 


where K = M?/4\ypm,*. We see that the second term 
is identical for both g; and ¢gy, but for an infinitely 
thin plate, where = the integration over .S 
drops out entirely—i.e., the contribution to the poten- 
tial of the upwash field cancels the direct contribution 
of that part of the ble*te which generated the upwash 
field. Therefore, alteration of the planform boundary 
in Sw: does not affect the potential at rw, sw in this 
case. On the other hand, if we have a blade with thick- 
ness only, then fr = fur =/, ¢1 = ¢u = ¢, and ¢ is 
given simply by the integration of f over Sy. 


(8) Calculation of the Pressure Coefficient 


The linearized pressure coefficient for steady, com- 
pressible isentropic flow is 


C, = —(2u/U) (36) 


where the Cartesian velocity components are (U + u), 
v, w, and u,v, w< U. For the propeller problem we 
have seen that the stream direction of the undisturbed 
flow, which is incorporated in our potential, does not 
in general coincide with the x axis of our local coordi- 
nates. Hence. if we are to use Eq. (36) with U = a7, 
it is necessary to take the u perturbation velocity to be 
in the undisturbed free-stream direction, rather than 
the negative x direction, in the xy plane of our coordi- 
nates. This implies that v is everywhere normal to the 
undisturbed streamlines in the xy plane, while w is 
al.vays in the z direction. Over a large portion of the 
blade away from the tip, the differences arising from 
calculating u and v in this manner, which is consistent 
with our modified potential, as compared to taking 
simply the perturbation components in the x and y 
directions, respectively, are of an order to be neglected. 
However, near the tip, where the perturbation velocity 
may assume very large values, our more precise defini- 
tion becomes increasingly important. 

For the zero thickness problem, the tip presents an 
additional complication since it is a subsonic edge. 
The velocity component parallel to the tip edge must 
be continuous across it,*® and remains finite, but the 
components normal to the edge become infinite. This 
leads to the familiar phenomenon of the leading-edge 
suction force, which can be calculated by the usual 
procedure,*® and a breakdown in the validity of our 
linearized theory. We must conclude then that a truly 
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line and the blade tip s = s:(r). From this and the 
equality of the 7 integrals, we can incorporate the in- 
tegration over f* in the upwash field into an equivalent 
integration over a known modified source distribution 
in the region Sw. The results, derived in detail in 
reference 47, are 


D ff fi u(r, — K(sw — rw +s 
dyds — 
— Sw)(r — rw) 


2V/(s Sw)(r rw) 


consistent determination of C, cannot be obtained with 
our theory over a part of the tip region very near the 
tip itself. Nevertheless, we will proceed as in the 
ordinary supersonic wing theory with the fact in mind 
that in this region the solution may not represent the 
actual physical behavior. 

If 6/60 represents the directional derivative along a 
streamline passing through the origin of the local co- 
ordinates, we find to the propex order that 


6/80 = (0/dx) — [(M2/M72)x](0/dy) (37) 


Then for negligib'e inflow ve'ocities, we can use Eqs. 
(35) and (37) to find the perturbation velocities, m;, uy; 
in the undisturbed stream direction from 


u1,1 = —(2/D) X 
1/Or w) w/60) + (O¢r,11/O5 w) (6sw/60)] (38) 


which together with Eq. (36) determines the pressure 
distribution. It should be noted that the possible ef- 
fect of inflow velocities are simply incorporated in the 
determination of f; and fi as given by Eqs. (25) and 
(26). They are assumed to be known approximately 
from a first or higher iteration. 


(9) Example 


In order to illustrate some quantitative results from 
our theory, let us consider a constant chord blade 
having zero thickness and twisted such that the angle 
of attack at each section is inversely proportional to 
its relative, free-stream Mach number, Vy. We will 
also assume that the inflow velocities are negligible. 
Then, 


= —act/Mr 

a(2r,/D, —) = 0 (39) 
Ui, Wi = 0 

Cre(2r,/D) — Crx(2r./D) 


where «x is a constant of O(/). In addition, we shall 
specify a symmetric tip, conforming to the undisturbed 
streamline passing through the axis of the blade at r, = 
9/2. This will be analogous to the rectangular wing 
problem in rectilinear, supersonic flow in that u van- 
ishes at the tip. From Eq. (28) we see that in the 
tip-fixed coordinates the equation for the tip, 7 = 7(&), 
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Fic. 13. Chordwise perturbation velocity over the high- 
pressure surface of a twisted blade of zero thickness as a function 
of chord position and span location. 


is given by 
T(t) = + 


Using Eqs. (21), 1 = {11 = ¢., and from Eq. (25) with 
w; = 0, we find f1(&, 7) = 2aapo, and similarly, from Eq. 
(26), fir(é, = —2aay. 

Over the main body of the blade where the upper and 
lower surfaces do not communicate, it is convenient to 
use Eq. (27) to obtain a relation for (x, 0, 0) at each 
r,. This expression can then be differentiated directly 
with respect to x by using the following technique. 
In the usual terminology of finite parts* 


Ox J x: Vx —y x1 Ox Vx —y 


This suggests the change of variables ¥ = (x — y), so 
that 


x L(y) _2 L(x 7) 
Ox J x Vx ox Jo V5 

oL(x — 5) dy 
0 


L(x) 
x1 ox 


As shown in reference 36, the result of this differentia- 
tion can be integrated analytically for the present case 
where f1(£, 7) = 2aao in order to obtain an expression 
for uzy/aa. The results for various values of 27,/D 
are shown in Fig. 13 for M = 0.707, M = 2.50, and 
x = 0.12. We see that, except for a very slight dip, 
uy/aao remains approximately at the Ackeret value over 
the front half of the section. Over the rear half it con- 
stantly rises, reaching about 20 per cent above the 
Ackeret value at the trailing edge. In principle, this is 
what we might anticipate. That is, near the leading 
edge the forward modified Mach cone and the potential 


dy = 


(40) 


of the contributing sources are essentially the ordinary 
supersonic ones. However, as the distance from the 
leading edge increases significantly, they begin to de- 
viate considerably, the forward modified Mach cone 
opening up toward the hub and increasing the number 
of contributing sources. Of course, the opening up of 
the cone is somewhat nullified by the rotation of the 
cone toward the tip for decreasing x. Note that the 
average “1/dao decreases as 27,/D increases due to the 
decrease in angle of attack, ao/ Mr. 

The above solution for uj becomes somewhat simpli- 
fied if we use the potential as expressed in the modified 
Mach coordinate system of Eq. (34), which we derived 
for the solution of the tip problem. This avoids using 
elliptic integrals, but only at the expense of complicat- 
ing the limits of integration. 

To consider the tip region of the blade, we return to 
Eq. (35) in which a single coordinate system fixed at 
r, = 9/2 is used for all field points in the region. 
Noting again that the second term in the potential 
vanishes for a blade of zero thickness, we can write the 
potential ¢g; as an integration over Sw, with flr, s) = 
Fi(—, n) = 2aao. To evaluate u;, we use Eq. (38) where 
6rw/60 and 6sw/60 are found directly from Eqs. (32) 
and Eq. (37). As shown in reference 47, the quantities 
O¢1/Orw and O¢g1/Osw are obtained by first carrying 
through the indicated integration on r, then differentiat- 
ing the result with respect to rw and sw by means of 
Eq. (40), and finally carrying out the resulting integra- 
tion. The limits appropriate to the region Sy, are ob- 
tained to the proper order by use of Eqs. (32). 

Fig. 14 shows the variation of u;/aao over the tip 
region of the symmetric blade for the same case as con- 
sidered in the noncommunicating problem. The two 
solutions are joined along the inner border of the tip 
region as shown. 


Conclusions 


With the theory developed, it is possible to calculate 
the pressure distribution over a blade with both thick- 


Fic. 14. Perturbation velocity in the local free-stream direc- 
tion for the high-pressure surface of a twisted blade of zero thick- 
ness, as a function of position within the blade tip region; M = 


0.707, = 2.50,« = 0.12 uz = w; = 0. 
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ness and camber, except for certain regions influenced 
by the hub. Solutions for the propeller tip region are 
restricted to configurations for which no vortex sheet 
appears off the tip. This restriction, however, could 
be removed by a similar extension of the analogous 
Evvard theory. 

It is found that for the case of high aspect-ratio 
blades, Ackeret theory appears adequate for determin- 
ing optimum configurations as well as computing overall 
thrust and torque. For blades of moderate aspect- 
ratio, Ackeret theory is no longer adequate for optimiz- 
ing blade geometry, and it is for this application that the 
present theory should prove particularly useful. It is 
anticipated that blade-interference velocities, as well as 
wake-induced velocities for subsonic advance speeds, 
can be included through the use of an iterative proce- 
dure in any optimization analysis. 
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Oscillatory Aerodynamic Coeflicients for a 
Unified Supersonic-Hypersonic Strip Theory 


WILLIAM. P. RODDEN* anp JAMES D. REVELL** 


North American Aviation, Inc. 


(1) Summary 


This investigation presents a derivation of the oscillatory 
aerodynamic coefficients for wings with supersonic leading edges 
from the second-order, nonlinear, unsteady, supersonic flow theory 
of Van Dyke. The theory is considered applicable throughout 
the supersonic-hypersonic regime at Mach numbers normal to 
the leading edge and reduced frequencies for which B2/k > 4. 
The coefficients are modified for sweep, and a finite-span correc- 
tion is suggested to increase the accuracy of strip-theory flutter 
analyses. The limiting values of the coefficients in steady flow 
are also discussed. 


(2) Symbols 


V = free-stream velocity 

p = air density 

x,y = Cartesian coordinates, normal and parallel 
to free-stream direction 

29 = Cartesian coordinates, normal and parallel 
to leading edge 

t,n = Cartesian coordinates at tip 

L = local lift, positive down 

M = free-stream Mach number; local pitching 
moment, positive leading edge up 

M = Mach number normal to leading edge, 

M = Mcosa 

w = circular frequency 

é = local chord 

c = local chord normal to leading edge 

Ct = tip chord 

b = local semichord 

b = location of pitching axis aft of leading edge 


Xo = location of pitching reference axis aft of 
leading edge 


Ay = strip width 

h = local deflection of reference axis 

a = local angle of attack 

& = angle of attack normal to leading edge 

k = reduced frequency, k = wh/V 

q = free-stream dynamic pressure, g = (1/2)pV? 

ce = lifting pressure coefficient, positive down- 
ward 

g = local semi-thickness of airfoil section 

T = airfoil maximum thickness ratio 

Tte = airfoil trailing-edge thickness ratio 

Sm = chordwise location of point of maximum 
thickness 

r = 

¥ = specific-heat ratio of air 

B = (M?— 1)1/2 
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B = (MP? 

N = [(y + 1)/2](M//B)? 

N = [(y + 

A = sweep angle of leading edge 

nN = sweep angle of trailing edge 

ai; = coefficients in series expansion of pressure 
coefficient 

Ly, Ly = oscillatory aerodynamic lift coefficients due 
to plunging and pitching 

M,, Ma = oscillatory aerodynamic moment coefficients 
due to plunging and pitching 

d;, do = distances from tip of inboard and outboard 
sides of strip 

e = distance from leading edge at strip center- 


line to tip leading edge measured in free- 
stream direction 

thickness integrals 

coefficients in series expansion for oscillatory 


In, Ins Kz, L, = 
An, B,, Cu, D,, = 


aerodynamic coefficients 


F(Bn/é) = pressure factor inside of tip Mach cone 

Sn = finite-span correction factors 

Om = finite-span correction integrals 

En = coefficients in series expansion of F(8n/z) 
Runa = integrals in series expansion of Q, 

Cy. = local lift-curve slope, positive up 

Cue = local moment-curve slope, positive leading 


edge up 


(3) Introduction 


KF THE FLUTTER ANALYSIS of surfaces with super- 
sonic leading edges, a strip theory is desirable that 
will avoid the lack of conservatism of the coefficients 
of Garrick and Rubinow' as derived from linearized 
theory and which will be applicable at lower Mach 
numbers than piston theory.” The second-order theory 
of Van Dyke? offers such a possibility, and, in fact, it 
will be seen that both the linearized theory and the 
piston theory are special cases of his solution. 


The analysis of Van Dyke is the result of both an 
iteration and a frequency expansion in the solution of 
the nonlinear, unsteady potential equation. The 
iteration has been carried through to obtain the 
second-order velocity potential (first-order in thick- 
ness) and the frequency expansion has been extended 
to the cubic term. The final result has been presented 
as the series form of the local pressure coefficient on the 
airfoil. The present development transforms the 
solution of Van Dyke to obtain the pressure coefficients 
for pitching and plunging of the infinite swept wing. 
The pressure coefficients are then integrated to obtain 
the strip oscillatory aerodynamic coefficients, and a 
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where 
M 
a = (8/B°) (2 ~ 
+x Gis = (16M?/B*)(N — 1) 
ai = (8/B?) 2) 


(b) 


Re 


Fic. 1. Geometry of infinite wing. (a) Wing plan form. (b) 
Elevation. 


finite span correction is suggested in the tip region 
based on steady linearized theory. 

The use of the term “‘hypersonic’’ must be qualified 
as “hypersonic second order.’ The reason is that 
Van Dyke’s solution provides the first-order thickness 
correction and reduces to the corresponding piston 
theory results when 1/ ~ © in sucha way that Mr << 
1. As long as Mr is small, the nonlinear (second-order) 
potential equation upon which Van Dyke’s solution is 
based is a valid model. A third-order solution would 
necessitate the consideration of entropy changes; and, 
therefore, the velocity could no longer be written as 
the gradient of a potential function. Furthermore, 
the ‘‘smoothing’”’ procedure wherein the leading-edge 
initial conditions are idealized might have to be modi- 
fied to account for the bow-shock curvature in order to 
obtain a consistent approximation. In view of these 
problems, a low Mach number counterpart of the 
third-order piston theory terms would be difficult to 
calculate. The reader must, therefore, bear these 
comments in mind when the “hypersonic’’ comparison 
is discussed in Section (7). 


(4) The Pressure Coefficients for Pitching and 
Plunging Motions 


By doubling the upper surface pressure coefficient 
of Van Dyke [reference 3, Eq. (67)] we obtain the 
lifting pressure coefficient (positive down) for sym- 
metrical airfoils. The coordinate systems are shown 
in Fig. 1. 

If we consider the ordinate of the symmetrical air- 
foil to be g rather than eZ, we place e = 1 in Eq. (67) 
of reference 3 and obtain the following expression for 
the lifting pressure coefficient: 


= (Gor + Gog’) + i(k/é) (Gut + + + 
+ + + + 


das f + + + + 
+ + + 

& z 


dis = (8/B')(2 — M*)(M?N 1) 

dx = (4/B°)(M? + 2) 

dx = 16/B° 

23 = (8M7?/B*)(3(3M? — 2)N — 2(5M? — 3)] 
a = (8M?/B*)(4N — 5) 

= — 7172)N + 4(2M? — 3)] 
= (16M?/B*)(N — 1) 


a: 
Gx = + 2)N — 4] 
Gx) = + 4) 

Gx. = 


G33 = — 10472 — 4)N — 

— 2)(4M? — 1)] 

M? — 2)N] 

-1)- 
(12M* — 347° — 4)N] 

= + 1) + 8)N] 

= ((83M? + 2) — (3M? + 4)N] 

dss = (8112/5) + 1) — 

G39 = + 2) — + 4)] 


= (16 M?/B*) _ 
= 


= — 
N = [(y + 1)/2] 


If we note from Fig. 1(b) that 
b = % — 


then the pressure coefficient may be separated into 
its components due to pitching and plunging. 


Cp = (Cog + Co &o + Cp,h 


| 


where 
Coa = (Gor + dng’) + + Gag + 
Gyst + (R/E)?(Gnt + Gos gdt + 
Z + + 1(R/E)* (Gut + 


and 


Cos, = —Co, = i(k/€) (aye Guz’) + 
+ Gof + + 


The two-dimensional pressure coefficient due to pitch- 


ing and plunging motions has been found. In Section 
(5) the effect of sweep on these coefficients is considered. 


(5) The Swept-Wing Pressure Coefficients 


We consider now the transformations necessary to 
obtain the pressure coefficients for the infinite swept 
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wing from those obtained in Section (4) The geom- 
etry of the swept wing is shown in Fig. 2. 
From geometrical considerations the following trans- 


formations are obtained: 


= =xcosA 
Xo = Xo cos A 
€ = 2bcosA 
M = Mcosa 
& = a/cosA 
g’ = g’/cosA 
gd =qcos’A 
and since CA = 
then C, = C, cos? A 


We note that h, k, and g are invariant. After making 
the above transformations in our previous results, we 
obtain the pressure coefficients for the swept wing. 


Co = (Cra + + 
where 


Cra = (do + dog’) + 1(k/2b) X 
(Qux + dig + aisxg’) + + 


23 f gdé + + dux*g’) + 
1(k/2b)*(asix* + as Egdé + 
35X gdé + + 


and 


Cox, = —Con = 1(k/2b)(Qi2 + aug’) + 


x 


+ gdé + azexg + 


Qu = An cos A 
ay = A 
Qi = G@ecos A 
Qo, = Ay cos A 


a2, = A» cos A 
= cos A 
= COS A 


and all remaining aj; = diy. 

We have obtained the pressure coefficients for the 
infinite swept wing due to pitching and plunging mo- 
tions, and are now in a position to derive the oscillatory 
aerodynamic coefficients. 


(6) The Oscillatory Aerodynamic Coefficients 


The strip oscillatory aerodynamic coefficients are 
defined by the following equations: 


+ 
4 pw%b*(M,h/b + Maa) 


dL/dy 
dM/dy 


The lift and moment per unit width are found from the 
pressure coefficients 


Fic. 2. Swept-wing geometry. 


2b 
dL/dy = af 
0 


2b 
af + + dx 
0 


26 
dM/dy = q f (x — 
0 


2b 
qd f (x Xo) [(Cre + xoCp,,) + Cynh| dx 
Comparing the above two sets of relationships we find 


Ly, 


2b 
(1/8k?) 
0 


R 
Il 


2b 
(1/4k2) (1/26) f (Cra + xoCr;,)dx 
0 
2b 
M, = (1/4k?)(1/26) f (x — xo) Condx 
0 


M. 


2b 
(1/2k?) (1/26)? (x — X0)(Coa + XoCpx,)dx 


These expressions can be rewritten in the following 
more convenient form: 


Le = Lay — (x0o/b)Ln, 
M, = My, — (x0/b)Ln, 
Ma = May — (Xo/b)(My, + Lax) + (x0/b)?LZy, 


where the leading-edge values are 


2b 
Condx 
0 


2b 


(1/4?) (1/26) 
0 


8 
II 


2b 
(1/4k?) (1/26) f xCpdx 
0 


II 


2b 
= (1/2k?)(1/2b)? f x 
0 


Next we define a set of dimensionless thickness integrals 
as follows: 
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~ 
ll 


2b 
0 


2b 
Jn = dx,n = 1(1)4 
0 


2b x 

= f gdé,n = 1(1)3 
0 0 
2b 

0 0 


These integrals are calculated for a typical airfoil in 
Appendix (A). 

By carrying out the indicated integrations we ob- 
tain the leading-edge values of the oscillatory aero- 
dynamic coefficients. 


Lng Az + 1(A3/k + 
Ly = (Bi/k? + Bs) + i(B3/k + Bik) 
= Gs 1(C3/k + Cyk) 
Ma, = (Di/k? + Dz) + i(Ds/k + Dik) 


where 
Ay = —(1/8)(dx2/2 + + 
A; = —(1/8)(ai2 + Sidi) 
A; — (1/8) (32/3 + + 1303s) 


By = (1/4) (ao + 

= (1/4)(au/3 + + + 

Bs = (1/4)(a@u/2 + Siais + 

By = (1/4) + Lids; + + + 


= — (1/4) (d2x2/3 + + 

C3; = — (1/4) (a2/2 + 

CG = — (1/4) (a32/4 + Keds, + + 1403s) 

D, = (1/2) (am/2 + L292) 

= (1/2) + + + [4Q27) 

D; = (1/2) (an/3 + Jodi3 + 

= (1/2) + Leas; + Kya35 + + 


In the usual manner, the aerodynamic coefficients are 
referred to the quarter-chord location by taking x) = 
b/2. Therefore, the quarter-chord oscillatory aero- 
dynamic coefficients become 


Ly = Az + 1(A3/k + A,k) 
Le = (Bi/k? + By — A2/2) + 
i[((B; — A3/2)/k + (Bs — 
M, = (C2 A:/2) + A;/2)/k + 
(Cy — A4/2)R] 


M, = |(Di — B,/2)/k? + De 
(Bz + C2)/2 + A2/4] + 
i{[D; — (Bs + C;)/2 + As/4]/k + 
[Ds — (By + Cy)/2 + A4/4]R} 


We have derived general expressions for the two- 
dimensional oscillatory aerodynamic coefficients for 
the case of the supersonic leading edge modified for the 
effects of sweep and thickness. In Section (7) it will 
be demonstrated that the present results agree with 
all of the known limiting cases. * 


* Van Dyke has already demonstrated compatibility of his 
fundamental pressure equation with linearized theory and with 
Busemann’s steady state second order theory. 


(7) Limiting Values of the Aerodynamic 
Coefficients 


As a check on the previous development we will 
investigate the limiting values of the aerodynamic 
coefficients. We will consider the unsteady cases, 
first, of zero sweep and thickness, second, of high 
Mach number, and, finally, the steady case with sweep 
but zero thickness. 

In the case of zero sweep and thickness all of the 
thickness integrals vanish and the quarter-chord 
coefficients become 


— (1/8%) + a[(1/Bk) — 


—(1/Bk?) + (5M? + 1)/68° — 
i{ (M2? — 3)/26k + [M2(4M? + 1)/687|R} 


M, = —(5/68) — i[(1/28k) — (M2/B°)R| 


—(1/2Bk?) + (3M? + 1)/46® + 
i{ (17 — 7M*)/126% — [M?(41.M? + 14)/606"] R} 


Ln 


M, 


The above expressions are based on a solution ob- 
tained by a series expansion in terms of reduced fre- 
quency; therefore, the coefficients of the present theory 
are not directly comparable to those of Garrick and 
Rubinow which are the results of a closed form solution, 
and we must resort to numerical comparison. We 
make the comparison with the quarter-chord coeffi- 
cients tabulated by Luke.‘ The comparison is pre- 
sented in Table 1. It should be noted that the values 
of Luke are normalized to the factor 7 rather than 4 
as is used here. 

An inspection of Table | shows that our series solu- 
tion agrees quite well with the linearized theory at all 
reduced velocities (1/k) for the higher Mach numbers 
and at high reduced velocities for the lower Mach 
numbers. An estimate of this correlation is given by 
B?/k > 4 and, accordingly, this condition is considered 
to be the practical limit of the present formulation. 
It should be noted that this limit specifies the minimum 
Mach number normal to the leading edge rather than 
free-stream Mach number. 

We remark that, for very low aspect ratios, the root 
effect becomes predominant, and the use of the 
unswept strip theory is necessary to achieve agree- 
ment with three-dimensional linearized theory—i.e., 
Ci, = 4/8 for the wide delta wing rather than 4 
cos A/B. Because of the additional restriction on the 
reduced frequency range as 6 ~ 0, the streamwise theory 
permits the use of the theory for a wider range of k. 
These arguments are believed to be the stronger when 
considering low aspect ratio wings. It is believed 
that the cross-flow theory produces more conservative 
results, but no conclusive comparisons have been 
made, and the flutter engineer should consider both 
possibilities in any particular application. 

To obtain a correlation with piston theory we con- 
sider the limiting case of M*? > 1 and neglect terms of 
order 1/M? in the expressions for the aerodynamic 
coefficients. We observe the limiting values 
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Bo M 
+ 1)/2 
—4/M 
ado, —2(y + 1) 
ay, —8/M 
—> 8/M 
—ay 4(7 + 1) 


d2;, and a3; are all of the order 
The aerodynamic coefficients 


and that the values d;, 
1/M?, 1/M?, or 1/M‘4. 


become 


—i(1/Mk) (1 + (1/2)(y + 


These coefficients are found to be in complete agree- 
ment with an unpublished investigation of piston theory 
by the present authors. 

The final check is on the effect of sweep. 
sider the case of steady flow, zero thickness, and in- 
finite aspect ratio. The linearized steady-flow pres- 
sure coefficient inboard of the tip Mach cone is given 
by Donovan and Lawrence’ to be 


C, = —4a(M? — sec? A)~“/” 


We con- 


Our solution in this case is 


Ly, = Cp = dna = (dm cos A)a 
L, = —(1/Mk?)[{1 + (1/2)(y + — = —4a cos A(M? cos* A — 
i(1/2Mk) [1 + (1/2)(y + DM(4h — = 4a(M? — sec? A)~“/” 
M, = —i(1/2Mk)[1 + (1/2)(7 + DM(4he — and agreement is again found. 
It has been shown that the present solution agrees 
M,= —(1/2Mk*)[1 + + with all known special cases: the steady and oscilla- 
1(7/12Mk) [1 + (3/14) (y + 1) X tory linearized solutions and the nonlinear solution of 
— + piston theory. It therefore can be concluded that the 
Comparison of the Present Values of Oscillatory Coefficients for Zero Thickness to Those Obtained by Linearized Theory 
Present Ly Luke's Present La Luke's 
10/9 2.0 -8,80233 + 119.0345 -0.57607 - 1 1.60430 36.60203 - 1 82.19608 - 3.76512 + 1 0.07210 
+* -8,80233 + 1 3.3230 -2.22951 - 4 3.22) 11.82513 - 1 17.78817 = 15.37998 + i 5.64206 
-0 -8,80233 - 110.7269 -6.20787 = 111.91630 - 87.28250 + 1 37.72563 = 97.00688 + 4 41.58997 
20.0 -8,80233 = 138.978) -8.326,0 = 139.06280 -781.03580 + 1145.62548 -782.76708 + 1145.89556 
S/u 2.0) 2.37037 + 1 0.62551 -0.48397 - 1 1.34093 0.85596 = 1 3.66468 - 3.11567 - 1 0.11797 
-2.37037 1 1.62090 = 1 4.03373 - + 1 3.27877 = 16.590) + 1 3.87099 
=-2.37037 - i 9.84362 -2.15843 - 1 9.8909 - + 1 11.86160 79.54277 + 4 11.9198 
20.0 -2.37037 - 126. -2.33522 = 126.306 -527.1Ujh00 + 1 33.36685 -527.21)22 + 1 33.37157 
10/7 2.0 -0.94176 1 1.03710 -0.43769 - 1 1.36277 2.23116 1 = 2.91501 - 1 0,05655 
4.0 -0.94176 - 1 3.4591) -0.78469 - 3.50696 13.99351 + 1 1.12925 - 1).19913 + 4 1.18611 
-0.9)176 - 1 7.61075 -0.90022 1 7.61695 - 61.0290 + 1 3.27458 - 61.09612 + 3,28218 
20.0 -0.94176 = 119.51159 -0.93510 = 119.51165 -390.38880 + 1 8.89771 +-390.39786 + 1 8.89856 
5/3 2.0 -0.42188 - 1 1.17041 -0.264)0 - 1 1.24554 2.41111 - 1 0.28047 2.57391 - 0.2097 
ee -0.42188 - 1 2.83521 -0.37707 - 1 2.84550 = 11.41113 + 1 0.00038 - 11.45676 + 1 0.01005 
0.42188 - 1 5.91760 -0.41029 - 1 5.91892 - 7.41113 + 0.281h - 47.42313 + 1 0.2827) 
20.0 -0.42188 -0.42003 114.96710 -299.41113 + 1 0.9000 -299.40949 + 1 0.90321 
2 2.0 +£0.192h5 - 1 1.02640 -0.13845 - 1 1.04719 - 2.08488 - 1 0.31362 - 2.13110 - i 0.29712 
0.19245 = 1 2.24525 -0.1777h 1 2.24797 9.01308 - 1 9.02548 - 1 0.44336 
-0.19245 - 1 4.58673 -0.18873 - 1 4.58712 36.72589 - 1 0.80009 36.7291, = i 0.7995) 
20.0 0.19245 = 111.53418 -0.19179 - 111.53420 -230.71600 - 1 1.93662 -230.71856 - 1 1.93993 
M Present Mh Luke's Present Ma Luke's 
10/9 2.0 7.33527 + 145.76300 -0.01956 1 0.5225 39.99799 - 1 94, .08967 - 1.70169 - 4 1.0901 
-7.33527 + 119.78439 -0.61858 - 1 0.16611 27.60953 1 28,65233 2.67742 - 1 0.9691 
00 + 1 3.69796 -4.5310h 1 3.84350 21.94427 + 1 2.45886 = 33.09275 + 1 27.0962 
20.0 -7.33527 = 115.86465 -6.81333 - 116.42937 -368.82090 + 1111.98163 -370.8257, + 1112.30,06 
5/4 2.0 1.97531 + 1 3.0562 -0.05011 - 1 0.25502 2.27160 - 1 5.21962 - 1.04529 - 1 1.06021 
-1.97531 = 1 0.47189 -1.16553 1 1.4213 5.72840 + 1 0.98278 6.31931 + 1 1.68939 
00 1.97531 1 4.2359) -1.74351 - 1.4.56568 37.7280 + 1 7.67657 37.1341 + 1 7.77544 
20,0 1.97531 = 112.8937 -1.93679 = 113.00776 -261.725 + 1 23,18913 +-260.7443h + 1 23.19281 
10/7 2.0 0.78480 4 0.07519 -0,24826 - 1 0.43079 0.34925 - 4 1.0180, 1.11307 1 0,55198 
hed -0.78480 = 4 1,50789 -0.61340 - 1.55328 6.230h3 + 1 0.12862 - 6.46383 + 1 0.19792 
-0, Te - 1 3.69454 -0.73922 - 1 3.69726 - 29.7551, + 1 1.34311 29.51686 + 4 1.3526 
20.0 -0.78480 - 1 9.71146 -0.77676 1 9.70988 +-19).42810 + 4 4.11592 ~-194.43317 + 1 4.11549 
5/3 2.0 -0.35156 = 4 0.50281 -0.18158 - 1 0.50729 = 0.94629 - 1 0.5670h = 1.12854 - 1 0.45318 
-0.35156 1.37640 -0.30254 - 1 1.34719 - 5.4629 - 0.54133 - 5.49779 - 1 0.53014 
-0,35156 1 2.93820 -0,33882 - 1 2.91909 23.4629 - 1 0.78629 = 23.45906 - 1 0.7950 
20.0 -0.35156 - 1 7.47528 -0.34950 = 1 7.46757 -149.4h629 - 1 1.75827 -149.45342 - 1 1.75929 
2 2.0 0.16038 - 1 0.50328 -0.10171 - 1 0.47312 = 0.94621 - 4 0.47970 - 0.99785 = 1 0.602) 
-0.16038 1 1.11766 -0.14428 - 1 1.09375 4.41031 - 1 0.76909 - 15 - i 0.76655 
20 -0.16038 - 1 2.29088 -0.15629 - 1 2.27773 - 226672 = 1 1.44302 - 18,27072 - 1 1.44278 
20.0 -0.16038 - 1 5.76610 -0.159h) 1 5.76090 -115.26159 - 1 3.54094 -115.25718 - 1 3.54215 
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Fic. 3. Tip region geometry. 


present solution offers a unified two-dimensional super- 
sonic solution suitable for strip theory flutter analysis. 

The nature of the steady-state pressure distribution 
given by linearized theory inside of the tip Mach cone 
by Donovan and Lawrence suggests one final modi- 
fication to our oscillatory coefficients for the tip strip 
or strips within the Mach cone. Section (8) considers 
this correction. 


(8) A Finite-Span Correction 


The steady linearized supersonic solution for the 
swept wing inside of the tip Mach cone is given by 
Donovan and Lawrence® for constant spanwise angle 
of attack as 


Cp = 4aF(Bn/t)(M? — sec? A)~°/” 


where 


| 
= (2/m) arctan 


and the geometry in the tip region is defined in Fig. 3. 
The nature of the steady-state pressure coefficient in 
the tip region suggests the use of a correction factor to 
the oscillatory pressure coefficient such that the strip 
lift is found from 


L = dyS CyF(Bn/£)dx 
and the strip moment from 
M = qJS'dyS (x — x0)C)F(Bn/é)dx 


in which the factor F(8n/é) is taken as unity throughout 
the region inboard of the tip Mach cone. The integrations 
are performed over the entire strip. The use of this 
steady-state factor is justified by the fact that super- 
sonic flutter usually occurs at a low value of reduced 
frequency so that a quasi-steady analysis may be em- 
ployed. It should be remarked that the use of this 
correction on the nonlinear terms is a convenience and 
is justified only by the semiempirical nature of the 
correction. 

This use of the pressure factor results in correction 
factors as follows: 
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fi = 1 — (1/2bAy)Qi 
= 1 — (1/2b?Ay) (Qe + €Q:) 
fs = 1 — (38/8b*Ay)(Qs + 2eQ2 + 
ts =]— (1/4b4Ay) + 3eQ; + + 
fs = 1 — (5/32b'Ay) X 
(Qs + 4eQ; + 6e?Q; + 4e*Q2 + e4Q,) 

where 

—1 n tan A+c, 


and the distance e is measured from the centerline of 
the leading edge of each strip to the leading edge at 
the tip. These factors are used in the terms com- 
prising the aerodynamic coefficients such that modified 
values may be substituted for the two-dimensional 
values in the following manner: 


A; = SiAs C;’ SoCs 
BY’ =fiB, Dy’ = 
B;’ = Dy = 
Bi = = 


These factors can be considered to be the ratios of 
the three-dimensional to the two-dimensional volume 
and its first four moments of the strip pressure dis- 
tribution diagram about the strip leading edge in 
steady flow. 

To facilitate carrying out the integrations required 
for the tip correction factors, it is convenient to seek a 
polynomial representation for the function 1 — F(@n/é). 
A plot of the function shown in Fig. 4 for various 
values of (1 + tan A/8) shows a cubic representation 
to be satisfactory. By passing the cubic through the 
points Bn/é = 0, —0.2, —0.8, —1.0, we obtain the 
polynomial form 


1 — F(Bn/t) = Eo + Ei(Bn/é) + 
E(8n/é)? + E3(Bn/é)* 


where 
Eo = 1 
E, = (1/0.096)(0.8F, — 0.2/: + 0.096) 
E» = (1/0.096)(1.8F, — 1.2F, + 0.6) 
E; = (1/0.096)(F; — F. + 0.6) 
F, = (2/m) arctan [(1/2)(1 + tan A/)'””] 
= (2/m) arctan [2(1 + tan 


The cubic curve fits are also shown in Fig. 4. From the 
polynomial form we may write the required integrals 
in the form 


where 
—do n tan A+ ct 
Run = fan f 
—di —Bn 


The integrals R,,, are presented in Appendix (B). 
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We see from Fig. 3 that d; cannot exceed ¢,/(8 + 
tan A). Therefore in the evaluation of the factors for 
the innermost strip affected by the tip Mach cone the 
value of d; must be taken as c,/(8 + tan dX). 

We have seen that a low-frequency tip correction 
is readily derived for strip theory. It is anticipated 
that its use will increase the accuracy of the preceding 
strip theory in supersonic flutter analyses. 


(9) Concluding Remarks 


The present study serves to emphasize the importance 
of the nonlinear unsteady supersonic flow solution 
obtained by Van Dyke and demonstrates that it is 
suitable for strip theory flutter analysis by the addi- 
tion of the sweep and finite span corrections. The 
significance of this solution is seen in its general ap- 
plicability throughout the supersonic-hypersonic re- 
gimes, the minimum normal Mach number being de- 
termined by B?/k> 4. 

The present results are considered as complementary 
to piston theory which provides a simpler solution, 
which though restricted to higher Mach numbers, has 
been extended to account for control surface deflections 
(cf. Ashley and Zartarian*), and can account for third- 
order effects (second-order thickness correction). 
Because of its simplicity, it would appear to be the 
logical avenue of attack on the calculation of the 
effects of camber and more general airfoil thickness 
distributions. The analysis of Landahl® is designed 
to extend the lower Mach number range of piston 
theory and has the advantage of retaining most of the 
simplicity of the piston theory results, but is not ex- 
pected to cover the range of Van Dyke’s solution. 

Although the present analysis suggests a finite span 
correction at the tip, the effects of fuselage interference 
and the inboard Mach cone have not been considered. 
The inboard Mach cone provides another lower Mach 
number limit as the zone of influence extends over a 
greater portion of the lifting surface as the normal 
Mach number approaches unity. Neglecting this 
effect is conservative because the lifting pressure is 
lower inboard of the root Mach cone than it is in the 
two-dimensional flow region. It is suggested that neg- 
lecting these effects may not be overly conservative 
since aeroelastic applications are concerned primarily 
with the load distributions on the outboard portions 
of the surface. 

The present oscillatory coefficients can be utilized 
in flutter analysis by the collocation method if the 
matrix of aerodynamic influence coefficients is found. 
Reference 7 formulates the aerodynamic influence 
coefficients in terms of the quarter-chord strip coeffi- 
cients. A similar formulation can be made in terms 
of the leading-edge coefficients. The present coeffi- 
cients can also be used in static aeroelastic studies by 
noting the identities: 


Cia = —lim (4k7L,) = lim (2k?M,) 


k—0 ko 0 


in which the coefficients are referred to the local 
chord, Ci, is positive up and Cm, is positive leading 
edge up. 

It is recommended that further work be done in ex- 
tending the solution to arbitrary frequencies for various 
specific but representative airfoil shapes and for more 
general downwash distributions, including control 
surfaces and camber modes. Van Dyke has indicated 
the procedure to be followed and has illustrated the 
calculations for the case of a wedge for pitching and 
plunging motion where, for example, the trend of the 
thickness correction reverses at higher reduced fre- 
quency, for a wedge pivoted at its leading edge. Ex- 
perimental correlation is suggested to determine the 
range of practical usefulness of the method, as well as 
the validity of the suggested tip correction. A pre- 
liminary correlation of the finite span correction on 
some untapered swept plan forms has shown that the 
use of the tip correction reduces the calculated flutter 
speeds. This indicates that the reduction in tip load 
is more than offset by the forward movement of the 
center of pressure, yielding a more critical situation. 


Appendix (A)—The Thickness Integrals /, J,, 
K,, and L,, 


Expressions for the thickness integrals are derived 
for an idealized typical airfoil section. The idealized 
section will be composed of two parabolic arcs giving 
a sharp leading edge and a blunt trailing edge. The 
airfoil geometry is shown in Fig. 5. 
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Fic. 4. Lifting pressure correction function in tip Mach cone 
region. 
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Fic. 5. Idealized airfoil geometry. 
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The equation for the forward parabolic are through 
the points (0, 0), (x, (/2) and horizontal at the latter 
point is 

gi(x)/c (r/2)(x/Xm) (2 X/¥m) 
where 7 is the airfoil thickness ratio, (/c. The equa- 


tion for the aft parabolic are through the points 
(Xm, t/2), (€, tte/2) and horizontal at the former point is 


g(x)/e = (7/2)[1 — — — — 
Tre/T 
The airfoil slope is given by 
g1'(x)/c = — x/¥m) 
and 
= —7(1 — r)(x — Xm)/(€ — Xm)? 


The required integrals can now be calculated from the 
expressions for the local thickness and slope. These 
are tabulated in Table 2; the values for the biconvex 
airfoil (r = 0, x»/c = 0.5) are given also. 


Appendix (B)—The Integral R,,,,, 


No general expression for the integral R,,» is avail- 
able because of the numerous combinations of m and 
n. Therefore, a tabulation is given in Table 3. Spe- 


cial consideration is necessary in the case of small 
trailing-edge sweep for Ru, Rw, Rie, Re, Re, and R3;, 
since the basic formulas lose significant figures. Ac- 
cordingly, series expressions for these integrals based 
on the series expansion of the logarithmic terms are 
given at the bottom of Table 3. Ten terms of the 
series have been found to be adequate for |d; tan d/c;! 
< 0.25. 
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TABLE 2 
The Thickness Integrals J,,, Jn, Ky, and Ly, 


Integral Functional Value _[ Note: =(l-r)/(1-y)? J 


Biconvex Value 


= (1/e) ("grax = tr/2 0 

Ip = (1/c?) = (t/6)[ y-2(2-3yty3)] 
13 = (1/03) = 
Ty, = (1/elt) ("x3g'ax = (T/20)[ y3-2(h-Sy+y5)] -37/10 
Ig = (1/25) ("alg tax = /15 
= (2/62) gax = (t/6)[2+r(1-y)] t/3 
= (1/63) ("xgax = (t/2h)[6-y2-2( t/6 
J3= (1/e4) (Px2gax = T/10 
Jy = (1/65) ("x3gax = (t/120) [15-y4-z (10-2 y+15y2-y6)] t/15 
= (1/e3) ("ax = (1-y)2] t/6 
Ko = (1/elt) | ("eas = (-15y+20y2-10y3+y5)] 77/60 
K3 = (1/05) ("x2ax("gae, = 10-36y+hSy2-20y 3+y%)] 
= (1/04) (Pax = (T/120)[10-Sy2+2y3=z ( 3-10y+10y2-syt+2y5)] 


Lo = (1/05) ("xax ("¢ eat = (T/2h0) [15-Sy2+yt-z t/20 
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The Integral Rn», 


(1/2){ 204 (p+tand) } 


+(c¢/tand) 


p°{(1/2) (1/6 +1/tand) (442-42) (dy tanN/(e¢-d, tand))} 


(p3/2) {( 1/2) (4,2-4,2) (dy-d,)=( In[(c¢-d, tan) tand) 
-1/(cy-d, tan) ] } 


= (1/6){3e42 (di } 


= {(1/3) (dy 74,3) (Betand)-(c4/2) } 


{- (1/4) (ay } 


~(e¢/2tand) (4q-4,)} 


(1/12){ yo (dy dy 24,2) dy (p>+tan JA) } 


(B/2)§ (1/4) (ag (p?-tan®d) +(2/3) (443-453 eg 


+(04/3tand) (dy 3-453) (4, (dy } 


(1/20){ )-100 3tani dy 24,7) +100 4, (4, 
(B/3){ (1/5) (4, (B>+tanrr) =( 3/4.) 3/2) 
(B°/2)§-(1/5) (dy = (1/2) (ay 4 egtand+(e 42/3) (dy } 
B?{(1/5) (44 } 
(1/30){ 664( 44-4, )-150 +200, 3tan2A( dy 3-4, 
+60 dy )-( dy (p>+tan®d) } 
-(c,4/2) 
(B2/3){-(1/6) (p>+tan3d) +(3/5) (44 tan®-( 3/4) (4, Je 43/3) (43-4, 3)} 
(2/2) { (1/6) +(2/5) (847-895 

Rm for small 


(8/2) tan +( tamd¥ 3043) (dy 3-4,3) 
+(tan®@ J} 


= 4 (1/2) (442-492) (1/3043) (dy +( tan (dy J} 


+( (ay 5-d95)+( (dy o-a,6)+... J} 


(p°/3 ) Ady] -d,31n tand) +e43[( 1/3e¢>) (dy 324,3)+(tanVke ( ) 
+( tan@d/Se,°) (4,5-a,5)+. ee } 
B?{-(1/38) (443-403) 4043 (dy (a4 J} 


}} 
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Further Comments on ‘‘Limiting Circulatory 
Lift of a Wing of Finite Aspect Ratio’’ 


H. S. Ribner 

Professor, Institute of Aerophysics, University of Toronto, 
Toronto, Canada 

January 19, 1960 


Ts REFERENCE 1, McCormick obtains the limiting circulatory 
lift of a wing with elliptic loading as 


(Cr) maz = 1.21A, A = aspect ratio (1) 


In reference 2, Hancock takes exception to some of the arguments 
and obtains 


(Ci) maz = 0.855A (2) 


The present note confirms Eq. (2) for the case of a planar ulti- 
mate wake, and it points out an alternative higher result for the 
case where the wake rolls up into a vortex pair. The basis is 
reference 3, which apparently was not known to either author. 

From an exact momentum balance, the lift coefficient of a 
wing with elliptic loading is obtained as [see Eq. (21) of reference 
3] 

= (1/2)7A sin — (1/4)7A sin? (3) 

where ¢ is the inclination of the wake far back, assumed plane. 
The maximum value, which occurs for sin e = 12/3, is 


(Ci)maz = V20A/3V3 
= 0.855A 
Thus, Hancock’s value is confirmed for the case of a planar ulti- 
mate wake. 
For the case where the trailing vortex sheet is fully rolled up 
into a vortex pair, the assumptions of reference 3 yield a different 
resuit, namely, 


Cr = (1/2)xA sin € — (1/4)7A sin? € sin e’ (4) 
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where e’ is the modified inclination angie of the rolled-up pair, 


. 


given by Eq 


sin ¢’/sin = 2/7? (5) 
The maximum, which now occurs for ¢ = 7/2, is 
(Ci) maz = (9 — ')A/2 (6) 
= 1.412A 


Eqs. (1) and (2), as pointed out in reference 2, are not appli- 
cable to a jet-flapped wing; likewise, neither is Eq. (6). How- 
ever, Hancock? has made the necessary extension to Eq. (2) and 
found the added terms to be small. The experimental data‘ 
cited in references 1 and 2 lie above Eq. (2) for the planar wake 
(and above its extended form), but below Eq. (6) for the rolled-up 
(The full predicted maximum of Eq. (6) is probably 
it is predicated on a wake angle of 90°). 
not 


vortex pair. 
unattainable since 
Thus, perhaps the rolling up of the trailing vortex sheet, 
previously considered in this problem, plays an important role. 
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The Problem of Strain Accumulation Under 
Thermal Cycling 


B. E. Gatewood 
Research Coordinator, Air Force Institute of Technology, 
Wright-Patterson Air Force Base, Ohio 


January 27, 1960 


Poazzes' and Sprague and Huang? have shown that it is possible 

for strain growth to occur in a beam structure under tempera- 
ture-load cycling. The various aspects of this problem as to 
criteria for convergence and divergence of the strain accumula- 
tion can be simply demonstrated by thermal cycling one element 
of a two-element structure. 

Consider two bars to be fastened together and assume no 
bowing or buckling. From page 165 of reference 3, the strains 
in the two elements are 
+ er + Cap 
= —(aT), + er + Cay 

E\A\(aT), + E,A+(aT 
— 


+ 


= —(aT); 


(1) 


where 7 is temperature, @ is thermal-expansion coefficient, E is 
secant modulus, A is area, éa, is applied strain, and subscripts 1 
and 2 refer to elements 1 and 2. With eq, as an elastic strain, the 
average applied stress on the cross section is 


Fa, = + E2A2)/(Ai + A2)] 


where FE, and £, are the elastic moduli. 

For a given value of the average applied stress, the strain can 
be calculated from Eq. (1), whence the stresses are obtained from 
the stress-strain curves for the material at the given tempera- 
ture. Figs. 1 and 2 show the results for two different average 
stresses on 2024-T86 aluminum alloy with A; = A2; AT = 300°F. 
temperature change on one element, and the same stress-strain 
curve is used for both elements. No creep is considered. The 
cycle used is: apply temperature (+7)), apply load (+ P,), re- 


(2) 
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move temperature (—7;), apply temperature (+72), (—7>), 


(+T7;), (—T3), ete. 

Fig. 3 shows the average stress, F,, against strain for various 
numbers of cycles, including elastic shakedown at which point 
the strain growth .or accumulation ceases. It should be noted 
that these curves in Fig. 3 may be used for design just as ordinary 
stress-strain curves. They give the allowable stress on the total 
cross section for any desired permanent set with due allowance 
for temperature effects on material properties, thermal stresses, 
and inelastic effects. This idea of determining an average allow- 
able stress for design when nonuniform stresses exist on the 
cross section is similar to that commonly used for allowable 
crippling stresses at room temperature. This average stress is 
also the easiest value that can be obtained from test data. 
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Torsion of an Elliptical Shaft With 
Diametrically Opposite Flat Sides 
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Engineering Sciences, Purdue University, Lafayette, ind. 


December 28, 1959 


SUMMARY 


The solutions for the torsional stress and torsional stiffness of elliptical 
shafts with flat sides which are presented in this paper have been ob- 
tained by means of the Galerkin technique. The solutions for both shear 
stress and torsional stiffness are given as functions of the ratio of semi- 
axes and of the ratio of flat-side distance with respect to major semiaxis. 
A comparison of the results obtained is made with previously given liter- 
ature. Close agreement in the solution is obtained. Numerical results of 
torsional stiffness and of maximal shear stresses for different shape ratios 
are tabulated. 


FORMULATION OF THE PROBLEM AND ITS SOLUTION 


TT PROBLEM Of the elastic torsional behavior of prismatical 

bars has been studied extensively by many investigators. 
Torsion of a circular shaft with diametrically flat sides has been 
treated by Carter and Oliphint, using relaxation methods.! 
Recently, Stanisi¢é and Shirely presented a solution of torsion of 
a cylindrical bar with a trapezoidal cross section using the 
Galerkin method.? The torsion of an elliptical shaft with dia- 
metrically flat sides is an important problem in modern aero- 
nautical structural analysis, since structural weight plays an 
enormous role. In this paper, the solution of this problem will 
be approached through the Prandtl stress function using the 


Galerkin method. The procedure of this method for determi- 

nation of stress and displacement fields has been outlined previ- 

ously by Stanisi¢ and Shirely.2 Therefore, only the results of 

mathematical analysis of this problem will be presented here. 
Assuming stress function in the Prandtl sense, 


n 
W(x, y) = Annn(x, y), (m = 1,2,3,...,”) (1) 
m =1 


the Galerkin form of the stated problem becomes 


feo yln(x, ydS = 0, (yr = 1,2,3,...,m”) (2) 


where 
y)] = + 2 


The 9m(x, y) are coordinate functions which satisfy the boundary 
conditions, ¥(T) = 0, on the boundary I; A,,’s are arbitrary 
constants which are to be determined from Eqs. (2); and S is the 
domain of integration. Hence, Eqs. (1) and (2) lead to 


DY Annm(x, »| +2 a(x, = 0 (3) 
Ss m=1 


Eq. (1) can now be written as 
W(x, = y)Pm(x, (4) 
where 
m(T) = 9 (5) 
and 
Prix, = Ay + Aox? + Asy? + +... + A,x?y" (6) 


Note that P,,(x, y) need contain only even terms since the sec- 
tion of Fig. 1 is symmetrical with respect to both coordinate 
axes. For P(x, y) = Ai, the Galerkin method yields 


A, = ff no(X, vas/ ff [V 2n0(x, y)] no(x, (7) 
J's 


However, by virtue of Eq. (5), it follows: 
nox, y) = (x? — h?)(b*x? + a®y? — a*b?) (8) 
Then, 


-2ff y)dS 
a 


— 16(h/a)? — 3] — (h/ay? + 
3[1 — 6(h/a)*] sin~! (h/a)} (9) 


and 
ff LV 2no(x, y)dS = (a%3/360) X 
S 


{(h/a)[45 + 165(b/a)? — 210(h/a)? — 1154(b/a)%h/a)? — 
440(k/a)* + 648(b/a)*%h/a)* + 80(h/a)® — 
304(b/a)(h/a)}] — (h/a)? + 15[—3 —11(b/a)? + 
16(h/a)? + 16(b/a)%(h/a)? — 48(h/a)! — 

sin~' (h/a)} (10) 


Egs. (9) and (10) yield A; and, hence, \,,,). 


The torsional 
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TABLE 1 


The Values of Shear Stresses and Torsional Stiffness 
S 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Tnaz/Gak 1.26 1.25 1.25 1.24 1.24 1.21 1.16 1.10 
D/Gh4 2.28 2.13 2.06 1.07 1.85 1.70 1.8 


rigidity is given by 


D= 26 ff W(x, y)dS (11) 


The stress tensor is given by 


0 0 GaWy(x, y) 
[ry] = 0 0 —GaW,(x, y) (12) 
GaWy(x, vy) —GaW¥,(x, y) 0 
Hence, 
tzr = 2GaA,a*y(x? — h?) \ (13) 


tTzy = + y? — (b/a)*(a? + f 
Thus, the displacement field is characterized? by 


y) = (2/3)Ai[(a? — 65?)x8y — 
a*xy® + 3(a%* + bh? — ath®)xy] + B (14) 


where 
—ays 
[ui] = (15) 
a®(x, y) 


and B is an arbitrary constant that corresponds to the rigid body 
displacement in the z direction. The displacement field, Eq. 
(15), is determined by Eq. (14) to within rigid body motion as 
stated by the first boundary value problem in static elasticity. 
The results for the stress and displacement field can be easily 
evaluated according to Eqs. (13) and (14) for any given geometry 
under consideration. 

For h = a, an elliptical cross section of semiaxes a and 3 is 
obtained. Then, 


—2 fom y)dS = (16) 
Ss 


f [V2no(x, qo(x, y)dS = + 4362) (17) 
Ss 


and 


Thus, 
D = xGa*bh*/(1.05a? + 1.296?) (18) 


Eq. (18) compared with exact solution® gives a close approxi- 
mation. Furthermore, if a = 6 = h, we obtain a circular section 
with torsional rigidity, 


D = nGa‘/2.34 (19) 


which also is in close agreement with exact solution. Note that 
torsion of a circular shaft with diametrically flattened sides, pre- 
viously published,! follows as the special case of this investiga- 
tion. Table 1 presents the values of stiffness and maximal 
shear stress for different shape ratios S = b/a = h/a. It is 
easily shown that the maximum shear stress occurs at the ex- 
tremities of the minor axis of the ellipse. 
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On the Effects of Flow Obstructions on Free- 
Convection Boundary-Layer Oscillations 


J. P. Holman, K. E. Stout, and E. E. Soehngen 
Fluid Dynamics Research Branch, Aeronautical Research Laboratory, 
Wright-Patterson Air Force Base, Ohio 


January 19, 1960 


| ese interferometric studies of the detailed nature of arti- 
ficially induced oscillations in free-convection boundary 
layers! generated interest in the effects which flow obstructions 
might have in stabilizing the induced oscillations. Consequently, 
three types of flow obstructions were introduced into the oscillat- 
ing boundary layer, and their respective proficiencies in damping 
the oscillations were investigated. 

A controlled oscillating free-convection boundary layer was 
produced on a vertical, isothermal flat plate by pulsing with an 
electrical signal a wire 0.001-in. in diameter inserted in the bound- 
ary layer. The details of the pulsing mechanism are described 
in reference 1. The flow obstructions were subsequently intro- 
duced into the boundary layer at various distances from the 
heated plate, and their damping effects on the amplitude of the 
wave were measured. The amplitude of the wave was defined as 


A = (Ymaz — Ymin)/2 (1) 


where Ymaz is the distance from the wall to the maximum point 
in a particular isotherm and y,,;, is the distance from the wall to 
the minimum point in the isotherm viewed in the interferograms. 
The amplitude is said to be evaluated at the mean distance 
= (Ymez + Yuin)/2 (2) 
Some interferograms of the waves before and after insertion of the 
obstructions are shown in Figs. 1, 2, and 3, and the data for an 
x-position of 13 inches are given in Table 1. 
As indicated in the figures, all three types of flow obstructions 
caused a reduction in the amplitude of the boundary-layer waves. 
The curved-surface obstruction is most successful in stabilizing 


Fic. 1. Effect of wedge-shaped obstruction. (Left) obstruction 
0.26 in. from plate. (Right) obstruction out of boundary layer. 
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Fic. 2. Effect of curved-surface obstruction. (Left) obstruction 
0.41 in. from plate. (Right) obstruction out of boundary layer. 


the oscillations due to the viscous action and channeling of the 
flow. The flat surface is equally as proficient, but it has the dis- 
advantage that it cuts off part of the flow. The wedge-shaped ob- 
struction does not stabilize the oscillations too well until it has 


| 
Fic. 3. Effect of flat-surface obstruction. (Left) obstruction 0.29 
in. from plate. (Right) obstruction out of boundary layer. 
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TABLE 1 
Amplitude Variations at X¥ = 13 In. for 1.45 cps 
Induced Oscillation 


Distance of | Wedge Curved Surface} Flat Surface 
Obstruction 
from Plate A Vin Ym A 
(in. ) (in. ) (in. ) (in. ) (in. ) (in. ) (in. ) 
Free stream 0.145 0.060 | 0 140 0.050 | 0.140 0.055 


0.213 0.103 | 0.200 0.080 | 0.215 0.108 
0.305 0.145 | 0.270 0.115 | 0.300 0.140 
0.430 0.210 | 0.358 0.163 | 0.430 0.210 
Wedge 0.39 | 0.198 0.053 | 0.128 0.023 | 0.140 0.020 
Curved 0.45 | 0.283 0.078 | 0.185 0.030 | 0.220 0.040 
Flat 0.49 | 0.390 0.110 | 0.250 0.035 | 0.315 0.065 
0.575 0.125 | 0.310 0.045 | 0.425 0.085 
Wedge 0.26 | 0.183 0 058 | 0.145 0.015 | 0.130 0.010 
Curved 0.41 | 0.263 0.083 | 0.205 0.025 | 0.200 0.020 
Flat 0.39 | 0.380 0.115 | 0.270 0.0385 | 0.285 0.035 
0.610 0.140 | 0.340 0.040 | 0.380 0.050 
Wedge 0.21 | 0.165 0.040 | 0.130 0 0.100 0 

Curved 0.21 | 0.243 0.063 | 0.180 0 0.160 0 

Flat 0.29 |} 0.358 0.103 | 0.235 0 0.230 0 

0.300 0 0.320 0 


been inserted into the boundary layer to the point where it is re- 
moving a significant portion of the flow. 

Calculations of heat transfer by the method of reference 2 in- 
dicated that the time-average value o¥ the heat-transfer coefficient 
for both the undamped wave and for the damped wave was al- 
most identical with the calculated steady-state value. Details of 
the experimental techniques used and an analysis of the specific 
behavior of the oscillating boundary layer under the influence of 
the flow obstructions is given in reference 3. Motion pictures of 
the oscillatory boundary-layer flow are available from the first 
author. 
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Use of Reference Enthalpy in Specifying the 
Laminar Heat-Transfer Distribution Around 
Blunt Bodies in Dissociated Air 


E. R. G. Eckert and O. E. Tewfik 


Heat Transfer Laboratory, Institute of Technology, 
University of Minnesota, Minneapolis, Minn. 


December 24, 1959 


SYMBOLS 
Cp = specific heat at constant pressure 
h = heat-transfer coefficient 
i = enthalpy 
L = characteristic body dimension 
q = heat flux per unit area 
q(0) = heat flux per unit area at the forward stagnation point 
(ro) = radius of cross section of body of revolution 
s = coordinate along the body surface measured from the forward 
stagnation point 
u = velocity 
e = density 
“ = viscosity 
k = coefficient of thermal conductivity 
Pr = Cpu/k = Prandtl number 
Nu = hL/k = Nusselt number 
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Subscripts 
e = conditions at outer edge of the boundary layer 
) = at the forward stagnation point at edge of the boundary layer 
= local recovery conditions 
total 


i = wall conditions 
= undisturbed flow conditions 
* = conditions at reference enthalpy 


oa & 


I REFERENCE 1, a method to predict the distribution of the 
ratio of the local heat flux to the heat flux at the forward 
stagnation point of an axially symmetric, highly cooled, blunt- 
nosed body in dissociated air was developed. The concept of 
“local similarity”? was utilized to determine the local value of the 
enthalpy gradient at the surface. The method was applied to a 
specific example of an air flow of Mach number 12 over a hemi- 
sphere-cylinder body and a flat-nosed cylinder at zero angle of 
attack, and the results were compared with the predictions of 
Lees’ method.? For the hemisphere-cylinder body, there is 
little difference between the predictions of the two methods, 
in spite of some important differences in assumptions between 
them. For the flat-nosed cylinder, the differences are larger. 
In reference 1, Figs. 3 and 4, the theoretical predictions were also 
compared with the results of shock-tube experiments, and the 
agreement appears to be fair. The purpose of this note is to test 
the concept of the reference enthalpy for the prediction of the 
heat-transfer distribution around blunt bodies by a comparison 
with experimental results. Besides its demonstrated applica- 
bility for laminar and turbulent air flows with zero pressure gra- 
dient,? the reference-enthalpy method has recently been shown 
to predict quite closely the heat transfer to the forward stagnation 
point of a blunt body from a dissociated air stream; and, thus, 
it seems reasonable to expect it to work in specifying the heat- 
transfer distribution around a blunt body also. In searching for 
a suitable relation in which to evaluate the properties at a “‘refer- 
ence enthalpy,’’ Lees’ expression for the local heat flux 


gue = 0.50 (Pr)~ 2/3) popo)!/? «ig: F(s) (1) 
where 
( pette/ (tte / Uc 
(pette/ posto) / to is| 
0 
(n = QO) for two-dimensional flow, » = 1 for rotationally 


symmetric flow) 


is chosen,? because in reference 5 its predictions were shown to be 
closer to experimental results than other predictions.*.7 In 
deriving Eq. (1), the specific heat was assumed constant and (pp) 
invariable across the boundary layer. These assumptions are 
no longer realistic at very high free-stream stagnation temper- 
atures, and, thus, Eq. (1) will be modified in the following analysis 
by the use of the reference-enthalpy concept to account for such 
variations. First, Eq. (1) will be transformed into an expression 
for the Nusselt number as a function of Reynolds number and 
pertinent dimensionless fluid properties. 


Nuy = hL/ky = 0.35(He/pw)Reo, (2) 
(pe/ te/Uo )(ro/L)" 


s/L 
(pe/ po) He/Mo)(Ue/ )( ro/L)" a(s/t) | 


Re, o = Pua L/ wo (3) 
h= Gulp» w/(to tw) (4) 


1/2 


In the above equations, it is proposed to evaluate the properties 
at areference enthalpy 7* given by* 


i* = (ze + tw)/2 + 0.22 (i, — te) (5) 
Thus, in terms of reference conditions, Eq. (2) becomes: 


Nu* = Lh/k* = 0.35 Re, H*(s)-(Pr)!4 (6) 
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erence-enthalpy method with other methods and experiments for 
a hemisphere-cylinder body. 


where 


= (7) 
(ro/L)" d (s/L | 
0 


In addition, it is recommended to replace the enthalpy ip in Eq. 
(4) by the local recovery enthalpy and to define the heat-transfer 
coefficient 
h = Gulp, w/(tr — tw) (8) 
Now the ratio of the specific local heat flux to the stagnation- 
point heat flux may be obtained from Eq. (6) in the form: 
Guw/Gu(O) = — — (9) 


The heat-flux distribution was computed with Eqs. (7) and (9) 
for the two body shapes selected in reference 1. The results are 
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plotted in Figs. 1 and 2 as full triangles and as a dashed and 
dotted line, respectively, while the full and dashed lines show the 
predictions of references 1 and 2, respectively; and the points 
correspond to the experimental data as shown in Figs. 3 and 4 of 
reference 1. It is evident that the reference-enthalpy method 
predicts a heat-flux distribution which is close to the predictions 
of reference 1 or reference 2 for the hemisphere-cylinder. For the 
flat-nosed body, reference 2 predicts a distribution which for 
s/L > 1 tends to be higher than either reference 1 or the present 
method or the experimental points. 

Thus, Lees’ method modified with the reference-enthalpy con- 
cept appears to predict closely the local heat-transfer coefficients 
along the surface of a blunt-nosed body in dissociated air flow. 
The simplicity of the method recommends it for engineering 
calculations. 
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SYMBOLS 
a = velocity of sound 
C = constant for friction velocity 
= tw/(1/2), pe, local skin-friction coefficient 
M = local Mach number outside the boundary layer 
n = velocity-profile parameter: u/ue = (y/6)(/™ 
= static pressure 
R = gas constant 
T = temperature absolute 
u = velocity inside the boundary layer 
ue = free-stream velocity at the edge of boundary layer 
u* = friction velocity 
x = distance along the body surface from the stagnation point 
y = normal from the wall 
6 = momentum thickness of boundary layer 
6 = thickness of boundary layer 
y = ratio of specific heats 
Pr = Prandtl number 
= 2c 
p = density of air 
7 = shearing stress 
v = kinematic viscosity 
Subscripts 
é€ = condition and properties of free stream outside the boundary layer 


behind the shock 
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= compressible flow 

= incompressible flow 

conditions at stagnation point behind shock wave 

= conditions of free stream ahead of shock 

am = arithmetic mean based upon Tgm = (1/2)(Te + Tw) 
= conditions and properties at the wall surface 


F® THE analysis, the following assumptions are made: (1) 
Pr = 1, and the wall is insulated; and (2) flow outside the 
boundary layer is isentropic and varies linearly with 7. For air, 
this is adequate in the range of T = —250°F. to T = 3250°F. 


MOMENTUM THICKNESS OF COMPRESSIBLE FLOW 


Approximating the velocity profile by a power law, the in- 
compressible boundary-layer momentum thickness 


0; = [n/(n + 1)(m + 


Similar to incompressible flow, expressing the friction-velocity 
law of two-dimensional compressible flow by 


u/u* = 
constant C, may be determined by using! 
Tw/( pam Ue”) 0.013/ 


and the Stewartson-Illingworth transformations? * of incom- 
pressible to compressible flow and the reference temperature‘ as 


Cm = + — + M2)/10)" 
[1 + 0.032 M%+ 0.29(y — 1)M2]-@+n/2 


Where C; is known constant for incompressible flow. Then ap- 
plying this, the compressible momentum thickness @, which is 
applicable for a body of revolution, is 


@ = A[10/(10 + + [(-y — 1)/4]M?} 
(1) 


where 

A =n/(n + :1)(m + 2)C" 
under assumption (1), 

Pe/ Pao = 
and from the energy equation, 


Me? /U = [2/(y — Pa /PoM w*)[1 — (Pe/Po)( po/ 
(2) 


Under assumption (2), 


vo(Te/To)/(0/ pe) = vol 


Hence, 


Ve 


Ue 


Ve 


where 


vo 


Therefore, Eq. (1) reduces to 


10 (n+) /2 


(be/ bo) (n+1)/2 
[1 — 


(3) 


DIFFERENTIAL EQUATION OF COMPRESSIBLE TURBULENT SKIN- 
FRICTION COEFFICIENT AROUND AN ARBITRARY Bopy OF 
REVOLUTION 


If pressure or Mach number distributions along the surface are 
known, then ¢; may be determined from the following momentum 
equation: 


(d0/dx) + 0[((H + 2 — M*)(1/ue)(due/dx) + 
(1/r)(dr/dx)] = tw/pete? (4) 
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Since wu, and 6 are determined as functions of p, and M behind the shock and as functions of property values ahead of shock wave from 
Eq. (2) 
1 due Ue d(ue/u.) [A d/dx)( Po) 


ue dx Ue dx 241 — 


and from Eq. (3) 


d6 10 n+1 2 4 
dx + [( 2 ( ) 10 + = ) ( 25 ) 


Diy (: ur) + 9. x 


(pe/po)~A 
TL = 


— — (pe/po)™ 


/ 0 j ( / 
[1 — (pe/po)% dx \ po 9 dx 4 [1 — (pe/po) 


H; = 1 + (2/n) 
H +1 = (%/7T.)(Ai + 1) 


For incompressible flow, 

and under assumption (1)® 

Therefore, the first term in its bracket of Eq. (4) is 
H +2 — M? = (1/n){3n +2 + [(n + 1)(y — 1) — n| M2} 

Finally, the momentum equation (4) can be written in a differential form as follows: 


d\/dx = —filn, Pes + fo(n, per M)A~ 


where 


2 d([p y-1 y-1 n+l ( 
film, Pe, M) E (1+ 4 ‘cake 


(pe/ py) 


+ — 1)/4) M2} — y Ppe/Po 2y 1 — (pe/po) 7 


1 
{38n +2 + [(n+1)(y¥ —1)—- 


and fo(n, pe, M) = 


A-B 10 


[2/(n + 1)] + M? 


(1 — ld 
1 + r dx 


1+ [(y — \p, 


This equation may be expressed as a function of p/p» or M only because 


M = {[2/(y — 1)] — 
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wm a plane shock, generated by a uniform compressive 
piston motion, enters a tube of varying cross section, the 
shock is perturbed, the shock strength is altered, and the subse- 
quent flow is nonisentropic. There are three distinct contribu- 
tions to the perturbation, namely, a permanent disturbance due 
directly to the area change; a transient disturbance, propagated 


upstream with sonic velocity relative to the flow behind the 
shock, due to reflections from the shock of the permanent dis- 
turbance; and an additional transient disturbance, propagated 
downstream with sonic velocity relative to the flow behind the 
shock, due to reflections from the piston. 

By the use of a linearization based on small area variations, 
this problem has been solved for the case of a slowly converging 
or diverging one-dimensional channel and an expression given for 
the pressure change behind the shock! and the results utilized to 
discuss converging cylindrical and spherical shocks.? Perturba- 
tions of a piston-driven shock due to mass, momentum, and heat 
sources have been considered by Mirels.* 

In the present paper, the effect of an arbitrary area variation, 
applicable, for example, to wall boundary-layer effect, and small 
nonuniformities. of the piston velocity are considered. With 
E(x), the cross-sectional area of the tube, it is convenient to let 


E(x) = Ey + E(x) 


where E, is the original uniform cross-sectional area and E,(x) 


is the small arbitrary perturbation. 
At t = 0, a piston is started impulsively from rest and moves 


with velocity 


(1) 
air, 

as 7 

P- 
is 
)/2 


Up = Uz + 


where t# is constant and 7(t) is a small perturbation. Let the 
original shock velocity be w. 

Along the shock, the sound-velocity perturbation and the par- 
ticle-velocity perturbation can be written in terms of the shock- 
velocity perturbation. Utilization of the nonisentropic perturba- 
tion of a constant state yields the solution for the state behind 
the shock in terms of the shock-velocity perturbation. But, 
since the particle-velocity perturbation is known on the piston 
path—viz., n(1), the following functional equation for the deter- 
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mination of the shock-velocity perturbation, w(t), is obtained: 
[ki + ke — — M)] + 
[coksv-! — ko + + = 
n(t) — + m))~! E,[wt/(1 — + 
(1/2)u2[E((m — E,[wt/A + 
us| — Ey [wet] (1) 
where the rest state in front of the shock is denoted by the sub- 
script 0 and the perturbed flow behind the shock by the subscript 
2; u and ¢ are the fluid and sound velocities; w = w — wm; 
My = w/o; M, = m = v = — 1), with c, the 
specific heat at constant volume; y the adiabatic index; 


ky = (1 + /(y + 1) 
(y + = 2M, + m(1 — 
and ks = ley — 1)/w){ — M2)/(y + 10 — + m? — 2 + [2 y — 1 + 


The solution of Eq. (1) is (see reference 1, p. 577): 


k=0 


(1/2)u2.[Ee(m — A, 


+ — [Eg(m? — Ey — + 


where 
1 4(1 — M,?) 27 1) 2M, m(1 — | (1 + My~*) 
(=17" 2yMo Uy +10 —  y+1 yt1* J (2) 
+ 2M, mi(1 — Mo-?) 1 4(1 — M,?) Ayvy-1) 
— -- m? — 2 
y¥+1 y+1 2yMo + 11 — y—1+2Mo-*f 
Eq. (2) gives the shock-velocity perturbation due to an arbitrary a ial 
cross-sectional area perturbation £,(x) and an arbitrary piston- 
velocity perturbation 7(f). C = constant of proportionality in temperature-viscosity relation; 
The coefficient Q(k) is a monotonic decreasing function of the 
: s = shock-wave heig 
shock strength and varies from the value 0 for a shock of zero K = hypersonic similarity parameter; for a fiat plate, K = Ma (d3/dz) 
strength to a value (0.457)*/1.52 for a shock of infinite strength. L = characteristic length 
The series (2) converges rapidly, and a few terms will suffice M = Mach number 


to give w(t). 

Since the pressure perturbation directly behind the shock is 
linearly related to w(t), the pressure perturbation is a monotonic 
decreasing function of the shock strength. 

This work was supported by a grant from the National Science 
Foundation. 
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On the Use of the Boundary-Layer Equations in 
the Hypersonic Viscous-Layer Regime* 
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SYMBOLS 


All flow quantities have been nondimensionalized with respect to appro- 
priate free-stream quantities—i.e., p = p*/po*, where p* is absolute pres- 
sure. Lengths are nondimensionalized with respect to L—i.e., 6 = 6*/L, 
x = x*/L, etc. 


* This work was supported by the Dept. of the Navy, Bureau of Naval 
Weapons, under Contract NOrd 7386. 


= static pressure 

Ry = Reynoldsnumber, Rp = poU ol 

x = dimensionless distance along flat-plate surface 
y = dimensionless distance normal to plate surface 
6 = viscous-layer displacement thickness 

€ = density ratio across shock, po */pshock*. 


[For >>1,¢€~(y — 1)/(y + 1).] 


Subscripts 

w = refers to conditions at wall 

o = refers to free-stream conditions 

6 = refers te conditions at y~ 5 
Superscripts 

0 = refers to zero slip at wall 

s = refers to finite slip at wall 

* = dimensional quantity 


bag BOUNDARY-LAYER equations have proved quite successful 
in the analysis of hypersonic viscous-interaction phenomena, 
both in the strong and weak interaction regimes. However, by 
virtue of their derivation, they could not be expected to hold near 
a sharp leading edge. Nevertheless, some recent studies!~* 
indicate that these equations appear valid much closer to the 
leading edge than would normally be expected. For example, 
Laurmann! shows that the use of the classical methods of bound- 
ary-layer analysis in computing the induced pressure on a flat 
plate agree with experimental data into the region in which the 
effects of slip would be expected to cause a significant difference. 
It is the purpose of this note to discuss this apparent extension in 
the validity of the boundary-layer equations. 

The viscous-layer regime‘ is defined as being slightly forward 
of the region in which the boundary-layer equations are valid. 
Since the boundary-layer equations are valid‘ if (1/e)6?< 1, 
the viscous-layer regime will be expected to include (1/«)é? < 1, 
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but not negligible. Furthermore, the regime is assumed to con- 
sist of a viscous region between a thin shock wave and the wall, 
but with small viscous stresses and conductive heat transfer near 
the shock. This regime may be analyzed by dividing the flow 
into an outer, nearly inviscid portion dominated by the char- 
acter of the shock wave, and an inner, viscous portion dependent 
upon the conditions behind the shock wave and at the wall. 
The shock wave is, of course, influenced by the mass displacement 
of the inner portion. The flow is further restricted to small 
shock-wave angles and K? > 1. 

The assumption of small shear stresses and conductive heat 
transfer near the shock permit an order-of-magnitude calculation 
of the flow properties behind the shock, assuming the classical 
shock conditions and adiabatic flow. Such a calculation will 
give the orders of the individual terms of the Navier-Stokes equa- 
tions in terms of the shock-wave slope and its derivatives. Since 
the region of interest is adjacent to the boundary-layer region, 
the shock-wave shape in the viscous-layer regime should be simi- 
lar to that of the boundary-layer region. The shock shape for a 
two-dimensional flat plate with a sharp leading edge is given by 
various boundary-layer analyses! as h, ~ x*/*, Also, in the vis- 
cous-layer regime, the shock-wave height and displacement thick- 
ness differ by the order of 6/K?, so it is assumed that in this re- 
gime 6 ~ x*/4, 

The results of the above analysis indicate that near the shock 
wave in the viscous-layer regime, the flow can be considered 
inviscid if 


CM.%/Ri< 1 


where L is a characteristic length defined as the distance from the 
leading edge to some position near the origin of the strong- 
interaction boundary-layer regime. The thickness of this outer 
portion can be no greater than the difference between the shock 
height and the displacement thickness—i.e. (Ay)o ~ 6/K?, as 
noted previously. The maximum change in pressure across this 
outer portion is found to be of order (Ap)o ~ 1/e. It is of in- 
terest to note that the outer-layer thickness decreases quite 
rapidly with an increase in AK, while the pressure change across 
the outer layer remains of the same order. The error involved 
in neglecting this pressure change is of order (A%e)~!. 

The region in which the shear stresses and heat conductivity 
can no longer be considered small will extend from the wall to, 
at least, the displacement thickness. The usual no-slip, insu- 
iated-wall boundary conditions and the orders of magnitude of 
the terms in the outer portion can be used to determine the orders 
of magnitude of the terms of the Navier-Stokes equations in the 
inner portion. The requirement that the inertia and viscous 
terms be of the same order results in 


52/e~ WC/RiMow (1) 


This is the same relation obtained for the strong-interaction 
boundary-layer regime. Furthermore, the above calculations 
indicate that the change in pressure across the inner viscous por- 
tion is of order 


(pu ps)/ps ~ V/C/RiM ~ &/e (2) 


Therefore, in the viscous-layer regime, the pressure at the wail 
without slip can be significantly greater than the pressure in the 
layer behind the shock. It must be remembered that this in- 
crease in pressure results from terms neglected in the usual 
boundary-layer analysis. The effect of slip has not been taken 
into account. 

The effect of slip on the wall pressure is given by Patterson’. 
The order of magnitude of the pressure change due to slip is 


(pe — po ~ —WC/RiMe ~ —(8/e) (3) 


The ratio of the pressure change across the inner portion of the 
viscous-layer regime to the pressure change at the wall due to slip 
is of order 


(pu bs) /(pe Pu) ~—]1 (4) 


This equation indicates that even for 6?/e < 1, but not negligible, 
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Fic. 1. The boundary- and viscous-layer regions on a flat plate 


in hypersonic flow. 


the pressure measured at the wa'l will be essentiaily the same as 
the pressure at the outer edge of the inner portion because of the 
opposing effects of slip and of the higher-order terms of the flow 
equations. For this reason, the calculation of wall pressure 
using the boundary-layer equations, without slip, will appear to 
be correct in the region in which the boundary-layer equations 
are becoming invalid. 

Eq. (4) also means that the wall pressure cannot be used to 
determine the location at which slip flow becomes significant. 
For example, the experimental data presented by Laurmann! 
would suggest that slip becomes important at 


~ 1 
[in the weak interaction regime (62/e) ~ V 
ever, according to Eq. (4), the slip may become of signif- 


icance at a much lower value of 7.2/R, but is masked by; the 
higher-order terms neglected in the boundary-layer analysis. 


How- 
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| TWO-DIMENSIONAL, steady flow there is a direct analogy 
between the magnetic vector potential A and the tempera- 
ture field 7, when the electric field in the flow is zero everywhere 
From Maxwell’s current equation and the generalized Ohm’s 
law for scalar conductivity, the following equation for the mag- 
netic field B is obtained: 

Vv X B = uo(u X B) (1) 
where uw is the velocity, u is the permeability, and ¢ is the electrical 
conductivity. 

For two-dimensional flow a vector magnetic potential or mag- 
netic field ‘‘stream function” is defined as foilows: 
B=V xX Ak (2) 


where k is the unit vector in the direction on which there is no 
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dependence. Eq. (1) can then be written in the form: 
V2A — po(u-)A = 0 (3) 

Replacing A by the temperature J and the magnetic viscosity 
(ua)! by the thermal diffusivity a = k/pc, we have the differ- 
ential equation for heat conduction in a moving fluid. 

The velocity field u must also satisfy a momentum equation 
which in magnetohydrodynamics contains the Lorentz force 
term. The effect of the magnetic field on the velocity is given 
by the magnitude of the magnetic pressure number 

Pp = (Bo?/u)/pU? 

If this is negligibly small, the flow velocity is unaffected by the 
presence of the magnetic field. On the other hand, the effect of 
the velocity on the magnetic field depends on the magnetic Reyn- 
olds number R, = woUL, which may be varied independently 
of Pg. An analogous type of reasoning is made in the heat- 
transfer problem in assuming that the fluid velocity is not affected 
by the temperature distribution. 

We next consider the boundary condition requirements for 
adjoining media in the two cases. In electomagnetic problems, 
the component of magnetic flux normal to an interface must be 
continuous, 


Bn = Ba 
Or, in terms of A, 
(0A = (0A (4) 


where 0/dt denotes the derivative tangent to the interface. Eq. 
(4) is equivalent to stating that A is continuous across a surface. 
The other boundary condition is 


where J, is the sheet current along the interface. In terms of the 
magnetic vector potential, this becomes: 


(0A — = Js (5) 


Similarly, for a temperature field, the temperature must be 
continuous, and 


k\(OT/dn), — kx OT/On)2 = Qs (6) 


where Q, is the surface heat generation and k is the thermal con- 
ductivity. Finally, within a solid placed in the flow, 


V2A = —yj: (7) 
= —g/k (8) 


where j; is the current density and qg is the heat generation per 
unit volume. Thus, if 


= 
Fe = BQs (9) 
je = Bq 


the analogy between the magnetic vector potential field and the 
temperature field is complete, where 8 is the constant of propor- 
tionality. 


APPLICATIONS 
A few problems in heat transfer whose solutions can be inter- 
preted immediately as those for the corresponding magneto- 
hydrodynamic problems are next discussed. 
The solution for a line source of heat perpendicular to a uni- 
form stream is given by:! 


T = (q/2mk)e U/e)r] (10) 


where Ko is the modified Bessel function of the second kind of 
zero order, 7 is the distance from the line source, and q is the heat 
emitted from the line per unit length per unit time. In the 
analogous magnetohydrodynamic problem,? we consider a line 
wire carrying a total current J. Then, Eq. (6) gives the solution 
for the magnetic field if g/k is replaced by yl, and T by A. 

As another example, consider the solution for a doublet in the 
x direction: 
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A = — 
(x/r)Ki[(1/2)uoUr]} (11) 
where /’ is the limit of uJd as the distance d in the x direction, 
between two wires carrying currents J and ( —/), tends to zero. 
Finally, consider an infinite strip —© < zs < © in the plane 
y = 0. Current J, per unit width flows in the z direction in the 
strip, and the surrounding fluid moves with velocity U in the 
direction of the x axis. The solution of the heat-transfer analog 
is discussed in Carslaw and Jaeger.! If we further make the 
assumption that the magnetic Reynolds number Rm is large, 
we may make the usual boundary-layer assumption that the 
term 0?B/dx? is negligible compared with the other terms in 
Eq. (3). Thus, the governing differential equation is 
/dOy? = po U(0A/Ox) (12) 
For a semi-infinite flat current sheet beginning at the origin, A 


for x > 0 can be written down from the corresponding transient 
heat-transfer solution on p. 75 of Carslaw and Jaeger :! 


A = exp [—y?/(4x/uoU)] — 

(1/2) |y| erfe (|y|/2Wx/yoU)} (13) 
Other heat-transfer problems and solutions involving heat sources 
may be interpreted as analogous magnetohydrodynamic prob- 
lems. For example, the flow normal to a rod with a uniform 
current within it corresponds to the flow over a rod with uniform 
heat generation within it. 
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Steady-State Melting of a Semi-I!nfinite 
Medium With Temperature-Dependent 
Properties 


Stephen J. Citron 
Assistant Professor, Division of Engineering Sciences, 
Purdue University, Lafayette, Ind. 


January 25, 1960 


ie PURPOSE of this work is to show that solutions for the 
steady-state temperature distribution, the steady-state 
melting rate and the amount of material melted and ablated 
from a semi-infinite medium subjected to a heat flux at the melt- 
ing face can be obtained when the medium has temperature- 
dependent thermal conductivity and specific heat. 

The present work will be concerned with the melting of a semi- 
infinite body acted on by a heat flux Q at the melting face under 
the assumption of instantaneous removal of the melt formed 
or the alternative assumption that the liquid layer does not 
appreciably reduce the heat input to the solid.!_ We shall allow 
the thermal conductivity and the specific heat to be temperature- 
dependent and assume that they ‘may be represented by power 
series of the form 


MN 
RT) = >> T/T*)" 
n=0,1 
(1) 
C(T) = >> 


n=0, 1 


where 7™* is the melting temperature of the material. 

If the position of the melting face is denoted by x = s(t) where 
x = Ois the position of the face at the start of melting, the deter- 
mination of the location of this moving boundary and the deter- 
mination of the temperature distribution within the solid requires 
the solution of the heat-conduction equation 


(0/02) [k(T)(OT/dz)] + pC(T)8OT/ds) = pC(T)OT/dr) (2) 
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TABLE 1 
Dimensionless Steady-State Melting Rates 
pLs;/ Q 
pLs;/Q Temperature Independent 

Temperature Dependent cc*) (Cr) 
Beryllium 0.239 0.196 0.303 
Copper 0.307 0.280 0.338 
Graphite 0.869 0.846 0.927 


given here in terms of the independent variables 2 and +r where 
z= x — s(t), r This equation is to be solved in the region 
0<2< © under conditions 


Tl. = T* (3a) 

* Tp (3b) 

T(z, r*) = T(z) (3c) 

Qir) = + (3d) 


where 7* is the time at which melting begins and L is the latent 
heat of fusion. 

To investigate steady-state melting, we set 07/07 = O in 
Eq. (2) and let Q(7) = Q, a constant, obtaining after making use 
of relations (1) 


Ni N 
d/dz | | + ps C,(1/T*)"dT/dz = 0 (4) 
0 0 


This equation may be integrated directly, yielding 


N 
(T\_ 


Determining the value of k(7)(d7/dz) at s = 0 from this equa- 
tion and substituting into Eq. (3d), we obtain for the steady- 
state melting rate (S;) in dimensionless form 


Ls 1 
(6) 


When the specific heat is constant, this reduces to 


pLi; 1 7) 
Q {1+ 

Table 1 was constructed to illustrate the differences which may 
occur in the value of the steady-state melting rate computed 
assuming the specific heat to be temperature-independent and 
computed using the actual temperature dependence. The vari- 
ation of the specific heat with temperature and the value of the 
latent heat of fusion were obtained from references 2-4. In the 
case of beryllium and graphite, a parabolic representation was 
used for specific heat vs. temperature while for copper a linear 
representation sufficed. The dimensionless steady-state melting 
rate for each material was then calculated from the variable- 
property solution of Eq. (6) and compared with that calculated 
from the constant-property solution of Eq. (7). In the constant- 
property solution, for comparison purposes, the melting rate was 
calculated for each material using two values of the specific heat: 
first the value at the melting temperature (C*) and then the value 
at room temperature (Cr). It is seen that the difference be- 
tween the constant- and variable-property solutions for the melt- 
ing rate can become appreciable. 

Having found the steady-state melting rate, we now determine 
the temperature distribution existing when steady-state condi- 
tions have been achieved. The determination of the steady- 
state temperature distribution 7;(z) requires the integration of 
Eq. (5) which may be written as 


41 n 
(Z) 


| 


1 T/T* 


9 #t+1 
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The quantity to be integrated is a rational function of 7/7* and, 
hence, the integration may always be carried out by the method 
of partial fractions’ and the temperature determined implicitly 
as a function of z. In the case of constant properties, Eq. (8) 
yields 


Ths) — T. = (T* — 


The amount of material melted as a function of time remains 
to be determined when the properties are temperature-dependent. 
We proceed by establishing bounds on the values of s(r) in a 
manner analogous to that carried out in® for the case of tempera- 
ture-independent properties. Integrating eq. (3d) with respect 
to 7, we obtain 


Tr oT 
Q(r — r*) f —k(T*) 3! dr + pLs(r) (9) 


* 
Substituting for 
—k( T*)(0T/dz) 
the expression 


—k(T*)(OT/dz)| 


oT oT 
pC(T) — dz + pC(T) — dz (10) 
0 Oz 0 Or 


obtained by integrating Eq. (2) with respect to z from 0 to © and 
using Eq. (3b) yields 


where, since 7, is a constant, the second term on the right may 
be rewritten as 


T Pa) 
ff eC(T) — (T — T,,)dedr = 
ade dr 
f eC(T\(T — T. Mas 
0 


ff (12) 


Making use of Egs. (1) and (12), Eq. (11) becomes upon direct 
integration of the first term on the right 


N 

Q(r — r*) + 


n+l 
[: J eC(T)(T(2, 7) — — 
0 


f — 
0 


a 
r*J0 


Substituting from Eq. (6) and equating the second integral on the 
right with the amount of heat absorbed by the body in the in- 
terval 0 < +r < before melting begins, we obtain 


s(r) = 1" a) pC(T)[T(s, r) — + 
2 Jo 


Eq. (13) differs essentially from Eq. (2.8) of reference 6 by the 
addition of the last term on the right which must be present if 
the specific heat is to be temperature-dependent. As T7(z, r) 
is nondecreasing and the specific heats of materials in general 
increase with temperature, it is seen that this term is positive- 
definite, that is 
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1 oc(T 
J pIT(z, r) — To] dzdr = M%r) >0 
Or 


Differentiating Eq. (14) with respect to 7, we obtain 
J, : C(T) ) T|d: (16) 
Mr) = 8 Zz, — T.,|dz 
f 


from which it can be seen that the largest melting rate possible 
is the steady-state one since C(7) and 07/07 are both positive 
quantities. 
As T(z, 7) increases with time and is positive-definite, 
we see that the quantity 
&(7r) = — + (17) 


obtained by using 7)(z) in place of 7(z, 7) in the first integral of 
Eq. (14) is greater than s2(7) = &;(r — +*), which is still greater 
than s(7) since ; is the largest melting rate possible. 


= &(r — > (18) 


Further, the quantity 
= 47 eC(T,)[T(z) — Ta (19) 
\ Jo 


obtained by using 7;(z) in place of T(z, 7) in the first integral of 
Eq. (14) and neglecting 1/2(7) will always be less than s(7) since 
M*%7r) > 0. We have therefore established the bounds s2(7) > 
s(r) > si(7) on the amount of material melted. |The difference 
between these bounds 


(Jo 


f Tae) ~ (20) 
0 


will be small if 7;(s) does not differ appreciably from 7)(z) as is 
often the case in the solution with constant properties.’ In 
practice, since steady-state conditions are achieved rapidly,* the 
upper bound will furnish a good approximation to the amount of 
material melted. 
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A Phenomenon in Unsteady 
Magnetohydrodynamic Flow Theory 


S. F. Borg 


Professor and Head of Civil Engineering, 
Stevens Institute of Technology, Hoboken, N.J. 


January 21, 1960 


SYMBOLS 


= applied magnetic-field vector 
= resulting magnetic-field vector, including induced effects 


= pressure 
v — (r/t), reduced velocity 


(0/0E, 0/07, 
= (0/dx, 0/dy, 0/02) 


r = radius vector (x, y, 2) 

t = time 

Vv = velocity vector 

7 = permeability 

(€,”,¢) = r/t, the similarity parameter, (x//, y/t, 


HE MAGNETOHYDRODYNAMIC fluid equations for a nonviscous, 
incompressible, non-heat-conducting and nonradiating liquid 
with very high conductivity are 


DB/Dt = (B-V)V (1) 
p(DV/Dt) + Vp = (1/u)(V X B) X B (2) 
V-V=0 (3) 


In the above, D/Dt is the substantial differential operator. In 
addition to the above, the fluid field must satisfy the Maxwell 
equations, but these can be solved independently after V, B, and 
p have been determined from the above. 

The B in Eqs. (1) and (2) is the resulting magnetic-field vector 
due to the disturbed action of the fluid which had an initial, 
applied field By. In other words, B includes the effects of the 
induced field. 

A problem of interest in several fields is that of a wedge (two- or 
three-dimensional, infinitely long, straight-sided) impacting 
with a constant velocity on a semi-infinite body of perfectly con- 
ducting fluid. A fundamental simplification is introduced in the 
solution of these problems, if they areformulated in terms of the 
similarity parameter, 


(& ¢) = r/t (4) 


If this substitution is made and if, in addition, the substitu- 
tion! 


Q = V —-(r/t) (5) 

is introduced, the equations (1), (2), and (3) take the form 
2B = V: X (CO X B) (la) 
+ + Vip = (1/u)(Vi X B) X B (2a) 
(3a) 


The quantity Q has the dimensions of a velocity; therefore, one 
may, in analogy to streamlines, define Q lines as the lines which 
are everywhere tangent to the direction of Q in the (&, 7, ¢) plane. 
In reference 1, it is shown that the Q lines are particle paths of 
the fluid in the (£, 7, ¢) plane and also that the Q lines extend to 
infinity from the neighborhood of the wedge. 

In Egs. (la) and (2a), it is assumed that B is either a constant 
or else B = B(E, n, £). 

Physically, the Eqs. (la), (2a), and (3a) mean that the time- 
dependent impacting phenomenon, which is represented by 
different maps, for different times, /, in the x, y, z plane, is repre- 
sented by a single map which holds for all times in the (é, 7, ¢) 
plane. 

We now specialize the problem to the case of a two-dimen- 
sional wedge (conditions the same in all planes perpendicular to 
thesaxis.) Also, we assume that 

By, = Bak 


where B, is constant throughout the field. 

In the present two-dimensional case, where the applied mag- 
netic field is normal to the plane of the flow, we may also take 
the induced and, hence, final field to be normal to the plane of 
fluid flow. Also, we may assume that the disturbances are not 
transmitted to infinity. 

Then, since V = V(x, y), Eq. (1) becomes 


DB/Dt = 0 (1b) 
which, in terms of Q and Vj, is given by 
Q-V.B = 0 (1c) 


that is, B is constant along Q lines. 
Because B is constant along Q lines and because the Q lines 


exte 
foll« 


forir 
cons 


Meck 


and 
the si 
gree ¢ 
Th 
in aer 
axis 1 
matec 
in ait 
used | 
earth 
force- 
or six 
permi 
The 
with t 
y axis 
the z 
The 
been o 


It 
side 
stan 
pres 
wha 
pres 
the | 
theo 
Fy 
1B 
An 
Fre 
Willi 
Aero 
i Bell 
Febri 
Lag 
m 
g 
é 
D 
L 
Yw 
H Pw 
“ k 
Ss 
h 
6 
Tz, 
I x2 
a 
B 
L 
M 
N 
A | 
= | 


us, 
1id 


READERS’ 


extend to infinity (where B = By, the constant applied field), it 
follows that B = constant everywhere. 

In other words, for the particular case of (a) an infinite, straight- 
sided, two-dimensional (x, y plane) wedge, (b) moving with con- 
stant vertical velocity (c) impacting on a semi-infinite incom- 
pressible perfect conductor (d) with applied magnetic field 
B, = Buk, where By, = constant, the magnetic field has no effect 
whatever upon the dynamic behavior of the wedge—that is, the 
pressure on and velocity of the wedge may be determined from 
the usual conservation equations of incompressible hydrodynamic 
theory, without regard for magnetic effects. 

From a similar analysis of the steady-state two-dimensional 
form of Eqs. (1)-(3), we may conclude the magnetic field B is 
constant along each streamline for these flows. 
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An Axis System for Five and Six Degree of 
Freedom Airplane Dynamic Problems 


William L. Wynn 


Aerospace Department, 
Bell Aircraft Corporation, Buffalo, N.Y. 


February 1, 1960 


SYMBOLS 
V = total velocity of airplane c.g. 
m = airplane mass 
g = acceleration of gravity 
= vertical flight-path angle 
= azimuth flight-path then pitch system 
D = aerodynamic drag 
L = aerodynamic lift 
Yw = aerodynamic side force perpendicular to the velocity vector 
Py = angle through which the side force, Vw, would have to rotate 
about the velocity vector to be parallel to the earth 
R = down range 
Zz = side distance 
h = altitude 
¢ = rate of roll about a body axis 
6 = rate of pitch about a body axis 
¥ = rate of yaw about a body axis 
Iz, ly, 1z = body-axis inertias 
Ixz = product of inertia 
a = aerodynamic angle of poe, then 
B = aerodynamic sideslip angle 
z = aerodynamic rolling moment 
M = aerodynamic pitching moment 
N = aerodynamic yawing moment 


A COMBINATION of an earth-wind force-equation axis system 
and a body moment-equation axis system is suggested as 
the simplest airplane representation to use for five or six de- 
gree of freedom airplane dynamic problems. 

The six degree of freedom airplane equations have appeared 
in aeronautical literature and textbooks in various forms. Body- 
axis moment equations are invariably used (or at least approxi- 
mated) so that the airplane inertias do not change with changes 
in airplane orientation. However, various axis schemes are 
used for the force equations, such as wind axes, stability axes, 
earth axes, as well as body axes.!. This note suggests yet another 
force-axis system that appears to offer some advantages for five 
or six degree of freedom airplane dynamic problems that do not 
permit the use of linearized equations. 

The force equations are written in an earth-wind axis system 
with the x axis along the velocity vector of the airplane c.g., the 
y axis perpendicular to the x axis and parallel to the earth and 
the z axis mutually perpendicular to form a righthanded system. 

The force equations of motion are thus (thrust terms have 
been omitted for simplicity) 
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Fic. 1. Equation schematic. 


mV = —D — mg sin y (1) 
mVy = Leos — Yy,sin — mg cos y (2) 
cos y = Y,, cos + L sin (3) 


where ©, is the angle through which the side force, Y,, would 
have to rotate about the velocity vector to be parallel to the 
earth. Likewise, earth distances are given by 


R = Vcos y cosé (4) 
8 = Vecos ysiné (5) 
h = Vsiny (6) 
Writing the moment equations in the usual way 
— + — — 166 = L (7) 
1,6 + Uz — — — = M (8) 
— + (I, — 12)66 + = N (9) 


relations are required for a, 8, and ®, as functions of the flight 
path and body motions as determined from the force and mo- 
ments equations above. 

The necessary kinematic equations are: 
a= 6- 

(¢ cos a + Wsin a) sin B + ¥ cos &, + Ecos y sin &, (10) 


cos B 


= dsina — ycosa+écos y cos — ysin & (11) 


$, = (¢dcos a + y sin a) cos B + (6 — &) sin B + é sin y 
(12) 

The use of the equations in a five or six degree of freedom simula- 
tion is illustrated schematically in Fig. 1. For a five degree of 
freedom simulation, Eqs. (1) and (6) are ignored, and, hence, 
the altitude and velocity are (assumed) constant. As illustrated 
in Fig. 1, the desired feedback quantities, such as Euler angles, 
distances, etc., are computed as required and fed back to the 
pilot or autopilot as the case may be. The aerodynamic forces 
may, of course, be originally represented in a body-axis system 
and then converted to the wind-axis system before inserting in 
(1)-(3), if desired. 

The relative advantages of using an earth-wind force-axis 
system, together with the required kinematic relationships, will, 
of course, depend on the particular problem at hand. However, 
the following general advantages appear to be important: 

(1) The extremely long trigonometric terms, or the complica- 
tion of directions cosines, usually necessary for the gravity terms 
in the force equations and to compute earth distances, are at 
least partially avoided; the use of Eqs. (10)-(12) instead resulted 
in a significant equipment saving on one five degree of freedom 
simulation problem at Bell Aircraft Corporation. 

(2) The airplane representation is independent of the gyro 
orientation used for the pilot display or for autopilot control. 

(3) The simplifications that can be obtained if certain angles 
can be assumed small or if certain motions are assumed to be 
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zero are almost immediately clear with a minimum of further 
algebraic manipulation. 

(4) The force and distance equations are written in a form 
similar to that usually used for performance work; and, hence, 
results and checks can be compared more easily than if, say, the 
force equations were written in the body-axis system. 
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Comments on ‘‘Matrix Solution of the 
n-Section Column”’ 


M. Gololobov 
Zelena 2, Prague, Czechoslovakia 
January 27, 1960 


I SHOULD LIKE to draw attention to the fact that the method 
given in the article of Professor W. T. Thomson* is not valid 
for any end conditions as he states in the summary. 

Thomson did not take into consideration the shear force. 
His solution is, therefore, valid only for the case of a column with 
simply supported ends or a symmetrical column with both 
ends built-in. It is not valid for an unsymmetrical column with 
built-in ends, or a columin with one end built-in. In these cases, 
there is a shear force present as may be seen from Figs. 1 and 2. 
These two last cases do not Jend themselves to the solution 
given. Eq. (14) is, therefore, in error. 


* Thomson, W. T., Matrix Solution of the n-Section Column, Journal of 
the Aeronautical Sciences, Vol. 16, No. 10, pp. 623-624, October, 1949. 


Vi (1/8;) sin Bil; 

yi" 0 cos Bil; 

M, | | 0 sin Bil; 

Q 0 0 0 


and the overall matrix equation to 


Yn Ay Ay Ais Ars |] 
0 Ao Az, Aa |} yo’ 
M, 0 A 32 A 33 Az Mo 
Q 0 0 0 1 JjLe 
For the fixed-end columns, the boundary conditions lead to 
| Ais Aig 
= 0 
| Aes Ay | 


On a Solution to the Unsteady Laminar 
Boundary Layer 


H. A. Hassan 
Associate Professor of Aeronautical Engineering, Virginia 
Polytechnic Institute, Blacksburg, Va. 


Januaty 97, 1960 


ig IS SHOWN below that the transformation E = x/V 20, 7 = 
y/ V 2ut reduces the unsteady laminar boundary-layer equa- 
tions in two dimensions to an equation in which ¢ does not appear 
explicitly, provided that the free stream velocity V(x, ¢) can 


(1/E;1:8:2)(1 — cos B./;) 
sin Bil; 
cos Bil; 
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Author’s Reply 


William T. Thomson 


Professor of Engineering, University of California, 
Los Angeles, Calif. 
February 26, 1960 


Me GOLOLOBOy is correct in pointing out the error in Eq. (14) 
which holds only for the symmetric n-section column. 
For the case with shear Q, the section matrix should be modified to | 


— sin Bilis) |] yi-1 
— cos Bil;) 
(—1/8;) sin Bil; 

1 Q 


be expressed as (V v/2t) h(é). A power series solution is as- 
sumed for the resulting equation, and it is shown that for 


h(=) = ane" 


(a an integer) the solution can be expressed in terms of universal 
functions.!: 2 For certain values of a, the zeroth term of the 
assumed series, which satisfies the outer boundary condition ex- 
actly, is governed by well-known equations. 

The differential equation for unsteady, incompressible laminar 
boundary layer in two dimensions can be written as 


cy — yy = + Vo Veo. yyy (1) 
The boundary conditions are given by 
W(x, y, t) = ¥, (x, y, t) = ¥, (x, y, t) = 0 for y = 0 

and 

WV, (x, y, t) > V(x, as (2) 
Letting 
E = = y/V 20, t) = 

V(x, y, t) = 7) 


Eq. (1) reduces to 
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lasn— 


$a & + hh’ — (Eh)’ = 0 (4) The solution of Eq. (4), after introducing Eq. (6), will be 
and Eq. (2) gives assumed as 
¢=¢%.4 = ¢ = Oforn = 0 
and ¢n( (7) 


: . : in the region of convergence of 
In order that the solution can be expressed in terms of uni- 8 8 


-ersal functions, the following transformation will be introduced: a 
n= aft), = (h/f) 0) (6) n=0 

Carrying out the indicated transformation, using Eqs. (7) and 


where 
(8), and letting 


= for a < 0 


and 
in¢ 1 + (th’/h) = >» be", = 1+ @) 
0 


This transformation gives 
one finds for a < 0: 


+ (1/2)a0(1 + + aa(l ¢")=0 (9) 


and 
+ (1/2)al1 + — + + + (1/2)1 + Gn = Roi, forn > 1 (10) 
where 
= = an{(n + alg)” — 1) — [In + (1/2)1 + + 
n—1 n—-1 | 
ao > {(k + — [k + (1/2)1 + @)] — ge > + (1/2) 1 + + 
k=1 
n—1 n—1 n—k 
ak + + aal(m + + — [m+ k + (1/21 + forn } (11) 
b=1 k=1m=1 | 
and | 
n+a-—1 
= + — bie’ nt+a—1—1 — (n+ — — (1/2) 1+ ane’ forn > — a 
k=0 J 
For a = 1, one obtains 
+ We)’ — 1) + (aoe) + + al — ¢.”) = 0 (12) 
and + + en” + (mn + — + ain + = Roa, forn > 1 (13) 
where 
n—1 n—1 
= + 1)[ — 1 — gage”) + ao (+ — — + + 
k=1 k=1 
n—1 n—1 a—k 
> ak + — — b,( — 1) + adm + k + — (14) 
k=1 k=l k=1m=1 
And, for a > 1, one finds 
+ (1 + ale’ — 1) + = 0 (15) 
and 
+ (Lt at men’ + = Ria, (16) 

Eq. (9) reduces to the Blasius equation for a = 0, is integrable 
where in closed form for a = —1, and reduces to the Falkner-Skan 
Ris = Qn-1 = —bnlga’ — 1) — equation for the other values of a by letting ¢ = Bx, \ = Bn where 

n—1 
nS B = [2/a(1 + a)]"? 
k= ‘ ‘ 
; Eq. (12) has no known solution, while Eq. (15) can be solved in 
and | (17) closed form. However, to satisfy the boundary conditions, one 
TTi~—e a+i—e—k ’ has to restrict the valves of a in Eq. (15) to even integers. The 
Rn = Qn + Zz, > atk +a+m) X solution of Eq. (15) is given as 
k=0 k=0 
= 1+ (—1)/2)-! [2/2 (@/2)!/a!] X 
(—n2/2)] /dn®: ¢ 
«> d* {exp (—7?/2)]/dn*; a 
The boundary conditions for the functions ¢, are ; The functions ¢n, (n > 1) for all values of a are governed by 
linear equations and can be expressed as linear combinations of 


universal functions for a given value of a). This can be accom- 
plished by letting 


¢(0) = ¢,(0) =0, | (18) 
= gn'(0) = = 0, n >If 


; 
(2) 
(3) 
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= oft + agi + 
+ + + dege + + (20) 
+ Pi,1 + afin, 1 + 1; ete. 

where 6, = 1 for all m and B, = b,-1 for allu > 1. Using Eq. 
(20), the method employed in reference 1 can be used to find the 
governing differential equations for the functions f, g,... , ete. 
All these functions satisfy the same boundary conditions imposed 
on gn, (n > 1), namely, Eq. (18). 


$8 
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On Supersonic Flutter of Long Panels 


John W. Miles 
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Me JOHN HEDGEPETH has brought to our attention the recent 
work of Movchan! on supersonic flutter of flat panels of 
finite width, and, in particular, an apparent paradox regarding 
the stability criterion as the length (in the direction of flight) of 
the panel tends to infinity. We shall show that this paradox is 
a consequence of an improper application of the piston-theory 
approximation for the aerodynamic forces and that each of the 
two different results obtained by Movchan for an infinitely long 
panelisinerror. The same error arises also in Stepanov’s treat- 
ment? of cylindrical shells. 
Movchan starts from the differential equation (after setting 
N, = = kh = = Oin his more general formulation ) 


DV4w + phw (1) 
where w(x, y, ¢) denotes the transverse displacement of the panel, 
D = (1/12)E’h® = (1/12) [EA3/(1 — v?)] (2) 


its flexural stiffness, p its density, 4 its thickness, and g the aero- 
dynamic loading. Letting the panel move in the positive x 
direction with the velocity c and with both sides exposed to a gas 
of static pressure po, sonic speed co, and specific-heat ratio y, he 
poses the piston-theory approximation 


qd = (2ypo/co)(cwz — wr) (3) 


Assuming the panel to be simply supported at x = 0,@and y = 
0, b and choosing a solution to Eqs. (1) and (3) in the form 


w = X(x/a) sin (nry/b) exp (wt) (4) 


he deduces that the critical ‘‘preflutter velocity”’ is 


cn’ = X 

[5 + + (5) 
We remark that the corresponding instability is a pure divergence 
and that w = Oforc = c,’. Setting m = 1 and taking the limit 
a/b — ©, he concludes that the minimum flutter speed for an 
infinitely long panel is 


= a/b—> (6) 

Movchan then goes on to say (in the words of his translator) 
that: 

“On the other hand, if in the initial problem it is taken into 
account beforehand that the panel is infinitely long and no condi- 
tion is set down for the infinitely remote end of the panel, except 
the condition of arbitrary smallness of the initial perturbation, 
it is possible to prove for the panel the existence of flutter motions 
when the velocity exceeds the magnitude 


= rWVE'/3p(h/b) (7) 


The value [(7)] may be smaller than the value [(6)], which 
attests to the inapplicability of the formula [(7)] for unbounded 
panels. The mentioned example shows that results obtained by 
means of the examination of panels, cylinders, and so forth, 
theoretically infinite in the direction of the unperturbed motion, 
are by no means always applicable to cases of finite dimensions, 
even if these dimensions are sufficiently large.”’ 

Movchan does not give his derivation of Eq. (7), but if we as- 
sume the panel to be unbounded in the x direction, a sufficiently 
general motion (in virtue of Fourier’s theorem) is 


w = exp [ik(x + Vt)] sin (nzy/b) (8) 


Substituting this in Eqs. (1) and (3), we obtain the eigenvalue 
equation 


+ (nw/b)?]? — phk?V? = 2ik(ypo/co(c — V) (9) 


Eq. (9) yields unstable solutions if Im{ Vv} < 0, which condition 
leads to the critical speed 


c= V = WE'/l2p[k + (nx/b)*k—Jh (10) 


the minimum value of which (for 2 = 1 and k = 7/b) is given by 
Eq. (7). 

We now examine the assumptions implied by the piston-theory 
approximation. The linearized differential equation for the 
perturbation velocity potential is 


+ byy + = (1/c0*)[(0/dt) — (11) 


This leads to Eq. (3) if and only if the terms ¢,, and d,, on the 
left-hand side are neglected. Noting that dé; = Oatc = c,’, that 
= O(a~*d) and that dy, = O(b~*6), we conclude that the 
validity of such an approximation in conjunction with the solu- 
tion (4) in the neighborhood of c = c,’ demands 


(c/co)? > 1 + (na/b)? (12) 


It follows that Eq. (5) cannot be uniformly valid with respect to 
na/b and, hence, that Eq. (6) is not consistent with its ante- 
cedent hypotheses. Similarly, the validity of Eq. (3) in con- 
junction with the solution of Eq. (8) requires 


— V>coll + (13) 


which condition obviously is violated for c = V, as in Eq. (10). 
We add that Stepanov’s results for a cylindrical shell,? which 
differ from those obtained by the writer,‘ are in error for the same 
reason. 

More adequate aerodynamic theories for traveling wave mo- 
tions similar to that of Eq. (8) are discussed in references 3-5. 
The more basic question as to whether traveling wave instability 
based on linearized theory is significant for panels of finite length 
depends essentially on the wave lengths at which instability is 
predicted. A necessary condition for the preceding problem 
would be ka > 1. The question then would remain whether, 
in the presence of nonlinear effects, an unstable (on linearized 
theory) traveling wave would experience sufficient growth be- 
fore reaching the downstream support, where it presumably 
would be reflected to form a standing wave. Presently available 
theory provides no reliable answer to this question, and it seems 
likely that experiment must be ultimate arbiter. 
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On Theoretical Considerations of Plastic-Flow 
Criteria 


W. Komkov 

Assistant Professor of Mechanical Engineering, University of Utah, 
Salt Lake City, Utah 

January 22, 1960 


INTRODUCTION 


- MANY FIELDS of engineering, simplifications of certain general 
theory afford the engineers a useful tool for analysis of many 
practical problems. The linearization of Navier-Stokes equa- 
tions in hydrodynamic applications may serve as an example of 


ou, 1 ou, Ou2 \?2 ou; \? 
M= {—4+-=|(— = 
+(=)'} 
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such treatment. The value of having a general theory, even 
though in itself it is too difficult for practical use, is quite obvious 
since comparisons of the approximate equations, or approximate 
criteria with the most general case, give us an insight into the 
limitations of the approximate theory. 

This paper constitutes a step in the direction of introducing 
such general theory into the fields of elasticity and plasticity. 
The linear elasticity, including Hooke’s law and the various 
criteria for plastic flow behaviors, should then become simplifica- 
tions of this general theory. 


ELEMENTARY THEORY 
Let x1, x2, x3; be the coordinates of a point in an undeformed 
body and x, x2, x3; the coordinates of the same point after de- 
formation. Then, we introduce: 
Ou; Ou; Ole Ou; 
Ox2 Ox) Ox) OX2 Ox) OX; OX; 


Ox; Ox) Ox) Ox; Ox) Ox; 

1 Ou; Ou» 1 ou; Ou; Ou2 Ou» Ou; 

- -(— + — — — }, ete. 


For simplicity’s sake, let us restrict ourselves to a two-dimensional stress system, although as can be seen, the problem discussed for a 


three-dimensional case is only more difficult in the arithmetical sense. 


Ou, 1 Ou; \? \? 
M = — - — — 
1 1 Ou; Ou, 
2 (= 2 (= Ox; 


The matrix M is the strain matrix, which upon discarding all 
terms of order higher then the first would become: 


2 \Oxe Ox; OxX2 


In our discussion, however, we cannot ignore the higher-order 
terms, and 


2 4 Ox, Ox 
1 / Ou, Ous 1 Om 


du, \? \? 


CRITERIA OF PLAstic FLOW DEVELOPMENT 


M= 
(3) 


Ou2 
4 


Assumption: 
Let us assume that the stress tensor is some nonlinear function 
of the strain tensor and of its derivatives; for example, let us have 


TH = + + A22€22 + + €12 + buen? + 
duenen + +... + +...,ete.... (5) 
We have assumed here that stress is a function of strains and of 
time derivatives of strain, but not a function of position deriva- 


19 | fdm\? (2 : 1 
1 


on ft 12 4 (amy 
Ox? 2 OXe Oxe2 


tus, 2 Dus Dus 


Ox) 0X2 (OX, Ox; OXe2 


The matrix M becomes upon elimination of variables u; and x3 
Ox2 Ox, Ox2 OX) Ox; OX2 (2) 


dx, J" 21 \dxe Ox2 


tives of strain. For greater generality, the position derivatives 
should have been included at this stage. Assuming that the 
inertia forces are small compared to stresses, we could obtain the 
standard equilibrium equations: 
(O7n/0x1) + (Or712/0x2) = 
(6) 
(O712/0x1) + (O712/0x2) = 


Should the inertia forces be large, and should mass forces (X) 
be also included in the system, we would have 
0) 


(O711/0x1) + (Or12/Ox2) + — 
(7) 
of 


(O722/0x2) + (Or12/0x1) + p(O%u2/dt?) — 


PLastic-FLow CRITERION 

Substituting the values of 711, 712, 722 from Eq. (5) and writing 
the strains in terms of displacements according to Eq. (4), we ob- 
tain two nonlinear partial differential equations in terms of dis- 
placements “; and u:. (Note that it would make no difference if 
a three-dimensional stress system were considered. ) 

The criterion for plastic stress condition in a region D is that 
the system of partial differential equations in u,; and “2: must be 
hyperbolic or parabolic in D. 

If in the region D the system of equations is elliptical, then an 
elastic stress condition exists in D. 

As an example of the application of the criterion to a simple 
stress system; we restrict relationship (7) to only linear terms in 
e and obtain a system of equations as follows. The first equation 

Oke 


is: 
re) Ou; Ou; 
Ox; (2 Oxe2 Ox OX2 


1d [ /du,\? du \? 
+ 2 OX [(2") + ] ) 
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The second equation may be obtained by interchanging x; and xe in 


O72 
02x 0x2" 


072; 


O71) 
Ox) 0X2 


x2?” 


O72: 2 


the condition that the system be hyperbolic or parabolic is that 


the determinant equation: 

A, B G D £E F, 

A» Py Cy FE, F, 

ae dx, dx, 0 0 0 0 

= 0 
0 dx» dx; 0 0 0 | 0 (I ) 
0 0 0 dx, dx. O | 


dx, | 

has at least one real root (A; is the sum of all terms with p in 

Eq. (8), B; all terms with 7, etc.), giving us: 
A, ay (1 + S) + (1/2) ap (1 + t) 
B, = (1/2) + S) + doo lt 
= + (oe (1 + etc. 


(11) 


Table for Thermal Stresses in Rings* 


Marvin Forray 
Design Specialist, Structures, Applied Research and Development, 
Republic Aviation Corporation, Farmingdale, L.I., N.Y. 


February 2, 1960 


_ PLANE-STRESS, linear, elastic solution of a traction-free 
circular ring of inner radius a and outer radius } subjected 
to a general two-dimensional temperature distribution was con- 
sidered in references 1 and 2. Theory and formulas were de- 
veloped for the determination of the polar-coordinate stress com- 
ponents corresponding to temperatures of the form T = To(r/b)* 
cos or T = 7,(r/b)* sin where K and n are restricted 
to nonnegative integers. The purpose of this note is to give a 
tabular set of quantities for the direct calculations of the stress 
components corresponding to values of K and m from 0 through 
3 for a/b = 0.2, 0.5, and 0.8, respectively. This range for the 
parameters K and n should be adequate for most practical appli- 
cations. A numerical example is given for a particular tempera- 
ture distribution to demonstrate the application of the table. 

The stress components corresponding to a temperature dis- 
tribution of the form 


T = Ti(r/b)* cos (1) 
may be expressed in the form 
orr/EaTy = Bx, n cos 0 
= sin nO (2) 
ovo/EaTy) = Dx, n cos nd 


while corresponding to 
T = T)(r/b)* sin (3) 
the stress components are 
Orr/EaT) = Br, nsin 
o7/EaTy = —Ck,n cos (4) 


The quantities Bx,n, Cx,n, and Dx,, are functions of a/b, 
r/b, K, n and may be found from Table 1 (see p. 479). 

For example, given 
T = T)(r/b)*1 — cos 0) + = To(r/b)? — To(r/b)? cos 0+ Ty 


and a/b = 0.5, then, 


* This research was supported by the USAF under Contract No, AF33- 
(616)-6066, monitored by Mr. C. Richard, WCLSS-T30, WADC. 
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the first equation. If we denote by 


Oxo?’ Ox’ 


Grr/EaT, = Bo,» — Bz, cos | 
o9/EaTy) = Co,o — C2,1 sin @ ? (5) 

The nvmerically greatest stress component, for example, is 

700 | 


= |—— = |—0.375 — 0.1733) = 0.55 
EaTo+/b=1, 0=" 


|_760 

To max 
This is a compressive stress in the transverse direction at the 
outer boundary. 
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A Method of Tra nsforming Concentrated 
Surface Forces Into Continuous Surface Forces* 


William Zuk and Mohammed Abdul Majeed 
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[ SUCH FIELDS of applied mechanics as strength of materials 
and elasticity, the presence of concentrated forces often 
creates some difficulties as such forces represent mathematical 
discontinuities. A number of existing methods employing 
integration by segments,! Fourier series,? Fourier integrals,* unit 
impulse functions,* and Dirac delta functions’ may be used. To 
this list of known methods, another method is suggested by the 
authors, which is less mathematically complex than some of the 
other methods and which yields a single closed-form solution 
valid for all regions. To be able to express a function more 
compactly is desirable, not only for general simplicity but also 
for possible computer use wherein less storage capacity is needed. 

As an elementary example of this method, consider a simply 
supported beam with a single concentrated load as shown in Fig. 
1. Fig. 2 shows the Cauchy distribution curve’ whose area 
(representing the concentrated force) equals P as s > 0. (s is 
the standard deviation. ) 

Based on the controlling differential equations, 


d‘y/dx* = p(x)/EI 

d*y/dx* = —V/EI 

d*y/dx? = —M/EI 
analytical expressions for shear (V), moment (7), and deflec- 
tions (y) as a function of x may be had by elementary integration: 
Thus, 

d‘y (x) P 1 
= 
sm 1+ [(x — a)?/s] 


Upon integrating, 

Vizy = —EI(d*y/dx*) = —(P/r) tan~ [(x — a)/s] + G 
Integrating again, 
Miz) = —EI(d*y/dx*) = 


—(P/x)((x — a) tan—[(x — a)/s] — (s/2) In 
{1 + a)?/s?}}) + 


* This research was sponsored by the Virginia Council of Highway 
Investigation and Research. 
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(9) 
“cc 
Table for Thermal Stresses in Rings 
(5) TABLE 1 
a r 
0.2 0 O 0.2 0 3 6 -.3173 +1428 -.2576 0.5 0 3 -.0370k =-.03897 -.1009 
10 .2 0 01956 -.1678 +9 002267 202612 
1210 0 +2680 8 -21916 -.03935 1.0 +2364 
ok 9 206326 -.1260 1993 
5 #1333 205556 1.0 ° 25399 13 5 +2437 
6 21136 --02L70 ~=,001371 204046 
the 06042 -01715 =07806 1714 -0463b 8 20243 -.0701 
203072 002418 oh --1593 208596 605862 9 -.01527 01326 
1.0 ° 5 -.1722 -.1139 1.0 1649 
6 -.1520 01361 
20 .2 ° 1800 7 -.1166 -.09160 -.1160 23 5 1050 
3 +3036 8 -.07572 -.0993k 041 7h 6 201647 ~ 002650 
ok, 1571 22025 9 203534 09830 27 =-0089L,3 -00375 - 03808 
ace 5 21575 1.0 ° 315k 8 -.01090 =.03453 
6 1422 01778 29 01422 200423, 
e171 ° 23 2 ° 1730 1.0 08225 
8 208438 3 002315 202425 
1.0 5 05869 008976 003619 
6 - -054,86 -.01460 -205208 0.86 0 ° 
30 0 0 4085 7 = 204437 -.0319k -.05103 10 . ° 1037 
3 1110 22785 6 03852 002464 79 204922 
01438 22087 9 -.01478 = 202952 03559 +9 005533 ° -.001829 
5 1502 1.0 ° 21362 +95 6003862 0 -. 
6 1421 205843 1.0 ° ° =.09632 
7 1228 0 204926 33 ° 0 
8 209295 01884 20 ° +1800 
205226 23647 0.5 0 0 ° ° ° 285 007923 ° 208960 
1.0 0 0 5835 10 ° ° 22778 0099 0 
6 03457 295 6007091 0 08960 
ci 7 204218 ° 203560 1.0 ° ° -.1800 
3 ° 6 -035L2 ° 205765 
ol 01675 23026 9 22030 - 30 ° ° 22350 
5 el! -21870 1.0 ° ° 02222 285 122k 
7 209527 202262 20 5 ° ° 375 2009779 0 -.1201 
8 - -12k0 6 04,890 ° 2161 1.0 ° = 22530 
9 = 203233 222h6 06245 207255 
1.0 ° 8 205485 206985 o1 ° ° -.1183 
9 203285 -02179 -.005002 005002 = 205823 
11 ° ° 1.0 ° 0 +9 =.006060 - 006060 =002070 
als 21 .2 21202 295 =.00L159 -.004159 05113 
en 3 203258 203258 21177 30 23916 1.0 21021 
ok 21163 6 205310 ° «2476 
cal 04,808 +1024 07069 +1038 21862 0 
6 204696 207364 8 206479 ~-06013 21 6 208883 
ing 7 02929 204051 02529 003901 -003901 204795 
8 203170 203170 203088 1.0 004909 00499! 
nit 9 01788 +01 788 -.1070 295 003496 2003496 Oly 
T 1.0 ° -21990 5 ° --3778 1.0 208902 
6 004379 04379 22095 
he ° 01255 7 ~.05077 205077 -.07812 31 8 0 ° 
3 203608 203608 1372 8 03660 -85 .006912 2006912 
he ol 205210 205210 «1510 9 202285 -.02285 6008863 2008863 201199 
5 206058 206058 #1483 1.0 295 2006432 = 207229 
on 6 7 1.0 ° 21668 
7 205893 205893 11 0 ° 
ore 8 204770 204,770 202555 21 ° ° 02 ° 201601 
Is 9 202838 202838 202114 202114 285 0002976 -0006951 =,002729 
so 1.0 0 +3323 02720 02720 0001211 -. ~-008135 
8 C2432 202432 -.002980 295 =.0003678  =.0006623 -.001L57 
= 2 ° 1.147 9 201486 1.0 ° 201507 
<3 2008964, 203603 1.0 0 -01733 
ol - 209150 - 207704 -.1833 12 6 ° 0 201185 
ig 5 21240 -202771 ° ° 22250 285 0002002 -000L671  =.001722 
6 -01193 - 208500 --1955 6 03043 03043 +1736  =.00007815 -. -.005L8L, 
ea -.09717 -210L5 -.1162 7 204104, 204104, +1060 295 -.000LL75 -.001062 
8 -.06718 = 209373 20021: 8 203839 203839 2009609 1.0 ° 01037 
1S 9 -.03402 -.05783 +1550 9 -.1236 
=201529 
2 02 5 0 285 -.0002700 -.0006277 002024, 
3 2005426 6 ~006185 202131 201315 200005473 
ok 203032 20359 205552 2005336 2002346 -.06916 295 «0003384, 0006083 2001742 
5 ~.00375 08584 8 -.01323 -.05020 1.0 -.01463 
6 04,620 -.029L2 208516 -.008798 -.01470 201487 
7 -203978 -204118 205908 1.0 ° 1137 04019 
8 202892 203986 --009200 000751 002315 
201534 -.02600 20639 12 +1000 29 -.0003015 
1.0 ° 600 6 2003102 01150 95 31 =,.002207 003284, 
2C- 7 2002649 -001760 1.0 ° 203363 
2 2 .2-2 ° ° 8 -.006015 =.007015 =.02932 
n +232 ° ° 21630 23 8 0 
3 = 204122 202839 1.0 ° 206665 .0003807 0018LL - 
201288 -.02L37 01622 29 =,0002134 
5 202210 2004188 d 22 ° ° 295 =.0007645  -.001767 6002889 
6 202526 01205 205408 32 -.08135 1.0 0273 
202405 202264, 04828 6 -003166 -.01038 ~.00007708 
8 201927 -0258L, 201947 7 001995 -.002L5 203205 0 ° 201722 
201123 01894, 203756 8 200551 200606, 03046 285 0002147 2001034 =.00239h 
1.0 21285 2004928 ~008166 -.003125 «9  ,0002362 -.00010,8  -.008109 
1.0 ° -.071L3 095 =.0004319  -.0009965 -.001777 
© 3 .2 ° ° 1.633 1.0 ° ° 201576 
3 22061 3910 202, 3 24061 
oh - 21926 6 -.002817 206325 ~.02870 33 ° 
5 -.03169 ~.2589 2007097 
C2 
ay 


‘the 
— 


P= p(x) 


Fic. 1. Simply supported beam. 


P(x) 


— 


L 


Fic. 2. Force-distribution curve. 


Since s is taken approaching zero to represent the concentrated 
force, terms multiplied by s vanish. Therefore, 


= —(P/r)(x — a) tan [(x — a)/s] + + 


The constants C; and C; are determined by the two boundary 
conditions 


(d*y/dx*)(x = 0) = 


and 
(d?y/dx?)(x = 1) = 0 
to be 
Ci = (P/2) — (Pa/L) and C, = Pa/2 
Thus, 


Viz) = —(P/r) tan [(x — a)/s] + (P/2) — (Pa/L) (1) 
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and 
Mw) = —(P/#Xx — 4) tan= [(x — «)/s] + 
(P/2L)(L — 2a) x + (Pa/2) (2) 

For 0 < x < a, the expression tan~! {[(x — a)/s] as s > O equals 
—(2/2),and +(2/2)fora<x<L. 

To obtain the deflection y(x), Eq. (2) is integrated twice. 
Omitting terms multiplied by s, 
El(dy/dx) = (P/2r)(x — a)? tan [(x — a)/s] — 
(Pa/2) x — (P/2L)(L — 2a)(x?/2) + C; 


and 


Ely = (P/6x)(x — a)? [(x — a)/s] — 
(Px?/12L)(L — 2a) — (Pa/4)x? + Cyxx + C, 


The constants C; and C, are evaluated from the boundary condi- 
tions 


= Oand yr=1) = O 


to be 


C3; = (Pa/12L)(4L? — 3aL + 2a?) 


and 


C, = —(Pa*/12) 


Thus, 
y(x) = (P/6xEI)(x — a)* tan [(x — a)/s] — 

(Px*/12EIL)\(L — 2a) — (Pa/4EI)x? + 
(Pa/12EIL)(4L? — 3aL + 2a?)x — (Pa*/12EI) (38) 


Eqs. (1)-(3) may easily be verified to give the proper values 
for shear, moment, and deflection for any value of x along the 
beam as s — 0. 

Although the example illustrated is a simple one perhaps more 
easily solved by another method, it is believed that use of the 
Cauchy curve for concentrated forces in general will yield easier 
closed-form solutions for many categories of problems in applied 
mechanics than methods heretofore used. 


REFERENCES 


i Timoshenko, Strength of Materials, Part 1, D. Van Nostrand Co., 1955. 

2 Timoshenko and Goodier, Theory of Elasticity, McGraw-Hill Book Co., 
Inc., 1951. 

3 von Karman, Th., and Biot, M. A., Mathematical Methods in Engineer- 
ing, McGraw-Hill Book Co., Inc., 1940. 

4Carslow and Jaeger, Operational Methods in Applied Mathematics, 
Oxford University Press, 1941. 

6 Hald, Statistical Theory with Engineering Applications, John Wiley, 1952. 


43 Schott, G. A., Electromagnetic Radiation, Cambridge Uni- 
versity Press, 1912. 

44 Rott, Nikolaus, Das Feld einer raschbewegten Schallquelle, 
A. G. Gebr. Leemann & Co., Zurich, 1945. 

46 Evvard, J. C., Use of Source Distributions for Evaluating 
Theoretical Aerodynamics of Thin Finite Wings at Supersonic 
Speeds, NACA Rpt. 951, 1950. 

46 Krasilshchikova, E. A., From Scientific Records of the Mos- 
cow State University, Vol. 154, Mechanics No. 4, 1951, with 
appendix condensed from a document Modern Problems of Me- 
chanics, Govt. Pub. House of Tech. Theor. Literature, Moscow, 


Theory of Supersonic-Propeller Aerodynamics 
(Continued from page 450) 


Leningrad, 1952. Translation in Finite Span Wings in Compres- 
sible Flow, NACA TN 1383, 1956. 

47 Hale, R. W., Aerodynamic Tip Theory of a Supersonic Pro- 
peller, M.S. Thesis, Graduate School of Aeronautical Engineering, 
Cornell, February, 1958, and OSR M.A. TN-58-335, ASTIA No. 
AD 154239, February, 1958. 

48 Sears, W. R., Ed., General Theory of High Speed Aero- 
dynamics, Vol. VI, Princeton University Press, 1954. 

49 Cohen, Clarence B., Influence of Leading-Edge Suction on 
Lift-Drag Ratios of Wings at Supersonic Speeds, NACA TN 1718, 
1948. 


480 j 
a 
~ 
: 
pix) s 
A 
a L-@ 
[ 
7 
: 
; 
t 
ay: 


1955. 
Cow 


ineer~ 


atics, 


Man 
desired 
script 
be dou 
consect 
for the 
ing, or 
therefo 
Typogt 
wish to 
mailing 
fat. 
protect 
Science 
will be 

Author 
is acce] 

and ini 
use of t 
quent d 
it diffic 
ticularl: 
services 
of the 


that pre 
REFE!I 


REQUIREMENTS for CONTRIBUTIONS 
to the publications of the 
INSTITUTE of the AERONAUTICAL SCIENCES 


The Institute of the Aeronautical Sciences invites both members and nonmembers from any 
country to submit papers for publication in the JoURNAL OF THE AERO/SPACE SCIENCES and 
AERO/SPACE ENGINEERING. The Institute, following the practice of other societies, does not pay for 


contributions. 


The following directions for the preparation of papers, if followed by authors, will save corre- 
spondence, avoid the return of papers for changes, minimize the work of preparation for the printer, 
and save the expense due to the charges made for “‘author’s corrections.” 


Manuscripts: The original typewritten copy of the paper is 
desired. To expedite review, an additional copy of both the manu- 
script and figures should be submitted. The manuscript must 
be double or triple spaced on one side of white paper sheets, 
consecutively numbered. There should be wide margins to allow 
for the marking of directions to the printer. Correcting, chang- 
ing, or adding to papers after they are in type is costly. It is, 
therefore, imperative that papers submitted be in final form. 
Typographical errors may be corrected on proofs, but if authors 

mwish to add material, they may do so at their own expense. In 
mailing, drawings may be rolled, but manuscripts should be sent 
flat. Send by first-class mail (register if you wish for your own 
protection) to the Editorial Office, Institute of the Aeronautical 
Sciences, 2 East 64th Street, New York 21, N.Y. All manuscripts 
will be examined by the Editorial Committee and by the Editor. 
s Authors will be advised as promptly as possible whether the paper 
is acceptable for publication. 

TrrLEs: The title of the paper should be brief. 


quent duplication of initials and surnames which sometimes makes 
it difficult to establish the identity of the author. This is par- 


of the organization with which the author is associated should 
be placed after his name on a separate line. The date on which 
the paper is received will be inserted by the Editor. The author's 
title or position should be indicated in a footnote. 
® SUMMARIES OR AsstrRacts: An abstract to be printed at the 
beginning should accompany each article. It should be adequate 
as an index and asa summary. It should contain a statement of 
Mmajor conclusions reached, since summaries in many cases con- 
itute the only source of information used in compiling scien- 
tific reference indexes. Abstracts printed in other journals, espe- 
jally foreign, in most cases consist of summaries from printed 
E The summary should explain as adequately as possible 
he major conclusions to a nonspecialist in the subject and should 
ontain from 100 to 300 words, depending on the length of the 


SUBHEADINGS: Subheadings should be inserted by the author 
fat frequent intervals. The work of editorial preparation will be 
simplified by the author’s provision of many subheadings. 

MattTer UsuaLty DELETED: Photographs or illustrations of 
little technical interest and not showing advances in gen 
practice. Too detailed tabular matter (general results of such 
tables may be included in the text). Lengthy descriptions of 
materials or processes or of preliminary experiments or theories 
that preceded final results; salient features only are of interest. 

REFERENCES: References should be numbered consecutively 
and grouped together at the end of the manuscript, with only 


aftangement should be as follows. 

Aerodynamic Theory, Ist Ed., Vol. 1, p. 23; Julius Springer, 
Berlin, 1934. For magazines: 'England, . R., Crawford, 
A, B., and Mumford, W W., Some Results of a Study of Ultra- 
Short-Wave Transmission Phenomenon, Proc. IRE, Vol. 20, 
No, 12, pp. 481-482, March, 1933. Please give author, title, 
dition, volume, number, page, publisher, and place and date of 
publication as indicated Omission of one required fact causes 
much extra editorial work and possible inaccuracies. 


ILLUSTRATIONS: I!lustrations should accompany manuscripts, 
and each should always be ref to in the text by number. 
Drawings or graphs should not be larger than 12 X 16 inches, 
and must be made with jet black India ink on white paper or 
tracing cloth, the latter being preferred. Do not use typewriter 
for lettering. The smallest lettering on 8 X 10 inch figures should 
be no less than !/, inch high. Cross-section paper (white with 
black lines) may be used, but it should not have more than 4 lines 
per inch. If finer ruled paper is used, the major division lines 
should be drawn in with black ink, omitting the finer divisions. 
In the case of finely ruled paper, only blue-lined paper can be 
accepted. Tracing paper and blueprints are not acceptable. 
Lettering and all markings must be large enough to be readable 
after reduction to single-column width (3°/\. in.).. Mail rolled or 
flat; never fold. Drawings that cannot be reproduced (including 
pencil drawings) will be returned to the author for redrawing, 
thus delaying publication of the paper. Photographs should be 
distinct and show clear black and white contrasts. They must 
be on glossy white paper. Avoid round or oval photographs. 


CAPTIONS AND LEGENDS: Legends or captions must accompany 
each drawing or photograph submitted. If written on the draw- 
ing or photograph they should be placed below and well outside 
the part to be reproduced. Each table should have a caption 
such as Table 1, Table 2, Table. 3, etc. Captions should be com- 
plete in themselves so as to make the data intelligible to the 
reader without reference to the text. A duplicate list of captions 
for figures should be included as the last page of the manuscript. 
Use ‘‘Fig. 1’’ (not Figure 1), ‘Figs. 3 and 4,” etc., in both the text 
and the numbering of illustrations. In the text, “Eq. (1)” or 
“Eqs. (1) and (2)” is used, not “Equation (1).” In captions, 
legends, and in table headings, write all words in full; do not 
abbreviate, except for “Fig.” and “Eq.” 


MATHEMATICAL WorRK: Formulas may be typewritten or 
carefully written in pen and ink, the writing to be large enough so 
that it can be marked for the printer. Considerable space for mark- 
ing should be allowed above and below all equations. All compli- 
cated equations should be repeated on separate sheets with plenty 
of space left for marking. The solidus should be used for sim- 
ple fractions appearing within the text. Make all expressions 
clear to the typesetter. Greek letters used in formulas should 
be clearly designated by name in the margin of the manu- 
script. The difference between capital and lower-case letters 
should be clearly distinguished and care taken to avoid confu- 
sion between zero (0) and the letter (0), between the numeral 
(one) and the letter (ell) and the prime (’), between alpha and 
a, kappa and k, u and mu, v and nu, n and eta. All sub- 
scripts and exponents should be clearly distinguished. Avoid 
complicated exponents and subscripts. Dots and bars over 
letters or mathematical expressions should also be avoided. 
When it is necessary to repeat a complicated expression, it should 
be represented by some convenient symbol. 


SYMBOLS AND ABBREVIATIONS: The symbols recommended 
in the American Standard Association, “Letter Symbols 
for Aeronautical Sciences,” ASA Y10.7—1954, should be used 
wherever practicable. All symbols should be clearly written and 
carefully checked. Standard abbreviations should be used, and 
it should be noted that most abbreviations are lower case, such 
as m.p.h., b.m.e.p., ih.p., b.hp., hp., ete. 


~ 
¥ 
f 
; 
4 
4 
ig land initials of the author should be written as he prefers. The 7 
fee use of the full name of an author is advocated because of the fre- bl 
fm meservices. All titles and degrees or honors are omitted. The name é 
of 
| 
| ‘ 
1€ a 
4 


CorRPORATE MEMBERS 


INSTITUTE OF THE AERONAUTICAL SCIENCES ae 


AC SPARK PLUG DIVISION, GENERAL MOTORS CORPO- 
RATION 


ACADEMY OF AERONAUTICS 
AEROJET-GENERAL CORPORATION 

AEROLAB DEVELOPMENT COMPANY 

AERONCA MANUFACTURING CORPORATION _ 
AERONUTRONIC SYSTEMS, INC. 

AEROQUIP CORPORATION 

AGAWAM AIRCRAFT PRODUCTS, INC, 

ALLIED RESEARCH ASSOCIATES, INC. 

ALLISON DIVISION, GENERAL MOTORS CORPORATION 
ALUMINUM COMPANY OF AMERICA 

AMERICAN AIRLINES, INC. 

AMERICAN BOSCH ARMA CORPORATION 


AMERICAN STEEL & WIRE DIVISION, UNITED STATES 
STEEL. CORPORATION 


AMOCO CHEMICALS CORPORATION 
AMPHENOL-BORG ELECTRONICS CORPORATION 
ANDERSON, GREENWOOD & CO, 

ARDE. ASSOCIATES 

ASSOCIATED AVIATION UNDERWRITERS 
AUTOMATION INDUSTRIES, INCORPORATED 
AVCO RESEARCH LABORATORY 

AVIEN, INC. 

BEECH AIRCRAFT CORPORATION 

BELL AIRCRAFT CORPORATION 

BENDIX AVIATION CORPORATION 

BOEING AIRPLANE COMPANY 

BOOZ, ALLEN & HAMILTON 


BULOVA. RESEARCH, AND. DEVELOPMENT LABORA- 
TORIES, INC, 


CANADAIR, LTD. 
THE CESSNA AIRCRAFT COMPANY 

CHANCE VOUGHT AIRCRAFT, INCORPORATED 
THE CHASE MANHATTAN BANK 

CHICAGO AERIAL INDUSTRIES, INC. 

THE CLEVELAND PNEUMATIC INDUSTRIES, INC. 
CONTINENTAL MOTORS CORPORATION 
CORNELL ‘AERONAUTICAL LABORATORY, INC. 
CURTISS-WRIGHT CORPORATION 

DANIEL, MANN, JOHNSON, & MENDENHALL 
THE DECKER CORPORATION 

DEL MAR ENGINEERING LABORATORIES 


OF THE 


Ti 
DOAK AIRCRAFT COMPANY, INC. 
DOUGLAS AIRCRAFT COMPANY, INC. 
DZUS FASTENER COMPANY, INC. 
EASTERN AIR LINES, INC. 
EATON MANUFACTURING COMPANY 
EDO CORPORATION 
ELASTIC STOP NUT CORPORATION OF AMERICA 
ENGINEERING SUPERVISION COMPANY 
ee CAMERA AND INSTRUMENT CORPORA- 


DIVISION, GENERAL MOTORS CORPORA- 


FAIRCHILD ENGINE AND AIRPLANE CORPORATION 
THE GARRETT CORPORATION 

GENERAL APPLIED SCIENCE LABORATORIES, INC. 
GENERAL DYNAMICS CORPORATION 

GENERAL ELBUTRIC COMPANY 

GENERAL PRECISION EQUIPMENT CORPORATION 
G. M. GIANNINI & CO., INC. 

GIANNINI PLASMADYNE CORPORATION 

THE B. F. GOODRICH COMPANY 

GOODYEAR AIRCRAFT CORPORATION 

GRUMMAN AIRCRAFT ENGINEERING CORPORATION 


HYDRO-AIRE CO. 


INSURANCE COMPANY OF NORTH AMERICA COM- 
PANIES 


INTERNATIONAL BUSINESS MACHINES CORPORATION 
THE INTERNATIONAL NICKEL COMPANY, INC, 


TTT LABORATORIES, DIVISION OF INTERNATIONAL 
TELEPHONE AND TELEGRAPH CORPORATION 


JANITROL AIRCRAFT, A DIVISION OF MIDLAND- 
ROSS CORPORATION 


JOHNS-MANVILLE SALES CORPORATION 
WALTER KIDDE & COMPANY, INC. 

KOLLSMAN INSTRUMENT CORPORATION 

LAVELLE AIRCRAFT CORPORATION 

LEAR, INCORPORATED 

C. W. LEMMERMAN, INC. 

LIBRASCOPE DIVISION OF GENFRAL PRECISION, INC, 
THE LIQUIDOMETER CORPORATION 

LOCKHEED AIRCRAFT CORPORATION 


LOEWY-HYDROPRESS DIVISION OF BALDWIN-LIMA- 
HAMILTON CORPORATION 


THE MARQUARDT CORPORATION 

THE MARTIN COMPANY 

McDONNELL AIRCRAFT CORPORATION 

MELETRON CORPORATION 
MINNEAPOLIS-HONEYWELL REGULATOR COMPANY 
NATIONAL CREDIT OFFICE, INC, 

NORTH AMERICAN AVIATION, INC, 

NORTHROP CORPORATION 


NUCLEAR DEVELOPMENT CORPORATION OF AME 

PAN AMERICAN WORLD AIRWAYS, INC. 4 
THE RALPH M, PARSONS COMPANY a 


PIASECKI AIRCRAFT CORPORATION 
RADIO CORPORATION OF AMERICA 
ASTRO-ELECTRONIC PRODUCTS DIVISION 
DEFENSE ELECTRONIC PRODUCTS ae 


RAMO-WOOLDRIDGE DIVISION, THOMPSON R 
WOOLDR4DGE INC, 


REPUBLIC AVIATION CORPORATION 
ROHR AIRCRAFT CORPORA‘ZION 
PAUL ROSENBERG ASSOCIATES 
RYAN AERONAUTICAL COMPANY Ey 
SANDBERG-SERRELL CORPORATION ae 
SHAFER BEARING DIVISION, CHAIN BELT COMPAMl 
SHELL OIL COMPANY 
SIMMONDS PRECISION PRODUCTS, INC. 
SOCONY MOBIL OIL COMPANY, INC, 
SOLAR AIRCRAFT COMPANY 
SOUTHWEST PRODUCTS CO. 

SPACE TECHNOLOGY LABORATORIES, INC. q 
R. DIXON SPEAS ASSOCIATES ; 


SPERRY GYROSCOPE COMPANY DIVISION OF SPER 
RAND CORPORATION 


STANDARD OIL COMPANY OF CALIFORNIA 
STANDARD OIL COMPANY (INDIANA) 


STANDARD-THOMSON CORPORATION 
CLIFFORD MANUFACTURING DIVISION 


STANLEY AVIATION CORPORATION 
THIOKOL CHEMICAL CORPORATION a 
THOMPSON RAMO WOOLDRIDGE INC. 
TRANS WORLD AIRLINES, INC. 
TURBO PRODUCTS, INC. 

UNION CARBIDE CORPORATION 
UNITED AIR LINES, INC. 

UNITED AIRCRAFT CORPORATION 
UNITED STATES AVIATION UNDERWRITERS, INC, 
VERTOL AIRCRAFT CORPORATION 
VITRO CORPORATION OF AMERICA 
WESTERN GEAR CORPORATION 
WESTINGHOUSE ELECTRIC CORPORATION 
WYMAN.GORDON COMPANY 

YOUNG RADIATOR COMPANY 


PHILLIPS PETROLEUM COMPANY 4 


LS 
f a 
4 
q 
4 
4 
ae 
ay | 
Bs 
4 
4 
i 
Be 
. a 
& 
% 
q 
~ i 


